Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,762,260	G→A	V26V (GTG→GTA)	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,260	0	G	A	100.0%	22.7 / NA	11	V26V (GTG→GTA)	rhsA	rhs element protein RhsA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (5/6);  total (5/6)

TGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCC  >  NC_000913/3762130‑3762325                                                                                                                                   |                                                                  tGTAGGGTTTAGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTTTTTAGGGGTCAGCCGGGGGACGTATTGGt                                                          >  2:394489/1‑140 (MQ=14)   tAGGGTTTAGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTTAGGGGTCAGCCGGGGGACGTATTGGTGc                                                        >  2:284548/1‑140 (MQ=14)        tttCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCa                                                   <  2:89694/140‑1 (MQ=255)        tttCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGACGTATTGGTGCCCCCa                                                   >  2:267850/1‑140 (MQ=255)        tttCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGCCGGGGGACGTATTGGTGcccccc                                                   >  2:317087/1‑139 (MQ=255)             ccTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGCCGGGGGACGTATTGGTGCCCCCACCCGc                                              >  1:3324/1‑140 (MQ=14)                gTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTg                                           >  2:46205/1‑140 (MQ=255)                 tCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTgg                                          >  1:598319/1‑140 (MQ=255)                               aCCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGc                            <  2:598319/140‑1 (MQ=255)                                 ccGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGccc                          <  1:558785/140‑1 (MQ=255)                                    aTAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGCCGGGGGACGTATTGGTGCCCCCCCCGGCGTGGCCTGTTCGGTGTGGCCcgg                       >  1:572479/1‑140 (MQ=14)                                              gATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACg             <  1:808247/140‑1 (MQ=255)                                                tGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTc           <  1:711531/140‑1 (MQ=21)                                                  agaAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCg         <  2:300868/140‑1 (MQ=18)                                                         aGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCTGGGCTCAGCCGGGGCACGCATTGGGGCCCCCCCCGGCGTGGCCCGTTCGGTGTGCCCCCGCGGAGGGGCCGCCGGcccccc  >  2:652222/1‑136 (MQ=11)                                                                                                                                   |                                                                  TGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCC  >  NC_000913/3762130‑3762325

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