Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	732,698	A→G	N1039S (AAC→AGC)	rhsC →	rhs element protein RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	732,698	0	A	G	100.0%	56.9 / NA	25	N1039S (AAC→AGC)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/15);  total (10/15)

GACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCG  >  NC_000913/732565‑732832                                                                                                                                      |                                                                                                                                       gACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAg                                                                                                                                  <  1:122400/140‑1 (MQ=21)     cGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAAcc                                                                                                                              <  2:940493/140‑1 (MQ=21)     cGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAAcc                                                                                                                              <  1:822143/140‑1 (MQ=21)      gCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAAccc                                                                                                                             >  1:434361/1‑140 (MQ=21)            ccGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATcc                                                                                                                       <  2:447779/140‑1 (MQ=21)             cGCAGGGTGGCAAAACTGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCa                                                                                                                      <  2:1113396/140‑1 (MQ=21)             cGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCa                                                                                                                      <  2:262814/140‑1 (MQ=21)               cAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCAGa                                                                                                                    >  2:505993/1‑140 (MQ=21)               cAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCAGa                                                                                                                    >  2:36726/1‑140 (MQ=21)               cAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCAGa                                                                                                                    <  1:1057599/140‑1 (MQ=21)                 gggTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCAGACg                                                                                                                  >  2:1040118/1‑140 (MQ=25)                 gggTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCAGACg                                                                                                                  <  2:511965/140‑1 (MQ=21)                    tGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCAGACGAtt                                                                                                               >  1:1380813/1‑140 (MQ=25)                       cAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCAGACGATTTAt                                                                                                            <  2:651168/140‑1 (MQ=21)                          aaCGGGTATGGCGGCGTGAACGGGACCTGATGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCAGACGATTTATCAg                                                                                                         >  1:465849/1‑140 (MQ=255)                          aaCGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGGACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGCCAGAACCCGCATCCAGACGATTTATCAg                                                                                                         >  1:1195058/1‑140 (MQ=255)                                tATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCAGACGATTTATCAGCCgggg                                                                                                   >  2:1204386/1‑140 (MQ=21)                                tATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCAGACGATTTATCAGCCgggg                                                                                                   >  2:1266973/1‑140 (MQ=21)                                  tGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAg                                                                                                 <  1:787408/140‑1 (MQ=17)                                  tGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAg                                                                                                 <  2:1170299/140‑1 (MQ=17)                                       gcgTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCa                                                                                            <  1:299876/140‑1 (MQ=11)                                       gcgTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCa                                                                                            <  1:1266973/140‑1 (MQ=11)                                       gcgTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAGCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCa                                                                                            <  1:1204386/140‑1 (MQ=11)                                                                                                                             cGATACAGAGCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGt      >  2:425925/1‑140 (MQ=11)                                                                                                                                 aCAGAGCGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCg  <  1:800312/140‑1 (MQ=18)                                                                                                                                      |                                                                                                                                       GACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGCCTGACCACGATACAGAACGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCG  >  NC_000913/732565‑732832

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