Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,033,271	C→T	V51M (GTG→ATG)	drpB ←	cell division protein DrpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,033,271	0	C	T	100.0%	75.0 / NA	27	V51M (GTG→ATG)	drpB	cell division protein DrpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/15);  total (12/15)

TTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATGTATAGAAAGCCCACAAAGCCAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTC  >  NC_000913/2033046‑2033507                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             tttCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCtcat                                                                                                                                                                                                                      >  1:539046/1‑250 (MQ=255)  ttCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCtcatc                                                                                                                                                                                                                     <  1:177471/250‑1 (MQ=255)      aaCTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCACATCCCGACATGGGGGCATCGGGAATCCGGCtcatcaaatc                                                                                                                                                                                                                >  1:232523/1‑247 (MQ=255)          ggtcctgccAGGGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGGCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCTCATCACCCAGATg                                                                                                                                                                                                            <  2:5552/248‑1 (MQ=255)                 ccAGTGTTTTAAGGCCTTACGCGCTCCGTTGCCAAAGAGTCGCCTCCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACTAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCAGAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAAcgcg                                                                                                                                                                                                     <  2:275547/251‑1 (MQ=255)                            aGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAACTCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAACTCCCGCCATGGGGGCATCGGGAATCCGGCTCATTACCCAGATGTAACGCGCCATTGcaccg                                                                                                                                                                                          >  2:385674/1‑247 (MQ=255)                                           cGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCTCATCACCCGGATGTAACGCGCCATTGCCCAGACGAAATAGCCacac                                                                                                                                                                           >  1:478314/1‑251 (MQ=255)                                               gCCAAAGAGTCGCCTCCTTAGCAGTTTTTATTCACAGCTTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAACCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCAt                                                                                                                                                                        <  2:405862/250‑1 (MQ=255)                                                ccAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATGt                                                                                                                                                                      >  1:155069/1‑251 (MQ=255)                                                              ccTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATGTATAGAAAGCCcaca                                                                                                                                                        <  1:3704/251‑1 (MQ=255)                                                                       tttttATTCATCGCGTCTGCTACGCGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACCATGCCCCGACCAGAGCCATAAATAACAATCCGACAAGTGCCCCTAACCCCTTTCCCGCGGTAGCGCCTAAATCGCTTTCAAATACCGACATGGTGGCATCGGGACTCCGGCTCATCCCCCAGATGTAACGCGCCATTGCCCAGCCGGAAAAGCCACACCCTGAATAGCAAGCCGACAAAGACAGtt                                                                                                                                               >  2:97078/1‑251 (MQ=255)                                                                          ttATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCggg                                                                                                                                                                                                                                  <  1:423072/165‑1 (MQ=255)                                                                          ttATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCggg                                                                                                                                                                                                                                  >  2:423072/1‑165 (MQ=255)                                                                             ttCATAGCGTCTGCTTCGTGCGGGACGAAGCCGGGTGTACCAGGCCATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGAGGTAGAGCCTAAATCGCTTTCAAATACCGCCATGGGGGCATCGGGAATCCGGCTCATCACCCACATGTAACGCGCCCTTGCCCAGCCGAAATAGCCACACCATGTATAGCACACCCACACAACCAGTTTgccgc                                                                                                                                         >  2:44928/1‑251 (MQ=255)                                                                                                              ggtttACCAGGCAATGCAGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGTGGTAGAGCCGAAATCGCTTTCAAAGCCCGACATGTGGGCATCGGGAATCCGGCTCATCCCCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATGTATAGAAAGCCCACAAAGCCAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTg                                                                                                        <  2:486057/247‑1 (MQ=255)                                                                                                              ggTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATGTATAGAAAGCCCACAAAGCCAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTg                                                                                                        >  1:189688/1‑251 (MQ=255)                                                                                                                   atcAGGCAATGCCGCCTAACTATGCCCCGACCCGATCCATAAATAAAAAGCCGACAAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTACAAATCCCGACTTGGGGGCATCGGGAATCCGGCTCATCACCCAGATGGAACGCGCCATTCCCCAGACGAAATAGCCACACCATGTATAGAAAGGCCACAAAGCCTGTTTGCCGCCAGGGCTTCGAGAGAGTCTTTCTTCCATCTTTGTCGaa                                                                                                   <  2:539046/249‑1 (MQ=255)                                                                                                                        gCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACTCCATGTATAGAAAGccc                                                                                                                                                           <  1:117760/190‑1 (MQ=255)                                                                                                                        gCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATGTATAGACAGccc                                                                                                                                                           >  2:117760/1‑190 (MQ=255)                                                                                                                             gctgcCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGTCATGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATGTATAGAAAGCCCACCAAGCCAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTcc                                                                                          >  2:450031/1‑250 (MQ=255)                                                                                                                                              cGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATGTATAGAAAGCCCACAAAGCCAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGa                                                                        <  1:230764/251‑1 (MQ=255)                                                                                                                                                           aaataaaAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATGTATAGAAAGCCCACAAAGCCAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACtt                                                            <  1:389223/250‑1 (MQ=255)                                                                                                                                                               aaaaaTCCGACCAGGCCCCCTAACCCCTTTCCCGCGGTAGAGCCCAAATCGATTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATGTATATAAAGCCCACAAAGCCAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTa                                                        <  2:204810/250‑1 (MQ=255)                                                                                                                                                                                           ttCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATGTATAGAAAGCCCACAAAGCCAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTcc                           >  1:548009/1‑251 (MQ=255)                                                                                                                                                                                              ccGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATGTATAGAAAGCCCACAAAGCCAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCcaca                        <  1:180845/251‑1 (MQ=255)                                                                                                                                                                                                                cGCTTTCAAATCCCGACATGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATGTATAGAAAGCCCACAAAGCCAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAg      <  1:476868/251‑1 (MQ=255)                                                                                                                                                                                                                     tCAAATCCCGACATGGGGGCATCGGGAATCCGGCGCATCACCCAGAGGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATGTATAGAAAGCCCACAAAGCCAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTc  <  1:384047/250‑1 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             TTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATGTATAGAAAGCCCACAAAGCCAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTC  >  NC_000913/2033046‑2033507

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