Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	730,644	A→G	R354R (CGA→CGG)	rhsC →	rhs element protein RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	730,644	0	A	G	100.0%	87.5 / NA	30	R354R (CGA→CGG)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/18);  total (12/18)

GCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTG  >  NC_000913/730518‑730779                                                                                                                               |                                                                                                                                        gcgcTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGtt                                                                                                                            <  2:401890/140‑1 (MQ=255) gcgcTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGtt                                                                                                                            <  2:282968/140‑1 (MQ=255)    cTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTAcc                                                                                                                         >  2:776315/1‑140 (MQ=255)         gCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGG‑GCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACg                                                                                                                   >  1:123551/1‑140 (MQ=255)                 ccccGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGACAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTAc                                                                                                                          <  2:56242/126‑1 (MQ=255)                 ccccGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGACAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTAc                                                                                                                          >  1:56242/1‑126 (MQ=255)                  cccGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGAc                                                                                                           <  2:1127338/140‑1 (MQ=255)                     gcggcgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGAcggg                                                                                                        <  1:673107/140‑1 (MQ=255)                            aCTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGa                                                                                                 >  1:174471/1‑140 (MQ=255)                            aCTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGa                                                                                                 <  2:835201/140‑1 (MQ=255)                                         tATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTaaa                                                                                    <  1:776315/140‑1 (MQ=255)                                              ccGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCgg                                                                               <  1:827839/140‑1 (MQ=255)                                              ccGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCgg                                                                               <  2:57847/140‑1 (MQ=255)                                                 cagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAg                                                                            >  1:1052120/1‑140 (MQ=37)                                                   gcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGc                                                                          <  2:1052120/140‑1 (MQ=25)                                                   gcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGc                                                                          >  2:207350/1‑140 (MQ=25)                                                           caGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTa                                                                  <  1:832203/140‑1 (MQ=17)                                                            cagggcgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTac                                                                 <  2:123551/136‑1 (MQ=16)                                                              gTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTacac                                                               >  1:848239/1‑140 (MQ=17)                                                              gTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTacac                                                               >  2:100015/1‑140 (MQ=17)                                                                  gcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTAt                                                           <  1:720174/140‑1 (MQ=18)                                                                  gcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTAt                                                           <  1:207350/140‑1 (MQ=17)                                                                  gcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTAt                                                           <  1:1112055/140‑1 (MQ=17)                                                                                     tAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGca                                        >  1:9115/1‑140 (MQ=21)                                                                                      aaaTACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcat                                       <  1:100015/140‑1 (MQ=37)                                                                                            cGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccat                                 >  1:921708/1‑140 (MQ=21)                                                                                             ggggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatc                                >  1:533700/1‑140 (MQ=21)                                                                                              gggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatca                               >  2:1079067/1‑140 (MQ=21)                                                                                                  cGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGa                           <  2:9115/140‑1 (MQ=21)                                                                                                  cGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGa                           <  1:715637/140‑1 (MQ=25)                                                                                                                           gccggccgGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTg  >  1:470957/1‑140 (MQ=255)                                                                                                                           gccggccgGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTg  >  1:1032910/1‑140 (MQ=255)                                                                                                                               |                                                                                                                                        GCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTG  >  NC_000913/730518‑730779

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