Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,366,694 T→C T9A (ACA→GCA)  insH‑10 ← IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,366,6940TC100.0% 102.9 / NA 32T9A (ACA→GCA) insH‑10IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/22);  total (10/22)

TTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGT  >  NC_000913/3366620‑3366826
                                                                          |                                                                                                                                    
ttttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACtt                                                                     <  2:1256037/140‑1 (MQ=18)
 tttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACt                                                                      <  1:1381246/138‑1 (MQ=17)
 tttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGGTCGCACCTTCCTTAGTACt                                                                      >  2:1381246/1‑138 (MQ=17)
  ttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTcc                                                                   >  1:1186983/1‑140 (MQ=18)
  ttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTcc                                                                   >  1:746043/1‑140 (MQ=18)
    cTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTcccc                                                                 <  2:92196/140‑1 (MQ=25)
     tGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTccccc                                                                <  2:1186983/140‑1 (MQ=37)
       gTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCaa                                                              <  1:1464232/140‑1 (MQ=255)
          tGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCaataa                                                           <  2:709083/140‑1 (MQ=255)
            aCGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAAtt                                                         >  1:951438/1‑140 (MQ=255)
               gcgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGc                                                      <  2:233034/140‑1 (MQ=255)
               gcgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGc                                                      <  1:232653/140‑1 (MQ=255)
                cgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCCCCTTCCTTAGTACTTCCCCCAATAATTGGCt                                                     <  2:707694/140‑1 (MQ=255)
                     aCTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGgtgt                                                <  2:537952/140‑1 (MQ=255)
                      ctgctgAATTCACTGTCGGCGAAGGTAAGTTGGTGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGtt                                               >  1:846059/1‑140 (MQ=255)
                          tgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCaataa                                                           <  2:1216754/124‑1 (MQ=255)
                          tgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCaataa                                                           >  1:1216754/1‑124 (MQ=255)
                            aaTTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCa                                         <  2:1204098/140‑1 (MQ=255)
                              ttCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCaaa                                       <  1:435921/140‑1 (MQ=255)
                                cACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCAAAAt                                     <  2:1255473/140‑1 (MQ=255)
                                  cTGTCGGCGAAGGTAAGTTGGTGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCAAAATgg                                   <  1:1054573/140‑1 (MQ=255)
                                   tGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGt                                  >  2:742487/1‑140 (MQ=255)
                                       ggCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCaag                              <  1:305631/140‑1 (MQ=255)
                                        gCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCaaga                             >  2:953755/1‑140 (MQ=255)
                                        gCGAAGGTAAGTTGATGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCaaga                             <  1:953755/140‑1 (MQ=255)
                                            aGGTAAGTTGGTGACTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTg                         <  2:846059/140‑1 (MQ=255)
                                                     gatgaCTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGAt                <  2:746043/140‑1 (MQ=255)
                                                      atgaCTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATa               >  2:225122/1‑140 (MQ=255)
                                                          cTCATGATGAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCgg           >  1:1056822/1‑140 (MQ=255)
                                                             atgatgAACCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAACAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAACTGATAGCGGTaa        <  2:863058/140‑1 (MQ=255)
                                                                   aaCCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGt  <  2:241354/140‑1 (MQ=255)
                                                                   aaCCCTGCTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGt  <  2:231048/140‑1 (MQ=255)
                                                                          |                                                                                                                                    
TTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGT  >  NC_000913/3366620‑3366826

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