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% 74.2 / NA 25C8Y (TGT→TAT) insAIS1 repressor TnpA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/13);  total (12/13)

ATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGC  >  NC_000913/291266‑291459
                                                          |                                                                                                                                       
aTGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGaa                                                         >  2:1015858/1‑139 (MQ=11)
     cGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTc                                                    >  2:1743020/1‑139 (MQ=21)
     cGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTc                                                    <  2:1038424/139‑1 (MQ=21)
        cAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGc                                                 <  2:675934/139‑1 (MQ=37)
               tttGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATg                                          <  1:367503/139‑1 (MQ=255)
                      ttACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACg                                   <  2:1367451/139‑1 (MQ=255)
                       tACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGt                                  >  1:1445655/1‑139 (MQ=255)
                              caccCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATg                           >  2:1583329/1‑139 (MQ=255)
                                ccccGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATg                                          >  1:1409589/1‑122 (MQ=255)
                                ccccGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATg                                          <  2:1409589/122‑1 (MQ=255)
                                ccccGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCCCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTc                         >  2:792009/1‑139 (MQ=255)
                                      cAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGtt                                        <  1:160906/118‑1 (MQ=255)
                                      cAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGtt                                        >  2:160906/1‑118 (MQ=255)
                                        gTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGccc                                               >  1:434149/1‑109 (MQ=255)
                                        gTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGccc                                               <  2:434149/109‑1 (MQ=255)
                                         tAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTg                <  2:1445655/139‑1 (MQ=255)
                                         tAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTg                <  1:1552759/139‑1 (MQ=255)
                                         tAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTg                <  1:1015858/139‑1 (MQ=255)
                                           gCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGc              >  1:214649/1‑139 (MQ=255)
                                              gAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGa           >  2:671080/1‑139 (MQ=255)
                                               aaCAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGAt          <  2:1674530/139‑1 (MQ=255)
                                               aaCAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGAt          <  2:214649/139‑1 (MQ=255)
                                                   ggaggGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAgg      >  2:978966/1‑139 (MQ=255)
                                                       ggATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCg                        >  2:1713245/1‑117 (MQ=255)
                                                       ggATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGc  <  1:1583329/139‑1 (MQ=255)
                                                          |                                                                                                                                       
ATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGC  >  NC_000913/291266‑291459

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: