Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,776,162:1	+G	11.6%	intergenic (+3/+9)	mtlR → / ← yibT	mannitol operon repressor/uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,776,162	1	.	G	11.6%	152.6 / 18.4	43	intergenic (+3/+9)	mtlR/yibT	mannitol operon repressor/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (21/17);  new base G (3/2);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

CACCGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAA  >  NC_000913/3776014‑3776308                                                                                                                                                      |                                                                                                                                                   caccGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCa                                                                                                                                                     >  1:25548/1‑149 (MQ=255)         ccACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGc                                                                                                                                            >  4:9481/1‑149 (MQ=255)                ttGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGc                                                                                                                                     <  2:38298/149‑1 (MQ=255)                    aCCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCct                                                                                                                                 >  4:93661/1‑149 (MQ=255)                              tCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAGACGT‑TGCTTACttt                                                                                                                                     <  1:100754/135‑3 (MQ=255)                              tCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAGACGT‑TGCTTACttt                                                                                                                                     >  2:100754/1‑133 (MQ=255)                                   tttGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAGACGT‑TGCTTACtt                                                                                                                                      >  4:1773/1‑128 (MQ=255)                                   tttGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAGACGT‑TGCTTACtt                                                                                                                                      <  3:1773/129‑2 (MQ=255)                                     tGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGa                                                                                                                <  1:60435/149‑1 (MQ=255)                                          gcaatgcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAg                                                                                                           >  1:82617/1‑149 (MQ=255)                                                 aaaTTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCcgacga                                                                                                    <  1:8909/149‑1 (MQ=255)                                                         cgcTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTc                                                                                            <  1:92843/149‑1 (MQ=255)                                                                aaCAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGt                                                                                     >  3:8208/1‑149 (MQ=255)                                                                 aCAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATATTCCAGACGtt                                                                                    >  2:30694/1‑149 (MQ=255)                                                                    ggCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTa                                                                                 <  1:47860/149‑1 (MQ=255)                                                                         tgCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGatat                                                                            >  4:52121/1‑149 (MQ=255)                                                                            gATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTAc                                                                                                                             <  4:33634/97‑1 (MQ=255)                                                                            gATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTAc                                                                                                                             >  3:33634/1‑97 (MQ=255)                                                                             aTCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACg                                                                        >  2:75768/1‑149 (MQ=255)                                                                                aacaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCa                                                                     >  3:20729/1‑149 (MQ=255)                                                                                aacaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCa                                                                     >  3:84464/1‑149 (MQ=255)                                                                                 acaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCAc                                                                    <  2:23019/149‑1 (MQ=255)                                                                                   aaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGc                                                                  <  4:84464/149‑1 (MQ=255)                                                                                        tCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCa                                                             >  1:36704/1‑149 (MQ=255)                                                                                                   actgaGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTcc                                                  <  2:25548/149‑1 (MQ=255)                                                                                                       aGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCcggac                                              >  2:87290/1‑149 (MQ=255)                                                                                                        gCTGATTTTGGAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTc                                                   <  1:20787/143‑1 (MQ=255)                                                                                                          tGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCcggacgga                                           >  4:75642/1‑149 (MQ=255)                                                                                                              ttCCAAAGTCAGCTTAAAACAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCGGCTTCGGCCGATAAAGCCGCCGATAATACTCAAGACGTTGTAGATATAACGGAACGCTTTCATCGGGTATTCCGGACGGAATcg                                       >  2:98289/1‑149 (MQ=255)                                                                                                                  aaaaTCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGtt                                   <  1:61194/149‑1 (MQ=255)                                                                                                                         gCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGgga                            <  1:57750/149‑1 (MQ=255)                                                                                                                          cTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGTCCGGAAtaccc                                     <  2:60181/139‑5 (MQ=255)                                                                                                                          cTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGTCCGGAAtaccc                                     >  1:60181/1‑135 (MQ=255)                                                                                                                               acaaaGCCTTTCAAAAGTAAGCA‑ACGTCTCCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGtt                                   <  2:444/134‑1 (MQ=255)                                                                                                                               aaaaaGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGtt                                   >  1:444/1‑136 (MQ=255)                                                                                                                                  aaGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGtttt                   <  4:27027/149‑1 (MQ=255)                                                                                                                                  aaGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGtttt                   >  1:68596/1‑149 (MQ=255)                                                                                                                                    gCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGtttttt                 >  2:15003/1‑149 (MQ=255)                                                                                                                                    gCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGtttttt                 <  2:61547/149‑1 (MQ=255)                                                                                                                                     ccTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTAGCTGCTTCGGCCGATAAAGACGACGATAATAATCCAGACGTTGTAGATATAACGGCACGCTTTAATCGGGTATTTAGGACGGAATAGCGTTACGAGGAGGAAGTTTTTTc                >  2:88083/1‑149 (MQ=255)                                                                                                                                         tctaaaGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCGGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGAt            >  4:7879/4‑149 (MQ=255)                                                                                                                                            gcaGTAAGCAGACGT‑TGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCTGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACt         >  1:47453/3‑149 (MQ=255)                                                                                                                                                tAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCg     >  1:100679/1‑149 (MQ=255)                                                                                                                                                   gCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGaa  <  3:93661/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CACCGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCA‑ACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCGTTACGGGGAGGAAGTTTTTTCAGATACTCCCGGAA  >  NC_000913/3776014‑3776308

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑