Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,659,576 T→C F14S (TTC→TCC)  clcB → voltage‑gated ClC‑type chloride channel ClcB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,659,5760TC96.0% 78.5 / ‑3.6 25F14S (TTC→TCC) clcBvoltage‑gated ClC‑type chloride channel ClcB
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (15/9);  total (16/9)
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.08e-01

TTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGATGCTGCTACTGGAG  >  NZ_CP009273/1659490‑1659661
                                                                                      |                                                                                     
ttttGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTCCCa                                                                                    <  1:339998/90‑1 (MQ=255)
ttttGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTCCCa                                                                                    <  1:417041/90‑1 (MQ=255)
                  gCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACTCTTCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAAc                                                                  <  2:286856/90‑1 (MQ=255)
                  gCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAAc                                                                  <  1:826125/90‑1 (MQ=255)
                      gtgtTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAgtc                                                              >  2:645264/1‑90 (MQ=255)
                      gtgtTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAgtc                                                              <  1:866805/90‑1 (MQ=255)
                      gtgtTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAgtc                                                              >  1:111643/1‑90 (MQ=255)
                                                gCATCATCTACACATCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTttgccgt                                    <  1:402454/90‑1 (MQ=255)
                                                   tcatcTACACATCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTttgccgttgc                                 <  1:779165/90‑1 (MQ=255)
                                                           acaTCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTc                         >  1:719894/1‑90 (MQ=255)
                                                           acaTCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTc                         >  1:604813/1‑90 (MQ=255)
                                                           acaTCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTc                         >  1:510078/1‑90 (MQ=255)
                                                               cTATCCGGATCTGCGCCCTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCa                     <  2:255872/90‑1 (MQ=255)
                                                               cTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCa                     >  1:467273/1‑90 (MQ=255)
                                                               cTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCa                     >  2:385552/1‑90 (MQ=255)
                                                                 aTCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATg                   >  1:617001/1‑90 (MQ=255)
                                                                      gATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGAtg              >  2:284693/1‑90 (MQ=255)
                                                                      gATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGAtg              >  2:623334/1‑90 (MQ=255)
                                                                      gATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGAtg              >  1:347159/1‑90 (MQ=255)
                                                                        tCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGAtgct            >  1:638632/1‑90 (MQ=255)
                                                                        tCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGAtgct            <  1:385552/90‑1 (MQ=255)
                                                                           gcgcACTATGTTCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGATgctgct         >  1:922498/1‑90 (MQ=255)
                                                                           gcgcACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGATgctgct         >  2:210318/1‑90 (MQ=255)
                                                                           gcgcACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGATgctgct         >  2:275373/1‑90 (MQ=255)
                                                                                  aTGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGATGCTGCTACTGGAg  >  2:891844/1‑90 (MQ=255)
                                                                                      |                                                                                     
TTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGATGCTGCTACTGGAG  >  NZ_CP009273/1659490‑1659661

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