Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	20,487	A→G	C8R (TGT→CGT)	insA ←	IS1 repressor TnpA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	20,487	0	A	G	100.0%	84.4 / NA	25	C8R (TGT→CGT)	insA	IS1 repressor TnpA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/17);  total (8/17)

CCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGA  >  NC_000913/20432‑20618                                                        |                                                                                                                                    ccGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCa                                                  <  1:295717/139‑1 (MQ=21) ccGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCa                                                  <  2:59300/139‑1 (MQ=21) ccGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCa                                                  <  2:463821/139‑1 (MQ=21)       gTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATgtgt                                            <  1:485373/139‑1 (MQ=255)          cTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATc                                         >  2:42698/1‑139 (MQ=255)              tGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACaa                                     <  1:128046/139‑1 (MQ=255)              tGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACaa                                     <  2:855726/139‑1 (MQ=255)                    tACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTAt                                          >  1:216017/1‑128 (MQ=255)                    tACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTAt                                          <  2:216017/128‑1 (MQ=255)                    tACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCa                               <  1:448985/139‑1 (MQ=255)                        caccacCCCGTCAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATAtt                           <  2:175427/139‑1 (MQ=255)                             ccccGTCAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTtg                      >  2:707516/1‑139 (MQ=255)                                cGTCAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTgg                   >  2:117145/1‑139 (MQ=255)                                 gTCAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACGATATCCATATTCTTTGTggg                  <  1:715517/139‑1 (MQ=255)                                  tCAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTgggg                 <  1:529764/139‑1 (MQ=255)                                   cAGTAGCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTggggg                <  2:134148/139‑1 (MQ=255)                                        gCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCt           <  1:368076/139‑1 (MQ=255)                                        gCTGAACAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGATGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCt           >  2:640217/1‑139 (MQ=255)                                            aaCAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATcc                                >  1:631649/1‑114 (MQ=255)                                            aaCAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATcc                                <  2:631649/114‑1 (MQ=255)                                             aCAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGgaga      <  1:880363/139‑1 (MQ=255)                                             aCAGGAGGGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGgaga      <  2:1032472/139‑1 (MQ=255)                                                ggaggGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTg   <  1:899088/139‑1 (MQ=255)                                                 gaggGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGa  >  1:701882/1‑139 (MQ=255)                                                 gaggGACGGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGa  >  1:633162/1‑139 (MQ=255)                                                        |                                                                                                                                    CCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATGTGTATCACAATATCCATATTCTTTGTGGGGGAGTCTGGAGATTGA  >  NC_000913/20432‑20618

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