Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,762,200	A→G	A6A (GCA→GCG)	rhsA →	rhs element protein RhsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,762,200	0	A	G	100.0%	80.4 / NA	29	A6A (GCA→GCG)	rhsA	rhs element protein RhsA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/13);  total (16/13)

TTTCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCC  >  NC_000913/3762076‑3762287                                                                                                                             |                                                                                        tttCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATTAAACAAGGCTGATAAAGGAGCGGAAAACCGGCGGCGCGTCa                                                                                 >  2:200421/1‑133 (MQ=18) tttCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATCAAACAAGGGTGAAAAATGAGCGGAAGACCGGCGGCGCGTCa                                                                                 >  2:166961/1‑133 (MQ=18) tttCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCa                                                                                 >  2:84733/1‑133 (MQ=255) tttCCTAATTATTTAACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAATGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCa                                                                                 >  1:421698/1‑133 (MQ=255)   tCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAgg                                                                               >  1:16544/1‑133 (MQ=255)    ccTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAATGAGCGGAATACAGGCGGCTCGTCAggg                                                                              >  1:239271/1‑133 (MQ=17)                    tGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTAt                                                              <  2:423194/133‑1 (MQ=255)                    tGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTAt                                                              >  1:476500/1‑133 (MQ=255)                         gATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAATGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCgg                                                         <  2:421698/133‑1 (MQ=255)                         gATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCgg                                                         >  2:254656/1‑133 (MQ=255)                             tctcACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGc                                                     <  1:426570/133‑1 (MQ=255)                              ctcACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCa                                                    >  2:126150/1‑133 (MQ=255)                                   ttCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGtt                                               >  2:66054/1‑133 (MQ=255)                                   ttCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGtt                                               >  1:236561/1‑133 (MQ=255)                                    tCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTc                                              >  1:459486/1‑133 (MQ=255)                                     cAGGCTTTCTTATAACTCTGTAGGATTTTGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCa                                             <  2:239271/133‑1 (MQ=11)                                              ttATAAATCTTTAGGGTTTCGCCTATCAGCAGGCAAATAACTCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGc                                    <  2:189912/133‑1 (MQ=11)                                              ttATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGc                                    <  2:16544/133‑1 (MQ=255)                                              ttATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGc                                    <  2:236561/133‑1 (MQ=255)                                               tataAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATATAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGcc                                   <  1:479714/133‑1 (MQ=255)                                               tataAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGcc                                   <  1:107295/133‑1 (MQ=255)                                               tataAATCTGTAGGGTTTAGCCGGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGcc                                   <  1:200421/133‑1 (MQ=255)                                                        gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTacgt                          >  1:506290/1‑132 (MQ=255)                                                        gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAACCATGTTCAGGGTTTACCCgggggaccg                          >  2:395273/1‑128 (MQ=11)                                                              tttCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGCCGGGGGACGTATTGGt                    >  1:31947/1‑133 (MQ=11)                                                                 cGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTGCAGGGGTGAGCGGGGGGACGGGTTGGTGcc                 >  1:163701/1‑133 (MQ=11)                                                                     tGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCa             <  2:506290/133‑1 (MQ=11)                                                                     tGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCa             <  2:53415/133‑1 (MQ=255)                                                                        cagcagACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGAGGATGTGTTGGTGACCCCACCg          >  1:203749/1‑133 (MQ=11)                                                                            agACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGt      <  1:320270/133‑1 (MQ=255)                                                                                aaaTAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGGATTGTTTACGGTTCAGCCGGGGGAAGTATTGGTGCCCCCCACTGCGTGGcc  >  2:316764/1‑133 (MQ=11)                                                                                                                             |                                                                                        TTTCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCC  >  NC_000913/3762076‑3762287

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