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

ATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGT  >  NC_000913/291266‑291456
                                                          |                                                                                                                                      
aTGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGaa                                                      >  2:1203240/1‑139 (MQ=11)
aTGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGaa                                                      >  1:788660/1‑139 (MQ=11)
          gTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGccc                                            >  2:1994584/1‑139 (MQ=255)
           tGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCg                                           >  1:2030132/1‑139 (MQ=255)
            gCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCgg                                          <  2:424605/139‑1 (MQ=255)
            gCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCgg                                          <  1:1879287/139‑1 (MQ=255)
               tttGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCGCCGAATTTTCAGCCCGCATg                                       <  2:1398473/139‑1 (MQ=255)
                 tGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGtt                                     <  2:1117190/139‑1 (MQ=255)
                    cGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCa                                  <  2:371184/139‑1 (MQ=255)
                    cGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCa                                  <  2:277201/139‑1 (MQ=255)
                    cGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCa                                  <  2:1451455/139‑1 (MQ=255)
                       tACGCACCACCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCAc                                 >  1:1315588/1‑139 (MQ=255)
                           caccacCCCGTCAGTAGCTGAACAGGAGGGATTGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCt                           <  1:829018/139‑1 (MQ=255)
                           caccacCCCGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCt                           <  2:1266433/139‑1 (MQ=255)
                                ccccGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGAt                                        <  1:1693769/121‑1 (MQ=255)
                                ccccGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGAt                                        >  2:1693769/1‑121 (MQ=255)
                                ccccGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGtt                                     <  2:21315/124‑1 (MQ=255)
                                ccccGTCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGtt                                     >  1:21315/1‑124 (MQ=255)
                                     tCAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAata                 >  1:745738/1‑139 (MQ=255)
                                      cAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGtt                                     <  1:1716415/118‑1 (MQ=255)
                                      cAGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGtt                                     >  2:1716415/1‑118 (MQ=255)
                                       aGTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATc               <  2:869507/139‑1 (MQ=255)
                                        gTAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCt              >  2:1665882/1‑139 (MQ=255)
                                         tAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGc                                              <  1:1034219/106‑1 (MQ=255)
                                         tAGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGc                                              >  2:1034219/1‑106 (MQ=255)
                                          aGCTGAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTgg            >  2:269317/1‑139 (MQ=255)
                                              gAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTAc                             <  1:1760544/118‑1 (MQ=255)
                                              gAACAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTAc                             >  2:1760544/1‑118 (MQ=255)
                                                aCAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATa      <  1:720550/139‑1 (MQ=255)
                                                aCAGGAGGGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATa      <  2:745738/139‑1 (MQ=255)
                                                 cAGGAGGGATAGCTGTTAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCTAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATaa     >  1:1373853/1‑139 (MQ=255)
                                                  aggaggGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATa      <  2:1315588/139‑1 (MQ=255)
                                                   ggaggGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACATACCTATGTCGAATATCTGGCCGATAAgg   >  2:1041606/1‑139 (MQ=255)
                                                    gaggGATAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGt  >  2:1070663/1‑139 (MQ=255)
                                                          |                                                                                                                                      
ATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCAT‑‑ACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGT  >  NC_000913/291266‑291456

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: