Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,067,300 C→T G8R (GGA→AGA)  insH1 ← IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,067,3000CT95.7% 134.9 / ‑1.9 47G8R (GGA→AGA) insH1IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (25/20);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

CTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCAC  >  NC_000913/2067222‑2067395
                                                                              |                                                                                               
cgttgcTGGTCTTACTGCTCTTACGGCTGAATTCACTTTCGTCGAAGGTAAGTTGATGACTCATGCTGAACCCTGTTCCATGGCTCAAGATGACCACCAt                                                                            >  2:598090/6‑100 (MQ=11)
                     tACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGt                                                       <  1:826872/100‑1 (MQ=33)
                                     gTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGa                                       >  2:2748725/1‑100 (MQ=1)
                                       cGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAac                                     >  1:998547/1‑100 (MQ=17)
                                        ggCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAaca                                    >  1:2535630/1‑100 (MQ=18)
                                         gCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAacac                                   >  2:1382013/1‑100 (MQ=25)
                                         gCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAacac                                   >  1:1789424/1‑100 (MQ=25)
                                          cGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAAcacc                                  <  1:2744759/100‑1 (MQ=38)
                                          cGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAAcacc                                  >  1:908431/1‑100 (MQ=25)
                                          cGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAAcacc                                  >  1:1949425/1‑100 (MQ=25)
                                          cGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGAGCTCATAGCAGGGACTTGTTCGCACCTTCCCTAGAAcacc                                  >  2:2291484/1‑100 (MQ=255)
                                           gAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAAcacca                                 <  2:2814487/100‑1 (MQ=255)
                                            aaGGTAAGTTGATGGATCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAAcaccac                                <  2:2874524/100‑1 (MQ=255)
                                            aaGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAAcaccac                                <  1:2152882/100‑1 (MQ=255)
                                            aaGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAAcaccac                                >  1:2593863/1‑100 (MQ=255)
                                                tAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAAtt                            >  1:3179028/1‑100 (MQ=255)
                                                tAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAAtt                            >  2:40598/1‑100 (MQ=255)
                                                tAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAAtt                            >  1:2886934/1‑100 (MQ=255)
                                                 aaGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAAttt                           >  2:2241078/1‑100 (MQ=255)
                                                  aGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTc                          <  2:1468411/100‑1 (MQ=255)
                                                  aGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTc                          >  2:2820026/1‑100 (MQ=255)
                                                    ttGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGc                        <  1:35242/100‑1 (MQ=255)
                                                       atgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGctct                     <  1:1536981/100‑1 (MQ=255)
                                                         gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGctctct                   <  2:1249350/100‑1 (MQ=255)
                                                          aCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCtctctc                  <  1:294227/100‑1 (MQ=255)
                                                           cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCg                 >  2:259563/1‑100 (MQ=255)
                                                           cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCg                 >  1:2073025/1‑100 (MQ=255)
                                                            tCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCgg                <  1:955363/100‑1 (MQ=255)
                                                            tCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCgg                >  2:427277/1‑100 (MQ=255)
                                                             cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGt               >  2:1219292/1‑100 (MQ=255)
                                                             cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGt               >  2:3039342/1‑100 (MQ=255)
                                                              atgatgAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTa              >  1:1214533/1‑100 (MQ=255)
                                                              atgatgAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTa              <  1:484475/100‑1 (MQ=255)
                                                               tgatgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTaa             >  2:2551109/1‑100 (MQ=255)
                                                                gatgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAg            >  1:1676846/1‑100 (MQ=255)
                                                                gatgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAg            >  2:280738/1‑100 (MQ=255)
                                                                 atgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAgg           <  1:402975/100‑1 (MQ=255)
                                                                    aaCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCAc        <  1:790508/100‑1 (MQ=255)
                                                                     aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCAc        <  2:2180717/99‑1 (MQ=255)
                                                                     aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACt       >  1:362873/1‑100 (MQ=255)
                                                                     aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACt       >  2:260893/1‑100 (MQ=255)
                                                                     aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACt       <  2:2072999/100‑1 (MQ=255)
                                                                      cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTg      <  2:1319956/100‑1 (MQ=255)
                                                                      cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTg      <  1:2352093/100‑1 (MQ=255)
                                                                       ccTGTTCTATGGCTCCAGATGACAAACATGATCTCATACCAGGGACTTGTTCGCACCTTCCCTAGAACACCGCAATTTCGCTCTCTCGGTAAGGCACt       <  2:490754/98‑1 (MQ=255)
                                                                        cTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTg      <  2:1303970/98‑1 (MQ=255)
                                                                          gTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTAGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTcac  <  2:2054596/100‑1 (MQ=255)
                                                                              |                                                                                               
CTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCAC  >  NC_000913/2067222‑2067395

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: