Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,467,800 A→T 9.3% S15S (TCA→TCT)  PP_3959 → voltage‑gated chloride channel family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,467,8000AT 9.3% 106.6 / 3.7 43S15S (TCA→TCT) PP_3959voltage‑gated chloride channel family protein
Reads supporting (aligned to +/- strand):  ref base A (21/18);  new base T (2/2);  total (23/20)
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.79e-01

GCAAGCGGATCGCCACTCCCGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCACCCTGGTCACCTCTTCTC  >  NC_002947/4467652‑4467947
                                                                                                                                                    |                                                                                                                                                   
gCAAGCGGATCGCCACTCCCGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCa                                                                                                                                                     >  7:81503/1‑149 (MQ=255)
      ggATCGCCACTCCCGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTgg                                                                                                                                               <  6:200155/149‑1 (MQ=255)
        aTCGCCACTCCCGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCGAAGTCAGCCATGTGTGGACAATGTCCGGCCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGc                                                                                                                                              >  7:112423/1‑148 (MQ=255)
              aCTCCCGGTGCGGCTACAACCTGTGTCGTGGAGGCCTGGGCGGGAAAACAAGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATg                                                                                                                                       <  4:180028/149‑1 (MQ=255)
                  ccGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCATTCTCGCCAGCGTGAACGGCgg                                                                                                                                   <  5:141626/149‑2 (MQ=255)
                     gTGCGGCTACAACCTGTTTCGGGTAGGCCTGTGCTGGAAAACACGCCATACTCACGGATTTCCAAGTCAGCCATGTGTTGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAg                                                                                                                                <  3:166839/149‑1 (MQ=255)
                      tGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGc                                                                                                                               <  8:81503/149‑1 (MQ=255)
                            tACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGc                                                                                                                                              <  1:111151/128‑1 (MQ=255)
                            tACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGc                                                                                                                                              >  2:111151/1‑128 (MQ=255)
                              cAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCg                                                                                                                       >  4:197957/1‑149 (MQ=255)
                               aaCCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGt                                                                                                                      >  4:226777/1‑149 (MQ=255)
                                    gtgtCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTg                                                                                                                  >  7:25739/1‑148 (MQ=255)
                                         gggTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCgg                                                                                                            <  3:197957/149‑1 (MQ=255)
                                           gTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCgggg                                                                                                          <  1:1550/149‑1 (MQ=255)
                                                         gAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTgggg                                                                                            >  2:70746/1‑149 (MQ=255)
                                                             acacGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCg                                                                                        >  3:220730/1‑149 (MQ=255)
                                                               acGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATAGGGGCATTGCTGGTGGGGCTCGtt                                                                                      >  5:255582/1‑149 (MQ=255)
                                                                 gCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTgc                                                                                    >  6:212518/1‑149 (MQ=255)
                                                                    aTACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGTCGTCAGCGACTGGCGTTCTGGATCGTGGCATTGCTGGTGGGGCTCGTTGCGCt                                                                                 <  8:25739/149‑1 (MQ=255)
                                                                        tCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCt                                                                             >  4:125157/1‑149 (MQ=255)
                                                                              cTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGcc                                                                       <  1:70746/149‑1 (MQ=255)
                                                                                   cAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCATTCTCGCCAGCGTGAACGGCgg                                                                                                                                   >  3:208262/1‑83 (MQ=37)
                                                                                   cAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCATTCTCGCCAGCGTGAACGGCgg                                                                                                                                   <  4:208262/84‑2 (MQ=37)
                                                                                    aaGTCAGCCATGTGTGGAGAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTGCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCt                                                                      <  6:57226/144‑1 (MQ=255)
                                                                                    aaGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCACGTTGCGCTGGCTTTCGCCTgagaa                                                                 >  5:57226/1‑145 (MQ=255)
                                                                                                  tgGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGtt                                                   >  1:95140/1‑149 (MQ=255)
                                                                                                    gACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCAGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCTCCTGGCTAGCCGACGCGGCGTTCg                                                 >  8:128524/1‑149 (MQ=255)
                                                                                                       aaTGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCtt                                              >  1:163158/1‑149 (MQ=255)
                                                                                                       aaTGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCtt                                              >  7:38924/1‑149 (MQ=255)
                                                                                                              gACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCaa                                       <  4:185921/149‑1 (MQ=255)
                                                                                                              gACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCaa                                       <  3:226777/149‑1 (MQ=255)
                                                                                                                 ccTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCaa                                       <  6:96129/146‑1 (MQ=255)
                                                                                                                 ccTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCaa                                       >  5:96129/1‑146 (MQ=255)
                                                                                                                 ccTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAgcg                                    >  6:1017/1‑149 (MQ=255)
                                                                                                                        gACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTg                             <  2:95140/149‑1 (MQ=255)
                                                                                                                            tGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAAt                         >  6:169351/1‑149 (MQ=255)
                                                                                                                                 cTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAgcg                                    <  7:219604/133‑1 (MQ=255)
                                                                                                                                 cTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAgcg                                    >  8:219604/1‑133 (MQ=255)
                                                                                                                                    cACTCCCCCGCCGTTCACGCTGGCTAGAATGGCGTCATCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCAc                   <  6:255582/147‑1 (MQ=255)
                                                                                                                                       tCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCACCCTgg              >  6:256434/1‑149 (MQ=255)
                                                                                                                                        ctgacgccaTTCTCGCCAGCGTGAACGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCACCCTGGt             >  7:221877/5‑149 (MQ=255)
                                                                                                                                                 ttCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCACCCTGGTCACCtcttt    >  2:241432/1‑148 (MQ=255)
                                                                                                                                                 ttCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCACCCTGGTCACCtcttc    <  7:128524/149‑1 (MQ=255)
                                                                                                                                                   cACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCACCCTGGTCACCtctttt   <  1:241432/148‑3 (MQ=255)
                                                                                                                                                    aCGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCACCCTGGTCACCTCTtctc  <  8:38924/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCAAGCGGATCGCCACTCCCGGTGCGGCTACAACCTGTGTCGGGTAGGCCTGTGCGGGAAAACACGCCATACTCACGGCTTTCCAAGTCAGCCATGTGTGGACAATGTCCGACCCTACCCGACCTGCACCTTCACTCCCCCGCCGTTCACGCTGGCGAGAATGGCGTCAGCGACTGGCGTTCTGGATCGGGGCATTGCTGGTGGGGCTCGTTGCGCTGGCTTTCGCCTGGCTAGCCGACGCGGCGTTCGCTTCGTTCAAGCGCTTGCTGGAATGGGCACCCTGGTCACCTCTTCTC  >  NC_002947/4467652‑4467947

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