Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	730,654	T→C	C358R (TGT→CGT)	rhsC →	rhs element protein RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	730,654	0	T	C	100.0%	71.2 / NA	22	C358R (TGT→CGT)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/10);  total (12/10)

GCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTT  >  NC_000913/730518‑730709                                                                                                                                         |                                                        gcgcTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGtt                                                      <  1:360387/140‑1 (MQ=255) gcgcTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGAGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGtt                                                      <  2:909750/140‑1 (MQ=255)    cTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTAcc                                                   >  1:706001/1‑140 (MQ=255)       tGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATGCGTTACCGtt                                                >  1:1303228/1‑140 (MQ=255)        ggCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATccgttaccgtta                                               <  1:716809/140‑1 (MQ=255)             gACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACa                                          >  2:496110/1‑140 (MQ=255)               cGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGc                                        >  1:488958/1‑140 (MQ=255)                 ccccGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGa                                      >  2:50895/1‑140 (MQ=255)                  cccGCGGCGAACTGGCGGCGGTTTATGACCTGAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGTATGGTGGCGCACCGTGACACGGGCCGACCGGAAATCCGTTAGCGTTACGACAGCgtc                                     >  2:221266/1‑138 (MQ=14)                   ccGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGAcg                                    >  2:1325144/1‑140 (MQ=255)                                 cggcggTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGaa                      <  2:706001/140‑1 (MQ=255)                                              ccGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCgg         <  1:1334278/140‑1 (MQ=255)                                              ccGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCgg         <  1:627571/140‑1 (MQ=255)                                              ccGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCgg         <  2:28505/140‑1 (MQ=255)                                                 cagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCa       <  2:554710/139‑1 (MQ=37)                                                 cagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCa       >  1:554710/1‑139 (MQ=37)                                                 cagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAg      >  1:173198/1‑140 (MQ=37)                                                 cagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAg      >  1:650261/1‑140 (MQ=37)                                                 cagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAg      >  1:400778/1‑140 (MQ=37)                                                   gcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCAGACGGGCCGACCGGAAATCCGtt                                                      >  2:1271445/1‑90 (MQ=255)                                                   gcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGtt                                                      <  1:1271445/90‑1 (MQ=255)                                                   gcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGc    <  2:227184/140‑1 (MQ=25)                                                    caATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCt   >  1:1105065/1‑140 (MQ=21)                                                     aaTACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCtt  <  1:178868/140‑1 (MQ=18)                                                                                                                                         |                                                        GCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTT  >  NC_000913/730518‑730709

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: ‑