Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 279,852 A→G C8C (TGT→TGC insA ← IS1 repressor TnpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913279,8520AG100.0% 108.2 / NA 30C8C (TGT→TGCinsAIS1 repressor TnpA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/20);  total (10/20)

TAGCGCTGATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAATATACTCAAC  >  NC_000913/279787‑279997
                                                                 |                                                                                                                                                 
tAGCGCTGATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAac                                                               >  1:137461/1‑150 (MQ=17)
    gcTGATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACaa                                                           >  2:545178/1‑150 (MQ=25)
     cTGATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAAt                                                         >  1:876428/1‑151 (MQ=37)
       gATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAAt                                                         <  1:344026/149‑1 (MQ=37)
       gATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAAt                                                         <  1:945038/149‑1 (MQ=37)
       gATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATgg                                                       <  2:137459/151‑1 (MQ=38)
        aTGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATggg                                                      <  1:290782/151‑1 (MQ=39)
        aTGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATggg                                                      <  2:282312/151‑1 (MQ=39)
        aTGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATggg                                                      <  1:675833/151‑1 (MQ=39)
         tGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGc                                                     >  2:367967/1‑151 (MQ=255)
          gTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGc                                                     <  1:930618/150‑1 (MQ=255)
          gTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCt                                                    <  1:765198/151‑1 (MQ=255)
               gCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTc                                               <  2:425136/151‑1 (MQ=255)
                cAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCg                                              >  1:277843/1‑151 (MQ=255)
                  gTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGc                                            >  1:648234/1‑151 (MQ=255)
                             gTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCt                                 <  1:709329/151‑1 (MQ=255)
                                 cGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTaaa                              >  2:926907/1‑150 (MQ=255)
                                 cGCACCACCCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTaaa                              <  2:758421/150‑1 (MQ=255)
                                    accaccCCGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGa                          <  2:60556/151‑1 (MQ=255)
                                         cccGTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGt                     <  1:792589/151‑1 (MQ=255)
                                            gTCAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCtt                  <  1:730204/151‑1 (MQ=255)
                                              cAGTAGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGc                >  2:964006/1‑151 (MQ=255)
                                                  aGCTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCaa            <  1:1045928/151‑1 (MQ=255)
                                                    cTGAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAat           >  2:639479/1‑150 (MQ=255)
                                                      gAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAata          <  2:345222/149‑1 (MQ=255)
                                                      gAACAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAata          <  2:725316/149‑1 (MQ=255)
                                                        aCAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAATATACt      >  1:308111/1‑151 (MQ=255)
                                                         cAGGAGGGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAATATACTc     <  1:232561/151‑1 (MQ=255)
                                                          aggaggGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAATATACTCa    <  2:515035/151‑1 (MQ=255)
                                                            gaggGGCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAATATACTCAAc  <  1:835969/151‑1 (MQ=255)
                                                                 |                                                                                                                                                 
TAGCGCTGATGTCCGGCAGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAATATACTCAAC  >  NC_000913/279787‑279997

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