Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 291,324 C→T C8Y (TGT→TAT)  insA ← IS1 repressor TnpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913291,3240CT100.0% 55.0 / NA 19C8Y (TGT→TAT) insAIS1 repressor TnpA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/10);  total (9/10)

TGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGC  >  NC_000913/291267‑291459
                                                         |                                                                                                                                       
tGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAAt                                                        >  2:179978/1‑139 (MQ=14)
         gTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAATACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGccc                                               >  2:1427618/1‑139 (MQ=255)
                    gTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGcc                                                >  2:611630/1‑127 (MQ=255)
                    gTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGcc                                                <  1:611630/127‑1 (MQ=255)
                        cGCACCACCCCTTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTAc                                <  2:1046597/139‑1 (MQ=255)
                           accaccCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTa                             <  2:204282/139‑1 (MQ=255)
                                  cGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGaa                      >  2:473473/1‑139 (MQ=255)
                                     cAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCTTCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAatat                   <  1:407356/139‑1 (MQ=255)
                                     cAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATg                                          <  2:932189/116‑1 (MQ=255)
                                     cAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATg                                          >  1:932189/1‑116 (MQ=255)
                                     cAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAatat                   <  1:179978/139‑1 (MQ=255)
                                       gTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCt                 >  2:772240/1‑139 (MQ=255)
                                        tAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGAt                                           <  1:1207204/112‑1 (MQ=255)
                                        tAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGAt                                           >  2:1207204/1‑112 (MQ=255)
                                             gAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAATACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGa           <  1:1427618/139‑1 (MQ=255)
                                              aaCAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGAt          >  2:235378/1‑139 (MQ=255)
                                              aaCAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGAt          >  2:1053136/1‑139 (MQ=255)
                                               aCAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATa         <  1:1053136/139‑1 (MQ=255)
                                                      ggATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGc  <  1:235378/139‑1 (MQ=255)
                                                         |                                                                                                                                       
TGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGC  >  NC_000913/291267‑291459

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: