Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,920,304 C→A 10.3% E51D (GAG→GAT)  prc ← tail‑specific protease Prc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,920,3040CA10.3% 115.5 / 4.5 39E51D (GAG→GAT) prctail‑specific protease Prc
Reads supporting (aligned to +/- strand):  ref base C (19/16);  new base A (2/2);  total (21/18)
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.64e-01

TGTCGGCAGCGGTGAAGTAGCTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCATGCTGATTCGACGTGGAGGCAATTGATCAC  >  NC_002947/1920159‑1920452
                                                                                                                                                 |                                                                                                                                                    
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                                                                                                 <  3:249192/131‑1 (MQ=255)
                                                                    gAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTgcc                                                                                                 >  4:249192/1‑131 (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                                            >  5:289061/4‑118 (MQ=255)
                                                                                                                                      accgCGACGAGATCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTggaa                                            <  6:289061/115‑1 (MQ=255)
                                                                                                                                      accgCGACGAGATCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCATGCTGATTCGACGTGGAgg             >  6:280888/4‑149 (MQ=255)
                                                                                                                                         gCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCATGCTGATTCGACGTGGAGGc            <  6:66395/147‑1 (MQ=255)
                                                                                                                                                 ctcGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCATGCTGATTCGACGTGGAGGCAATTGATCAc  >  8:66660/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGTCGGCAGCGGTGAAGTAGCTGCGCGCCGGGTCCAGCAGCTTGATGTAGCTGTCGTAGATAATGACCGAACGGGCATCGTCCAGGGGCGGCTTGCTGTAGTGGTGGCGCTTGAGCAATTCCACCACGTTGAGGCTGGCGACGATCTCGTCGCGGTCGGGCTGCAGGCTGTCCCACTTGTTGGCGGCCGCCGCATTGCCGCCGAGCGTGAGGCTGCCCAGGCCGATCATCACAGCCAGTGCGGTGCTTGGAAGGAAATGCTTCATGCTGATTCGACGTGGAGGCAATTGATCAC  >  NC_002947/1920159‑1920452

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