Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	92,573	Δ1 bp	coding (1161/1767 nt)	ftsI →	peptidoglycan DD‑transpeptidase FtsI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	92,570	0	A	.	76.5%	12.4 / 1.0	17	coding (1158/1767 nt)	ftsI	peptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base . (13/0);  total (16/1)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGGCGATGCGTAACCGTACCATCACCGACGTGTTTGAACCGGGCTCAACGGTTAAACCGATGGTGGTAATGACCGCGTTGCAACGTGGCGTGGTGCGGGAAAACTCGGTACTCAATACCATTCCTTATCGAATTAACGGCCACG‑AAATCAAAGACGTGGCACGCTACAGCGAATTAACCCTGACC‑GGGGTATTACAG‑AAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCG‑TTAGCGATGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑TTACTCACG‑TTTTGGACT‑GGGAAAAGCGACC‑AA‑TTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAAC‑AACGGTGGT‑CTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGA  >  NC_000913/92285‑92868                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                               agagGCGATGCGTAACCGTACCAT‑ACCGACGTGTTTGAACCGGGCTCAACGGTTAAACCGATGGTGGTAATGACCGCGTTGCAACGTGGCGTGGTGCGGGAAAACTCGGTACTCAATACCATTCCTTATCGAATTAACGGCCACG‑AAATCAAAGACGTGGCACGCTACAGCGAATTAACCCTGACC‑GGGGTATTACAG‑AAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCG‑TTAGCGATGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑TTACTCACG‑TTTTGGACT‑GGG‑AAAGCGACC‑AA‑TTTGGGGTTGGTCGgag                                                                                                                                                                                                                                                                                  >  1:718724/1‑312 (MQ=255)                                                       tAAACCGATGGTGGTAATGACCGCGTTGCAACGTGGCGTGGTGCGGGAAAACTCGGTACTCAATACCATTCCTTATCGAATTAACGGCCACG‑AAATCAAAGACGTGGCACGCTACAGCGAATTAACCCTGACC‑GGGGTATTACAG‑AAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCG‑TTAGCGATGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑TTACTCACG‑TTTTGGACTGGGGAAAAGCGACC‑AA‑TTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATc                                                                                                                                                                                                                                                                >  1:718727/1‑279 (MQ=255)                                                            cGATGGTGGTAATGACCGCGTTGCAACGTGGCGTGGTGCGGGAAAACTCGGTACTCAATACCA‑TCCTTATCGAATTAACGGCCACG‑AAATCAAAGACGTGGCACGCTACAGCGAATTAA‑CCTGACC‑GGGGTATTACAG‑AAGTCGAGTAACGTC‑GTGTTTCCAAGCTGGCG‑TTAGCGATGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑TTACTCACGTTTTTGGACT‑GGGAAAAGCGACCAAA‑TTTGGGGTTg                                                                                                                                                                                                                                                                                         >  1:9780/1‑247 (MQ=255)                                                                                                  cGGGAAAACTCGGTACTCAATACCATTCCTTATCGAATTAACGGCCACG‑AAATCAAAGACGTGGCACGCTACAGCGAATTAACCCTGACC‑GGGGTATTACAG‑AAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTTAGCGATGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑TTACTCACG‑TTTTGGACT‑GGG‑AAAGCGACC‑AA‑‑TTGGGGTTGGTCGGAGAACGCAGTGGCTTATATc                                                                                                                                                                                                                                                                >  1:718734/1‑234 (MQ=255)                                                                                                                   cAATACCATTCCTTATCGAATTAACGGCCACG‑AAATCAAAGACGTGGCACGCTACAGCGAATTAACCCTGACC‑GGGGTATTACAG‑AAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCG‑TTAGCGATGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑TTACTCACG‑TTTTGGACT‑GGGaaaa                                                                                                                                                                                                                                                                                                             <  1:718730/175‑1 (MQ=255)                                                                                                                   cAATACCATTCCTTATCGAATTAACGGCCACG‑AAATCAAAGACGTGGCACGCTACAGCGAATTAACCCTGACC‑GGGGTATTACAG‑AAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCG‑TTAGCGATGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑TTACTCACG‑TTTTGGACT‑GGG‑AAAGCGACCAAA‑TTTGGGGTTGGTCGGAGAACGCAGTgg                                                                                                                                                                                                                                                                        >  1:718736/1‑210 (MQ=255)                                                                                                                   cAATACCATTCCTTATCGAATTAACGGCCACG‑AAATCAAAGACGTGGCACGCTACAGCGAATTAACCCTGACC‑GGGGTATTACAG‑AAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCG‑TTAGCGATGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑TTACTCACG‑TTTTGGACT‑GGG‑AAAGCGACC‑AA‑‑TTGGGGTTGGTCGGAGAACGCAGTg                                                                                                                                                                                                                                                                         >  1:718737/1‑207 (MQ=255)                                                                                                                   cAATACCA‑TCCTTATCGAATTAACGGCCACGAAAATCAAAGACGTGGCACGCTACAGCGAATTAACCCTGACC‑GGGGTATTACAG‑AAG‑CGAGTAACGTCGGTGTTTCCAAGCTGGCG‑TTAGCGATGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑TTACTCACG‑TTTTGGACT‑GGG‑AAAGCGACC‑AA‑TTTGGGGTTGGTCgg                                                                                                                                                                                                                                                                                    >  1:9783/1‑196 (MQ=255)                                                                                                                                                                             cGAATTAACCCTGACC‑‑GGGTATTACAG‑AAGTCGAGTAACGTCGGTG‑TTCCAAGCTGGCG‑TTAGCGATGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑TTACTCACG‑TTTTGGACT‑GGGaaa                                                                                                                                                                                                                                                                                                              <  1:9787/115‑1 (MQ=255)                                                                                                                                                                                     cccTGACCGGGGGTATTACAGAAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTTAGCGATGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑TTACTCACG‑TTTTGGACT‑GGG‑AAAGCGACC‑AA‑TTTgg                                                                                                                                                                                                                                                                                              >  1:9789/1‑125 (MQ=255)                                                                                                                                                                                                 gTATTACAG‑AAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCG‑TTAGCGATGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑TTACTCACGTTTTTGGACT‑GGG‑AAAGCGACC‑AA‑TTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAAC‑AACGGTGGT‑CTGACATAGAGAGGGCCACCTTCTCTTTCGGCTAc                                                                                                                                                                                                         >  1:718741/1‑195 (MQ=255)                                                                                                                                                                                                      aCAG‑AAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCG‑TTAGCGATGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑‑TACTCACG‑TTTTGGACT‑GGG‑AAAGCGACC‑AATTTTGGGGTTg                                                                                                                                                                                                                                                                                         >  1:9791/1‑111 (MQ=255)                                                                                                                                                                                                                                               aGCGATGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑TTACTCACG‑TTTTGGACT‑GGG‑AAAGCGACC‑AA‑TTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAAC‑AACGGTGGT‑CTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACgtg                                                                                      >  1:718743/1‑265 (MQ=255)                                                                                                                                                                                                                                                   aTGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑TTACTCACGTTTTTGGACT‑GGG‑AAAGCGACC‑AA‑‑TTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAAC‑AACGGTGGT‑CTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCa                                                                                                                                            >  1:718744/1‑207 (MQ=255)                                                                                                                                                                                                                                                       cGTCCTTCAGCGTTAAGTAGATAACTTTACTCACGTTTTTGGACT‑GGGAAAAGCGACC‑AA‑ATTTGGGGTGGTCGgaga                                                                                                                                                                                                                                                                                 <  1:9793/78‑1 (MQ=255)                                                                                                                                                                                                                                                                  gTT‑AGTAGAT‑AC‑TTACTCACGTTTTTGGACT‑GGG‑AAAGCGACC‑AA‑TTTGGGGTTGGTCGGAGaa                                                                                                                                                                                                                                                                                >  1:9795/1‑64 (MQ=255)                                                                                                                                                                                                                                                                              ac‑ttactCACG‑TTTTGGACT‑GGG‑AAAGCGACC‑AA‑‑TTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAAACGGTGGT‑CTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGa  >  1:718745/1‑321 (MQ=255)                                                                                                                                                                                                                                                                                  tactCACGTTTTTGGACT‑GGGAAAAGCGACC‑AA‑TTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTC‑AAAAC‑AACGATGGTGCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGcc                                                                                                                                   >  1:9796/1‑191 (MQ=255)                                                                                                                                                                                                                                                                                  tactCACG‑TTTTGGACT‑GGG‑AAAGCGACC‑AA‑TTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAAC‑AACGGTGGT‑CTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTa                              >  1:718746/1‑290 (MQ=255)                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                               AGAGGCGATGCGTAACCGTACCATCACCGACGTGTTTGAACCGGGCTCAACGGTTAAACCGATGGTGGTAATGACCGCGTTGCAACGTGGCGTGGTGCGGGAAAACTCGGTACTCAATACCATTCCTTATCGAATTAACGGCCACG‑AAATCAAAGACGTGGCACGCTACAGCGAATTAACCCTGACC‑GGGGTATTACAG‑AAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCG‑TTAGCGATGCCGTCC‑TCAGCGTT‑AGTAGAT‑AC‑TTACTCACG‑TTTTGGACT‑GGGAAAAGCGACC‑AA‑TTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAAC‑AACGGTGGT‑CTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGA  >  NC_000913/92285‑92868

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