Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,806	G→A	P152P (CCG→CCA)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,806	0	G	A	96.2%	68.9 / NA	26	P152P (CCG→CCA)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/10);  minor base T (0/1);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCC‑AGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGT  >  NC_000913/71560‑72048                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      cGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGt                                                                                                                                                                                                                                                   >  1:66891/1‑251 (MQ=255)                                tttAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGgcgc                                                                                                                                                                                                                   >  1:144728/1‑251 (MQ=255)                                   aaaaaGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCCGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCgcggg                                                                                                                                                                                                                >  2:182661/1‑250 (MQ=255)                                             ttcATCGCTGAGCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGATCATGCGTTGCATCAACACATCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGatat                                                                                                                                                                                                       <  2:191078/248‑1 (MQ=255)                                                         cATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAAt                                                                                                                                                                                           <  1:41966/250‑1 (MQ=255)                                                               tGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCagagg                                                                                                                                                                                    >  1:85772/1‑250 (MQ=255)                                                                aaagaaAAACAAATCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATTTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCAGGTGCCAATGG‑GGTGGGAAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATAGTATCGACTGCCTGCTGTTTCCACCGTTTTTCGTCCTGCC‑TGGGATTCTGGCGGGCGCTGCGGTCGATATTCCAGTGGGAATGCAGAGCg                                                                                                                                                                                   <  2:245737/250‑1 (MQ=255)                                                                 aagaaAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTAAGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCgg                                                                                                                                                                                  >  1:166446/1‑251 (MQ=255)                                                                      aaaCAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGt                                                                                                                                                                             >  1:31144/1‑251 (MQ=255)                                                                                        tAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCC‑GACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTTCAGCCGGAATGTCGAGCGGCGAGGTTAAATGGTTTCTGCTggg                                                                                                                                                          >  1:196220/1‑250 (MQ=255)                                                                                             cTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGGAGAGGGGCGAGTTTAAATGGTTGCTGCTGGcgacg                                                                                                                                                      >  1:198240/1‑247 (MQ=255)                                                                                                                   caGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGAcggcgggatttgtggggggtggtg                                                                                                                                >  1:79230/1‑233 (MQ=255)                                                                                                                    aGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTAAGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGGTGCTGCTGGCGACGGCGGGATTTGGGT‑GGGTtggtgg                                                                                                                              >  1:70463/1‑251 (MQ=255)                                                                                                                       aTGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGGTCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGGTTTTTTGT‑GGGTTTGtggctg                                                                                                                           >  1:197781/1‑251 (MQ=255)                                                                                                                                    tGGTCGGTCGCTTTGTTGTCCCGACGCGTCCGCTGGTGCCAATGGCGGCGGGCA‑GTATGGATCTGCCGGTGGCAAAATTCATTTCGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑CGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGGGTTTAAATGGTTGCTGCTGGCGACGGCGG‑GATTTGGTTGGGTTGGTgggtggcttttgtggcgg                                                                                                              >  1:268400/1‑241 (MQ=255)                                                                                                                                      gtcggtcgCTTTGTTAGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTTCC‑TGGGATTCTGGCGGGAGCGGCGATCGATATTCCCGCCGGAATGCAGAGCGGCGAGGTTAAAAGGTTGCTGCTGGCGAAGGCGGTATTTTTGT‑GGGGTGGGGGCTGGCTTTGCTggggggt                                                                                                            >  1:217854/1‑246 (MQ=255)                                                                                                                                                  gttgGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCgg                                                                                                <  1:47816/251‑1 (MQ=255)                                                                                                                                                             gggaccGCTGGGGCCCATTGTGGCGGTGA‑TGCTGGAGATGGCGGTGGCGAAATTCATTACGCCGAATATTCTGGCCTGCCTGCGGTGGCCTCTGTTTTGCTCCCTGCCGAGGG‑TGCTGGCGGGCGCGGTGATCGGTATTCCAGCGGGAATGCAGAGCGGCGAGTGTAAATGGTTGCTGCTGGGGACGGCGGTGTGTGTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaaccg                                                                                     <  2:61791/247‑6 (MQ=255)                                                                                                                                                                 ccGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaaccgatc                                                                                  <  1:229292/250‑9 (MQ=255)                                                                                                                                                                                  ggCGGGAA‑TGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTg                                                                 <  1:15759/250‑1 (MQ=255)                                                                                                                                                                                        aa‑TGCTGGATCTGCCGGTGGCAAAAGTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCGGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTcgcgcg                                                          <  2:240090/251‑1 (MQ=255)                                                                                                                                                                                                     gccggtgatacaaTTCATTACGCCGAATATTATCGGCTGCCGGCTGTGGCCACCGGTTTACTTCCTGCC‑GGGGATACTGGCGGGCGCGGCGATCGATATTCCAGGCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTGTGGCGCAGCGGTAAAGCACCCGATCGTTTAAGTCATTGTTTGTCGCGCGGCCG‑Attgtt                                               <  2:268400/238‑1 (MQ=255)                                                                                                                                                                                                               aaaTTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTTTGGCGCAGCGGGAAAGCAACCCATCGTTTTAGTCATTATTTTTTGCGCGGGCG‑TTTTTTGTTGCttgcgcc                                   >  1:217954/1‑245 (MQ=255)                                                                                                                                                                                                                          gCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGGGCAGCGGGAAAGGAACCGATCGTTTTAGGCATTATTTTTCTCGCGGCCGGTTTGT‑GTGGGTGACGCCggtgttttttg                        >  1:196775/1‑241 (MQ=255)                                                                                                                                                                                                                                           gcctgcTGTGGCCACCGTTTTTCTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGgc        <  2:211933/250‑1 (MQ=255)                                                                                                                                                                                                                                                  gtgGCCACCGTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCtggt  <  2:144728/249‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      CGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCC‑AGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGT  >  NC_000913/71560‑72048

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