Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	731,919	A→G	G779G (GGA→GGG)	rhsC →	rhs element protein RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	731,919	0	A	G	100.0%	55.5 / NA	31	G779G (GGA→GGG)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/10);  total (21/10)

GAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGA  >  NC_000913/731782‑731960                                                                                                                                          |                                          gAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGTAGCATGAAACCGGAAATGCGTACAACGCGCAGGGGCt                                         >  1:478355/1‑140 (MQ=18)   aGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTgg                                       >  1:1038874/1‑140 (MQ=18)     ggCACCGGGTGACGGCGCACTAAGGCTATCACAGTAATGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCg                                     <  1:435077/140‑1 (MQ=17)       cACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGaa                                   <  2:273064/140‑1 (MQ=17)       cACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGaa                                   <  1:331357/140‑1 (MQ=18)          cGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCg                                <  2:1243675/140‑1 (MQ=18)            ggTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCt                              <  2:174736/140‑1 (MQ=18)            ggTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCt                              <  2:1384212/140‑1 (MQ=17)                       aCTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGTACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGAc                   <  2:1108442/140‑1 (MQ=17)                        cTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACa                  >  1:1078980/1‑140 (MQ=17)                           tGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGTGAACCGCTGTATACCGGACAGcc               >  1:1437614/1‑140 (MQ=14)                             gTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTg             >  1:1013213/1‑140 (MQ=14)                              ttATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGc            >  2:176902/1‑140 (MQ=14)                                  gACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGAGGCTGGCGAACGGCAGTATACAGGACAGCCTGCccgc        >  1:448301/1‑140 (MQ=14)                                   aCAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCccgcc       >  1:1408995/1‑140 (MQ=14)                                     aGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGt     >  2:980841/1‑140 (MQ=11)                                     aGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGt     >  1:301775/1‑140 (MQ=11)                                       tAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATCCCGGACAGCCTGCCCGCCGTgg   >  1:699402/1‑140 (MQ=11)                                       tAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCccgc        >  1:1487523/1‑135 (MQ=11)                                       tAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCccgc        <  2:1487523/135‑1 (MQ=11)                                       tAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTgg   >  1:36136/1‑140 (MQ=11)                                       tAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTgg   >  1:826182/1‑140 (MQ=11)                                       tAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTgg   <  1:980841/140‑1 (MQ=11)                                       tAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTgg   >  1:1531134/1‑140 (MQ=11)                                       tAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTgg   >  1:1374378/1‑140 (MQ=11)                                       tAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTgg   >  1:1157150/1‑140 (MQ=11)                                       tAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTgg   >  1:1105701/1‑140 (MQ=11)                                       tAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCCCTGCATGCCGGACTGCCTGCCCGCCGTgg   >  1:710746/1‑140 (MQ=11)                                        aaaGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGcc           >  2:1159202/1‑131 (MQ=11)                                        aaaGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGcc           <  1:1159202/131‑1 (MQ=11)                                        aaaGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGa  >  1:354930/1‑140 (MQ=11)                                                                                                                                          |                                          GAAGGGCACCGGGTGACGGTGCACTATGGTTATGACAGTAAAGGCCGCCTCGCCAGTGAACACCTGACGGTGCATCATCCGCAGACGAATGAACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGA  >  NC_000913/731782‑731960

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