Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,475,902 G→A 8.3% T11T (ACC→ACT)  PP_3965 ← transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,475,9020GA 8.3% 123.3 / 2.4 48T11T (ACC→ACT) PP_3965transposase
Reads supporting (aligned to +/- strand):  ref base G (25/19);  new base A (2/2);  total (27/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 = 1.13e-01

CACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGC  >  NC_002947/4475757‑4476011
                                                                                                                                                 |                                                                                                             
cacCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTg                                                                                                            <  3:6036/149‑1 (MQ=255)
 acCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTgg                                                                                                           >  6:144677/1‑149 (MQ=255)
 acCGACGCCATCGTGCACCAGCACTTTCATACGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCCCCCGGGCAAATGCGGTCTCCGTGCCACCGCGCATGTCCATCGGCTCGGTgg                                                                                                           >  3:257594/1‑149 (MQ=14)
    gACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGgcca                                                                                                        <  5:144677/149‑1 (MQ=255)
         cATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGgccagc                                                                                                      >  1:225813/1‑146 (MQ=255)
         cATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGgccagc                                                                                                      >  8:183945/1‑146 (MQ=255)
         cATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGccagccag                                                                                                   <  3:105461/149‑1 (MQ=255)
                cACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCAt                                                                                            <  6:163877/149‑1 (MQ=255)
                cACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCAt                                                                                            <  2:11507/149‑1 (MQ=255)
                        cTTTCATCCGGTTAGACCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCgg                                                                                    <  4:277790/149‑1 (MQ=255)
                         tttCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGa                                                                                   >  5:10206/1‑149 (MQ=39)
                           tCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGa                                                                                   <  2:225813/147‑1 (MQ=39)
                              tCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCAtt                                                                              <  6:10206/149‑1 (MQ=38)
                                       cccgggtgattgcggtgtttggtgcggcgaagccgcactgcgcttatctgtttgccaatcgccgggccaaTGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAg                                                                     >  3:299271‑M2/71‑149 (MQ=255)
                                       cccgggtgattgcggtgtttggtgcggcgaagccgcactgcgcttatctgtttgccaatcgccgggccaaTGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAg                                                                     >  5:217120‑M2/71‑149 (MQ=255)
                                        CCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGc                                                                                      <  8:26340‑M1/131‑1 (MQ=255)
                                        CCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGc                                                                                      >  7:26340‑M1/1‑131 (MQ=255)
                                        CCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGc                                                                    >  7:183945‑M1/1‑149 (MQ=255)
                                          GGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGccc                                                                  >  5:64744‑M1/1‑149 (MQ=255)
                                                tGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATa                                                             <  8:201878/148‑1 (MQ=37)
                                                   cggtgttcggtgcggcgaagccgcactgcgcttatctgttcgccaatcgccgggctaaTGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAAtt                                                         >  1:101264‑M2/59‑149 (MQ=255)
                                                     aaCAGATAAGCGCAGTGCGGCTTCGCCGCACCGCACACCGCAATCACCCGGGCTAATGCTGTCTAGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCCGCCAGATGGCCTCGCTGCGGATCATTGAGCGAAAGCCCGAATAAATTgc                                                       >  8:259173/1‑149 (MQ=11)
                                                      tgtttggtgcggcgaagccgcactgcgcttatctgtttgccaatcgccgggccaaTGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTgcg                                                      <  4:299271‑M2/94‑1 (MQ=255)
                                                        aGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGca                                                    <  7:259173/149‑1 (MQ=34)
                                                         gATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGca                                                    >  6:335809/1‑148 (MQ=34)
                                                          aTAAGCGCAGTGCGGCTTAGCCGCACCGAAAACCGCCAACACCCGGGCTAATGCTGTCTGGGTGCCACCGCCCATGTCCCTAGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATAATTGAGCGAAATCCCCAATAAATTTCGcaa                                                   >  6:227124/1‑147 (MQ=11)
                                                           tAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCAGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACaa                                                 >  5:318466/1‑149 (MQ=21)
                                                           tAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCAGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACaa                                                 >  5:322927/1‑149 (MQ=21)
                                                                 cAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCg                                           <  4:204116/149‑1 (MQ=34)
                                                                        gCTTCGCCGCACCCAACACCGCAATCACCCGGCCTAATGCTGTCTCGGAGCCAGCGCCCATGTCCATCGGCTCGGTGGCCAGCCAGAAGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATcc                                     >  6:113876/1‑148 (MQ=11)
                                                                                  aCCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCGTGTACATAGGCTCGGTCGCCAGCCACATCCCAACGATCCGCATCATACAGCGAAATCCCGAATAAATTGCGCACAACCCTCCGAATCCGACGCTGGcc                           >  2:150818/1‑148 (MQ=11)
                                                                                  aCCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCa                          >  2:76209/1‑149 (MQ=34)
                                                                                    cGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCa                          >  5:152919/1‑147 (MQ=33)
                                                                                    cGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCa                          <  6:152919/147‑1 (MQ=33)
                                                                                     gAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCAGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCAttt                       <  6:322927/149‑1 (MQ=17)
                                                                                     gAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCAGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCAttt                       <  6:318466/149‑1 (MQ=17)
                                                                                      aaCACCGCAATCAGCCGGGCTAATGGTGTCTCGGTGCCAGGGCGCGTGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGATAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTc                      <  1:150818/149‑1 (MQ=11)
                                                                                           cGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGt                 <  5:335809/149‑1 (MQ=32)
                                                                                            gccaatcgccgggccaaTGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGt                 <  6:217120‑M2/131‑1 (MQ=255)
                                                                                              aaTCAGCCGTGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATTCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGAt              <  8:127521/149‑1 (MQ=12)
                                                                                                  aCCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACt          >  3:59271/1‑149 (MQ=32)
                                                                                                   cccGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCatc          >  4:196291/1‑146 (MQ=17)
                                                                                                         ctaaTGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCgg   >  7:61412‑M2/5‑149 (MQ=255)
                                                                                                          taaTGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGTGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGc  >  5:293669‑M2/4‑149 (MQ=255)
                                                                                                          taaTGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGc  >  5:278470‑M2/4‑149 (MQ=255)
                                                                                                           aaTGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCccgtc  >  5:181049‑M2/3‑146 (MQ=255)
                                                                                                           aaTGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGc  >  4:87817‑M2/3‑148 (MQ=255)
                                                                                                           aaTGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGc  <  3:257401‑M2/146‑1 (MQ=255)
                                                                                                                                                 |                                                                                                             
CACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGGGCTAATGCTGTCTCGGTGCCAGCGCGCATGTCCATCGGCTCGGTGGCCAGCCAGATGGCATCGATGCGGATCATTGAGCGAAAGCCCGAATAAATTGCGCACAACCCTCGGGATCCGAGGCTGGCCATTTCACCGTGATCACTTGCCCGGC  >  NC_002947/4475757‑4476011

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