Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,824	A→T	P158P (CCA→CCT)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,824	0	A	T	96.0%	75.0 / ‑6.6	27	P158P (CCA→CCT)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base T (13/11);  minor base C (1/0);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCC‑AGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTA  >  NC_000913/71591‑72055                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                          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)                                                                                                                                                                                                                       ccaccgTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGtt   <  2:85772/248‑1 (MQ=255)                                                                                                                                                                                                                        caccgTTTTACTTCCTGCC‑TGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTa  <  2:197781/249‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                          TTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAA‑TGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCC‑AGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGT‑GGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCG‑TTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTA  >  NC_000913/71591‑72055

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: ‑