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,5760TC87.5% 53.1 / 1.9 24F14S (TTC→TCC) clcBvoltage‑gated ClC‑type chloride channel ClcB
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (11/10);  total (12/12)
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.27e-01

TTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGATGCTGCTACTGG  >  NZ_CP009273/1659490‑1659659
                                                                                      |                                                                                   
ttttGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTCCCa                                                                                  <  2:141928/90‑1 (MQ=255)
                      gtgtTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAgtc                                                            <  2:230073/90‑1 (MQ=255)
                      gtgtTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAgtc                                                            >  1:565297/1‑90 (MQ=255)
                         tttCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAgtcgtc                                                         <  2:784609/90‑1 (MQ=255)
                             aTTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTa                                                     >  1:393961/1‑90 (MQ=255)
                                          gcacATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCtt                                        >  2:771603/4‑90 (MQ=255)
                                             cATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTttgc                                     >  1:195451/1‑90 (MQ=255)
                                               tGCATCATCTACACATCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTttgccg                                   >  1:537794/1‑90 (MQ=255)
                                                         acacATCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGtt                         >  2:382660/1‑90 (MQ=255)
                                                          cacaTCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGttt                        <  2:519084/90‑1 (MQ=255)
                                                             aTCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGt                     >  1:297435/1‑90 (MQ=255)
                                                             aTCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGt                     <  2:537794/90‑1 (MQ=255)
                                                             aTCTATCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGt                     <  2:264084/90‑1 (MQ=255)
                                                                 aTCCGGATCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATg                 >  2:747778/1‑90 (MQ=255)
                                                                        tCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGAtgct          <  1:125945/90‑1 (MQ=255)
                                                                        tCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGAtgct          <  2:434934/90‑1 (MQ=255)
                                                                        tCTGCGCACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGAtgct          <  2:796782/90‑1 (MQ=255)
                                                                           gcgcACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGATgctgct       >  2:272802/1‑90 (MQ=255)
                                                                           gcgcACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGATgctgct       >  2:696362/1‑90 (MQ=255)
                                                                           gcgcACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGATgctgct       >  1:350637/1‑90 (MQ=255)
                                                                           gcgcACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGATgctgct       >  1:120891/1‑90 (MQ=255)
                                                                              cACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGATGCTGCTACt    <  1:794870/90‑1 (MQ=255)
                                                                              cACTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGATGCTGCTACt    <  1:45140/90‑1 (MQ=255)
                                                                                cTATGTCCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGATGCTGCTACTgg  <  1:724679/90‑1 (MQ=255)
                                                                                      |                                                                                   
TTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTCATGCGATGCTGCTACTGG  >  NZ_CP009273/1659490‑1659659

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