Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,920,303 T→G 12.5% I52L (ATC→CTC)  prc ← tail‑specific protease Prc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,920,3030TG12.5% 105.7 / 6.3 39I52L (ATC→CTC) prctail‑specific protease Prc
Reads supporting (aligned to +/- strand):  ref base T (19/15);  new base G (2/3);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TGTCGGCAGCGGTGAAGTAGCTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCATGCTGATTCGACGTGGAGGC  >  NC_002947/1920159‑1920442
                                                                                                                                                |                                                                                                                                           
tGTCGGCAGCGGTGAAGTAGCTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCCTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGTAATTCCACCACGTTGAGGCTGGCGACGATCtcg                                                                                                                                         >  4:306131/1‑149 (MQ=255)
    ggCAGCGGTGAAGTAGCTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTcgc                                                                                                                                     >  1:304835/1‑149 (MQ=255)
      cAGCGGTGAAGTAGCTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCGTGAGCAATTCCACCACGTTGAGGCTGGCGACGAGATCGTcgcca                                                                                                                                   <  8:90980/149‑3 (MQ=255)
           gTGAAGTAGCTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCggg                                                                                                                              <  4:182553/149‑1 (MQ=255)
             gAAGTAGCTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGACCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCt                                                                                                                            >  8:58379/1‑149 (MQ=255)
              aaGTAGCTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGACCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTg                                                                                                                           <  7:58379/149‑1 (MQ=255)
                gTAGCTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGTAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCa                                                                                                                         <  3:306131/149‑1 (MQ=255)
                    cTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCt                                                                                                                     <  5:183627/149‑1 (MQ=255)
                               gTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCTCACttgt                                                                                                          >  5:198498/1‑149 (MQ=255)
                                  cagcagCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTgg                                                                                                       >  6:89792/1‑149 (MQ=255)
                                   agcagcTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTggc                                                                                                      <  5:89792/149‑1 (MQ=255)
                                       gcTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCTCACTTGTTGGCGgcc                                                                                                  <  6:198498/149‑1 (MQ=255)
                                                aGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGGGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCAtt                                                                                          <  2:304835/148‑1 (MQ=255)
                                                     tCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTgccgc                                                                                     <  6:286758/148‑1 (MQ=255)
                                                                   cGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGc                                                                      >  2:156413/1‑149 (MQ=255)
                                                                    gAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTgcc                                                                                       >  4:249192/1‑131 (MQ=255)
                                                                    gAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTgcc                                                                                       <  3:249192/131‑1 (MQ=255)
                                                                      aCGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCa                                                                   <  1:156413/149‑1 (MQ=255)
                                                                           cATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGcc                                                               >  8:223918/1‑148 (MQ=255)
                                                                            aTCGTCCAGGGGCGGGTTGCTGTAGTGGTGGGGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGa                                                             <  8:87385/149‑1 (MQ=255)
                                                                            aTCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGa                                                             >  2:217797/1‑149 (MQ=255)
                                                                                   aGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGAtcatca                                                       >  2:156093/1‑148 (MQ=255)
                                                                                                tGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGGGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTg                                         >  8:96221/1‑149 (MQ=255)
                                                                                                tGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGCTCATCCCAGCCAGTGCGGTg                                         >  6:222835/1‑149 (MQ=255)
                                                                                                 gTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGc                                        >  5:66395/1‑149 (MQ=255)
                                                                                                       gtggCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTggaa                                  <  1:217797/149‑1 (MQ=255)
                                                                                                                  gCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCtt                        >  1:216260/1‑148 (MQ=255)
                                                                                                                    aaTTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCAt                     >  8:198765/1‑149 (MQ=255)
                                                                                                                      ttCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCATGc                   >  7:133357/1‑149 (MQ=255)
                                                                                                                         caccacGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCATGCTGa                >  8:59948/1‑149 (MQ=255)
                                                                                                                         caccacGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCATGCTGa                >  1:80924/1‑149 (MQ=255)
                                                                                                                          accacGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCATGCTGAt               >  2:257245/1‑149 (MQ=255)
                                                                                                                                  gAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCATGCTGATTCGACGt        <  2:216260/148‑1 (MQ=255)
                                                                                                                                     gcgggCGACGTGCTCGTGGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATGGCCGCCGAGCGTGAGGCTGCCCAGTCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCATGCTGATTCGACGTGGAg    <  7:169715/146‑1 (MQ=255)
                                                                                                                                     gCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGTTTCATGCTGATTCGACGTGGAg    <  7:198765/149‑1 (MQ=255)
                                                                                                                                      accgCGACGAGATCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTggaa                                  <  6:289061/115‑1 (MQ=255)
                                                                                                                                      accgCGACGAGATCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTggaa                                  >  5:289061/4‑118 (MQ=255)
                                                                                                                                      accgCGACGAGATCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCATGCTGATTCGACGTGGAgg   >  6:280888/4‑149 (MQ=255)
                                                                                                                                         gCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCATGCTGATTCGACGTGGAGGc  <  6:66395/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                           
TGTCGGCAGCGGTGAAGTAGCTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCATGCTGATTCGACGTGGAGGC  >  NC_002947/1920159‑1920442

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: ‑