Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	524,426	A→G	D389G (GAC→GGC)	rhsD →	protein RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	524,426	0	A	G	100.0%	97.5 / NA	28	D389G (GAC→GGC)	rhsD	protein RhsD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/16);  total (12/16)

CGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCAGCGGGTATGACGC  >  NC_000913/524300‑524552                                                                                                                               |                                                                                                                               cGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCAtcaccgtc                                                                                                                   >  2:267496/1‑140 (MQ=255)     ccGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCAtcaccgtcaccg                                                                                                               >  2:679247/1‑140 (MQ=255)      cGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGa                                                                                                              <  1:918390/140‑1 (MQ=255)                cgcACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAAc                                                                                                    >  2:426452/1‑140 (MQ=255)                cgcACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAAc                                                                                                    <  1:679247/140‑1 (MQ=255)                    ccGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGc                                                                                                <  2:865473/140‑1 (MQ=255)                       ttACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTg                                                                                             >  1:117607/1‑140 (MQ=255)                       ttACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTg                                                                                             >  1:1015410/1‑140 (MQ=255)                       ttACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTg                                                                                             <  1:183591/140‑1 (MQ=255)                       ttACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGATGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTg                                                                                             <  2:826464/140‑1 (MQ=255)                         aCGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAg                                                                                           <  1:251271/140‑1 (MQ=255)                         aCGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGACACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAg                                                                                           <  1:810305/140‑1 (MQ=255)                            cgGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGtg                                                                                        <  1:262925/140‑1 (MQ=255)                                aaGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGtgctgc                                                                                    <  1:267496/140‑1 (MQ=255)                                     cGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTg                                                                                             <  1:332298/126‑1 (MQ=255)                                     cGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTg                                                                                             >  2:332298/1‑126 (MQ=255)                                          aTGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAgg                                                                          <  1:759932/140‑1 (MQ=255)                                             cgctaccgctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCgg                                                                       >  1:715143/1‑140 (MQ=255)                                                        cgacgaTACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGc                                                                    <  2:176015/132‑1 (MQ=255)                                                        cgacgaTACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGc                                                                    >  1:176015/1‑132 (MQ=255)                                                                            tggAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGaa                                        <  2:235385/140‑1 (MQ=255)                                                                                         cccGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCa                           >  1:327905/1‑140 (MQ=255)                                                                                                     aaGCTACCGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCa               <  2:765929/140‑1 (MQ=255)                                                                                                            cGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCAGCGGGTa        >  1:169234/1‑140 (MQ=255)                                                                                                            cGCTATCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCAGCGGGTa        >  1:158248/1‑140 (MQ=255)                                                                                                                aTCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCAGCGGGTATGAc    >  1:86523/1‑140 (MQ=255)                                                                                                                aTCTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCAGCGGGTATGAc    <  2:327905/140‑1 (MQ=255)                                                                                                                  cTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCAGCGGGTATGAcgc  >  1:203791/1‑140 (MQ=255)                                                                                                                  cTTTATGAGCAGGGCCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCAGCGGGTATGAcgc  <  2:203791/140‑1 (MQ=255)                                                                                                                               |                                                                                                                               CGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCAGCGGGTATGACGC  >  NC_000913/524300‑524552

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