Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	134,256	G→T	intergenic (+44/‑132)	acnB → / → yacL	AcnB/PF06062 family protein YacL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	134,256	0	G	T	96.3%	89.0 / NA	27	intergenic (+44/‑132)	acnB/yacL	AcnB/PF06062 family protein YacL
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (12/14);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGCGCGAATGTCTTCCTGGCTT‑CTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGTAAAGGTGCGTATGTCCATGGGGCATGAAGTGGTCGGGCACTGGTTTAATGAAGAGGTGAA  >  NC_000913/134012‑134485                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      ggCGCGAATGTCGGCCTGGCGG‑CTGCGGAGCTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCGGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGGGAGCTGAACTTCAACCAGCCTTCGCAGTGCACCGAGAAAGCCGATGGGGTGATTTGCCAGACTGCGGGTTAAAAAGGCAGCGGACGCGATGCGCATAAACGACACAATGCCCGGGTAATgag                                                                                                                                                                                                                                 <  1:127502/251‑1 (MQ=255)    gcgAATGTCTTCCTGGCTTGTTGCG‑AACTGGCGGCTGTTGCGGCGCGGATTGGCAAACTGGCGGCGCCTGAAGAGTTCCAGGCCTACGTGGCGCATGTAGAAAAAACAGCCGTTGATTCTTTCCCTTTTCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGGTGGGGTGATTTTCCCGACTGAGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGTTTAATGAGAtt                                                                                                                                                                                                                              <  2:2702/251‑1 (MQ=255)       aaTGTCTTCCTGGCTT‑CTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTccc                                                                                                                                                                                                                           <  1:174481/251‑1 (MQ=255)          gTCTTCCTGGCTT‑CTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCgg                                                                                                                                                                                                                         <  2:79707/250‑1 (MQ=255)              tCCTGGCTT‑CTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGAGTCCCGGGc                                                                                                                                                                                                                       <  2:18965/248‑1 (MQ=255)                                cTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAAcc                                                                                                                                                                                                    >  2:225906/1‑250 (MQ=255)                                                                   gACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCg                                                                                                                                                                <  1:9830/251‑1 (MQ=255)                                                                                  ccAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATg                                                                                                                                                 <  2:136767/251‑1 (MQ=255)                                                                                      aCCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGgcgc                                                                                                                                             >  2:106949/1‑251 (MQ=255)                                                                                              ggCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTGTTATTTCTAATCCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTTCGATAATTACATCAATGGCGCAATGCGAt                                                                                                                                     <  1:270200/251‑1 (MQ=255)                                                                                                   aGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCgg                                                                                                                                >  1:51711/1‑251 (MQ=255)                                                                                                    ggTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGt                                                                                                                               <  1:132525/251‑1 (MQ=255)                                                                                                               cAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGAAATGCGATTTCGGTGCATTGCCggg                                                                                                                    >  1:200326/1‑251 (MQ=255)                                                                                                                          aCTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAac                                                                                                         <  1:305754/251‑1 (MQ=255)                                                                                                                               ccGTTATCTGAACTTCAACCAGCTTTCTCAGTACCCCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTGAAAAAGTCAGCGCACGCGCTGCGCTTAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCAGCGCCGTTCCGCTGTTTTTCTCCGGTATGGTTTCGATAAATACATCAACGGCGCAATGCTATTACGCTGCCTCGCCGGGAGCAGAGGTacacacc                                                                                                    <  2:141391/251‑2 (MQ=255)                                                                                                                                     tCTGAACTTCAAGCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCTGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTCTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTACGGCTGCGATAATTACATCACTGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGAt                                                                                              >  2:309818/1‑251 (MQ=255)                                                                                                                                          aCTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACg                                                                                          >  1:250525/1‑250 (MQ=255)                                                                                                                                                                     agaAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCt                                                                                                                                                                   >  1:236084/1‑150 (MQ=255)                                                                                                                                                                     agaAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCt                                                                                                                                                                   <  2:236084/150‑1 (MQ=255)                                                                                                                                                                                     gATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGTAAAGGTGCGTATGTcc                                              >  2:266898/1‑251 (MQ=255)                                                                                                                                                                                          tCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGTAAAGGTGCGTATGTCCATggg                                         >  1:223897/1‑251 (MQ=255)                                                                                                                                                                                               aCTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGTAAAGGTGCGTATGTCCATGGGGCATg                                    >  1:86484/1‑251 (MQ=255)                                                                                                                                                                                                   cGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGTAAAGGTGCGTATGTCCATGGGGCATGAAGt                                <  1:95058/251‑1 (MQ=255)                                                                                                                                                                                                              gTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTGTTTATTTCTAAACCATCGCCGTTCCGCTGGTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGTAAAGGTGCGTATGTCCATGGGGCATGAAGTGGTCGGGCACt                     <  2:135797/251‑1 (MQ=255)                                                                                                                                                                                                                        gcgcTGCGCATAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGTAAAGGTGCGTATGTCCATGGGGCATGAAGTGGTCGGGCACTGGTTTAATg            >  1:95336/1‑250 (MQ=255)                                                                                                                                                                                                                                 aTAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGTAAAGGTGCTTATGTCCATGGGTCATGAAGTGGTCGGGCACTGGTTTAATGAAGAGGTgaa  <  1:290254/251‑1 (MQ=255)                                                                                                                                                                                                                                 aTAAACGACACAATGCCCGGTTAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGTAAAGGTGCGTATGTCCATGGGGCATGAAGTGGTCGGGCACTGGTTTAATGAAGAGGTga   >  1:208090/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      GGCGCGAATGTCTTCCTGGCTT‑CTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTCCCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCCGGTAAGGCTGCGATAATTACATCAATGGCGCAATGCGATTTCGGTGCATTGCCGGGAGCAGAGGAACACACTATGGATTACGAATTTCTGCGCGATATTACCGGAGTGGTAAAGGTGCGTATGTCCATGGGGCATGAAGTGGTCGGGCACTGGTTTAATGAAGAGGTGAA  >  NC_000913/134012‑134485

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