Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,619,209	G→A	A6A (GCG→GCA)	rhsB →	rhs element protein RhsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,619,209	0	G	A	100.0%	154.6 / NA	47	A6A (GCG→GCA)	rhsB	rhs element protein RhsB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/22);  total (25/22)

TCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGC  >  NC_000913/3619073‑3619314                                                                                                                                         |                                                                                                          tctcTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGgcagcg                                                                                                        >  2:167011/1‑140 (MQ=255)         gTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAggg                                                                                                <  2:647998/140‑1 (MQ=255)          tAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAgggc                                                                                               >  1:110553/1‑139 (MQ=255)                 cTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGa                                                                                        <  2:590246/140‑1 (MQ=255)                 cTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGa                                                                                        >  1:759099/1‑140 (MQ=255)                   aaaTATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACg                                                                                      <  2:1150444/140‑1 (MQ=255)                        tatTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTa                                                                                 >  2:409773/1‑140 (MQ=255)                         atTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAggg                                                                                                <  2:667479/124‑1 (MQ=255)                         atTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAggg                                                                                                >  1:667479/1‑124 (MQ=255)                         atTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTAt                                                                                <  2:847322/140‑1 (MQ=255)                          ttCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGa                                                                                        <  1:110918/131‑1 (MQ=255)                          ttCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGa                                                                                        >  2:110918/1‑131 (MQ=255)                          ttCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATg                                                                               >  1:394592/1‑140 (MQ=255)                            ccTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGc                                                                             >  1:492447/1‑140 (MQ=255)                              tAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCgg                                                                           >  2:215046/1‑140 (MQ=255)                              tAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCgg                                                                           >  1:922272/1‑140 (MQ=255)                               aaaaaTCTTCAATTCATTGTGACCACAAGTTTATCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCgg                                                                           <  2:931125/139‑1 (MQ=255)                               aaaaaTCTTCAATTCATTGTGACCACAAGTTTATCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCgg                                                                           >  1:931125/1‑139 (MQ=255)                                aaaaTCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTAt                                                                                <  1:280823/133‑1 (MQ=255)                                aaaaTCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTAt                                                                                >  2:280823/1‑133 (MQ=255)                                aaaaTCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTa                                                                         <  2:759099/140‑1 (MQ=255)                                   aTCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCa                                                                      <  2:922272/140‑1 (MQ=255)                                   aTCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCa                                                                      <  2:601045/140‑1 (MQ=255)                                        cAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGtt                                                                 >  1:1074841/1‑140 (MQ=255)                                                tgtgACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTc                                                         >  2:334714/1‑140 (MQ=255)                                                              tttCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCgca                                           >  1:740304/1‑138 (MQ=255)                                                                     gCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGc                                    <  2:110553/140‑1 (MQ=18)                                                                          ttCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCa                               <  1:334714/140‑1 (MQ=18)                                                                          ttCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCa                               <  2:1074841/140‑1 (MQ=255)                                                                          ttCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCa                               <  2:492447/140‑1 (MQ=18)                                                                               aTGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGc                          <  1:409773/140‑1 (MQ=255)                                                                               aTGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGc                          <  1:167011/140‑1 (MQ=255)                                                                                   aGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTgg                      >  2:361971/1‑140 (MQ=255)                                                                                     aTACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGcc                    >  2:489399/1‑140 (MQ=255)                                                                                        cTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGt                 >  1:1114270/1‑140 (MQ=255)                                                                                         tGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGtt                >  1:511834/1‑140 (MQ=255)                                                                                         tGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGCGCATTGGTGCCCCCACCGGCGTGGCCTGtt                >  1:229534/1‑140 (MQ=14)                                                                                         tGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGCCGGGGGGCGCATTGGGGCCCCCACCGGCGTGGCCTGtt                >  2:1095136/1‑140 (MQ=14)                                                                                           tCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCTACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCg              >  1:895439/1‑140 (MQ=14)                                                                                           tCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCg              <  2:1114270/140‑1 (MQ=255)                                                                                             aTTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGgt            >  2:147336/1‑140 (MQ=255)                                                                                               tAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGgtgt          >  1:325447/1‑140 (MQ=14)                                                                                                aaaaTAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGtgtg         >  2:912913/1‑140 (MQ=255)                                                                                                   ataataGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGccc      <  1:215046/140‑1 (MQ=255)                                                                                                   ataataGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGccc      <  2:334403/140‑1 (MQ=255)                                                                                                      ataGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCTACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcgg   <  2:895439/140‑1 (MQ=11)                                                                                                      ataGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcgg   <  1:912913/140‑1 (MQ=255)                                                                                                       taGAAAAGGATTTTACGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCcggc  >  1:985164/1‑140 (MQ=255)                                                                                                                                         |                                                                                                          TCTCTTTTGTAACGTTCTAAATATATTCCTAAAAATCTTCAATTCATTGTGACCACAAGTTTTTCTTCGCTTTTTCGTATGAAGATACTGTCATTAAAATAATAGAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGTGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGC  >  NC_000913/3619073‑3619314

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