Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,391,872	G→A	R53C (CGC→TGC)	tldD ←	metalloprotease TldD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,391,872	0	G	A	100.0%	168.3 / NA	47	R53C (CGC→TGC)	tldD	metalloprotease TldD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/23);  total (24/23)

GCGCCGCTTGCGCACTCTGTTCCAGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCA  >  NC_000913/3391740‑3392008                                                                                                                                     |                                                                                                                                         gcgcCGCTTGCGCACTCTGTTCCAGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGtcttct                                                                                                                                   <  2:185070/140‑1 (MQ=255)   gccgcTTGCGCACTCTGTTCCAGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGt                                                                                                                                        >  1:922262/1‑133 (MQ=255)   gccgcTTGCGCACTCTGTTCCAGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTaa                                                                                                                                 >  2:898555/1‑140 (MQ=255)     cgcTTGCGCACTCTGTTCCAGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTaaaa                                                                                                                               >  1:50977/1‑140 (MQ=255)       cTTGCGCACTCTGTTCCAGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAAcc                                                                                                                             <  1:453458/140‑1 (MQ=255)         tGCGCACTCTGTTCCAGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCa                                                                                                                           <  1:1063621/140‑1 (MQ=255)         tGCGCACTCTGTTCCAGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCa                                                                                                                           <  2:1066339/140‑1 (MQ=255)                 ctGTTCCAGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCgg                                                                                                                   >  1:717611/1‑139 (MQ=255)                 ctGTTCCAGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGt                                                                                                                   >  2:629555/1‑140 (MQ=255)                    ttCCAGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGAt                                                                                                                <  1:890124/140‑1 (MQ=255)                        aGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCt                                                                                                            <  1:886005/140‑1 (MQ=255)                        aGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCt                                                                                                            <  1:609721/140‑1 (MQ=255)                        aGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCt                                                                                                            <  2:188741/140‑1 (MQ=255)                        aGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCt                                                                                                            <  1:172749/140‑1 (MQ=255)                          cgcCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCg                                                                                                          >  1:41600/1‑140 (MQ=255)                          cgcCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCg                                                                                                          >  1:351557/1‑140 (MQ=255)                          cgcCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCg                                                                                                          >  1:991841/1‑140 (MQ=255)                            ccAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGt                                                                                                                   >  1:906393/1‑129 (MQ=255)                            ccAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGt                                                                                                                   <  2:906393/129‑1 (MQ=255)                             cagcagGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCGCCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACt                                                                                                       <  2:615959/140‑1 (MQ=255)                             cagcagGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACt                                                                                                       <  1:629555/140‑1 (MQ=255)                              agcagGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTg                                                                                                      >  2:695281/1‑140 (MQ=255)                              agcagGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTg                                                                                                      >  1:1040807/1‑140 (MQ=255)                                        tttGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAgaga                                                                                            <  1:326634/140‑1 (MQ=255)                                         ttGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGAt                                                                                           <  1:346532/140‑1 (MQ=255)                                         ttGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGAt                                                                                           <  1:424421/140‑1 (MQ=255)                                                 cGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATaa                                                                                   >  1:171533/1‑140 (MQ=255)                                                           tCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGAc                                                                         >  2:444784/1‑140 (MQ=255)                                                                     tCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAg                                                               <  2:15619/140‑1 (MQ=255)                                                                      cACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGt                                                              <  1:444784/140‑1 (MQ=255)                                                                      cACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGt                                                              <  2:50977/140‑1 (MQ=255)                                                                                      cacaCCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGa                                              >  1:616205/1‑140 (MQ=255)                                                                                      cacaCCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGa                                              >  1:990510/1‑140 (MQ=255)                                                                                                   tGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATg                                 >  1:769136/1‑140 (MQ=255)                                                                                                   tGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATg                                 >  1:818404/1‑140 (MQ=255)                                                                                                   tGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACATTCGGCCAGTTGACCGAGGATCGCGAAc                                            >  1:88581/1‑129 (MQ=255)                                                                                                   tGATCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACATTCGGCCAGTTGACCGAGGATCGCGAAc                                            <  2:88581/129‑1 (MQ=255)                                                                                                     aTCGATGTTGTAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGtt                               >  1:1075038/1‑140 (MQ=255)                                                                                                          tgttgtAAGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAgg                          >  1:23718/1‑140 (MQ=255)                                                                                                                 aGAACCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTcg                   <  1:551836/140‑1 (MQ=255)                                                                                                                    aCCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTcgccg                >  1:723951/1‑140 (MQ=255)                                                                                                                    aCCATCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTcgccg                >  2:200724/1‑140 (MQ=255)                                                                                                                        tCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAg            >  1:1022596/1‑140 (MQ=255)                                                                                                                        tCTTTAATAATGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAg            <  2:723951/140‑1 (MQ=255)                                                                                                                             aataatGCAGTCTTCTAAAACCCTGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAAtt       <  2:283914/140‑1 (MQ=255)                                                                                                                                aatGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGtt    <  2:41600/140‑1 (MQ=255)                                                                                                                                  tGCAGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCa  >  2:662447/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                         GCGCCGCTTGCGCACTCTGTTCCAGCGCCAGCAGGCTGATTTGGTCAGCGTAAGCAAATCCGGTTTTTTCACCGCTGATTGCACGCACACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCA  >  NC_000913/3391740‑3392008

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑