Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,666,468 A→T 6.5% Q385L (CAG→CTG)  clcB → putative chloride:H(+) antiporter ClcB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,666,4680AT 6.5% 163.3 / 2.5 60Q385L (CAG→CTG) clcBputative chloride:H(+) antiporter ClcB
Reads supporting (aligned to +/- strand):  ref base A (23/32);  new base T (2/2);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

CGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTT  >  NC_000913/1666326‑1666600
                                                                                                                                              |                                                                                                                                     
cggaGCCTGGCATTATGGTTCCCTGATGGCGACGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCa                                                                                                                                       <  2:738618/140‑1 (MQ=255)
    gCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTa                                                                                                                                   <  1:517128/143‑1 (MQ=255)
       tGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGGTGATGCTATGTCAAGCGACCGGGGCGTCTCCGCCAccc                                                                                                                                >  2:1007370/1‑141 (MQ=255)
       tGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACGGTTGGCGGCGACCGCGCACGCGCCGATTATGTCGACGTTGATGATATGGGAAAGGACCGGGGGGTATCAGCGACTc                                                                                                                                >  1:218375/1‑143 (MQ=255)
          gATTTTGGTGCCCGGGTGGCGAAGAAATTACACTTTTACTCGGATGGACCGGGATGGCGACACTGTTGGGGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTcccc                                                                                                                             <  2:1013719/143‑1 (MQ=255)
           aTTATGGTGCCCTGATGGCGAAGAAATTACACTTTTCCTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCg                                                                                                                            <  1:683338/143‑1 (MQ=255)
                gtttCCCTGATGGCGAAGAATTTACACTTTTACTCGGATTGACCGGGAGGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCGGtcat                                                                                                                      <  2:114220/141‑4 (MQ=255)
                 gTTCCCTGAGGGCGAAGAACTTACACTTTTACTCGGATGGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTat                                                                                                                     <  2:237246/143‑1 (MQ=255)
                  ttCCCTGATGGCGAAGAAATTACACTTTTACTCGGATAGACCGGAGTGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTGGATGGTATGTGAAATGACCGGGGAGTGTCTGCTGCTCCGCGGgttgtt                                                                                                                    >  1:546967/1‑143 (MQ=255)
                   tCCCTGAGGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattg                                                                                                                   <  1:1016581/143‑1 (MQ=255)
                     tctcaTGGCGAAGAAACTACCCTTTTACTCGTATTGACCGGGATGGCGACACTGTGGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgat                                                                                                                 <  1:771363/139‑1 (MQ=255)
                     ccTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTattgat                                                                                                                 <  2:893248/143‑1 (MQ=255)
                        gATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTCCTCCCCGGTTTATTGATTGc                                                                                                              <  1:704638/143‑1 (MQ=255)
                        gATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGc                                                                                                              >  1:760168/1‑143 (MQ=255)
                          gggCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTTCCt                                                                                                            <  1:781757/142‑1 (MQ=255)
                              gaagaaATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCg                                                                                                                            >  1:735371/1‑124 (MQ=255)
                              gaagaaATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCg                                                                                                                            <  2:735371/124‑1 (MQ=255)
                              gaagaaATTACACTTTTACTCGGATTGACAGGGATGGCGACACTGTTGGCGGTGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGCGAGTATCAGCTACTCCCCGGCTTATTGCCTGCTTGCgc                                                                                                        >  2:784003/1‑142 (MQ=255)
                                     ttACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTAGCGGCGACCACGCACGCGCCGCCCATGTGGACGGTGATGCACACTGAAATGACCGGGGAGTCACAGCTACACCCCCGGTTCTTGATTGCCCGCGTaaccgcg                                                                                                 >  2:921655/1‑138 (MQ=255)
                                                  cGGATTGGACGGGATGGCGACACTGTTGGCGCCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgc                                                                                    <  1:648514/143‑1 (MQ=255)
                                                  cGGATTGACCGGGCTGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgc                                                                                    <  1:987288/143‑1 (MQ=255)
                                                   ggctggcccGCGATGGCGACACGGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                                                   <  1:162281/136‑1 (MQ=255)
                                                   ggATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTGATTTcgcg                                                                                   >  1:301546/1‑143 (MQ=255)
                                                   ggATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGAGGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTcgcg                                                                                   <  2:468113/143‑1 (MQ=255)
                                                     aTTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCAAGCGCTGATTATGTCGACGTTGCTTATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGa                                                                                 <  1:663034/143‑1 (MQ=255)
                                                         aCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGAGGtt                                                                             >  1:747290/1‑143 (MQ=255)
                                                          ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTa                                                                            <  1:168904/143‑1 (MQ=255)
                                                                tGGCGACACAGGGGGCGGCGAACACGCACGCGCCGATTATGACGACGTTGATGATAGGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGt                                                                      <  1:298484/143‑1 (MQ=255)
                                                                    gACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGAct                                                                  >  1:788405/1‑143 (MQ=255)
                                                                    gACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATCTCGCGGACGTTACACCGTGAct                                                                  >  2:188225/1‑143 (MQ=255)
                                                                      cacTGGAGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGActct                                                                <  1:982345/143‑1 (MQ=255)
                                                                       acTGTTGGCGGCGACCGCGCACGCGCCGATTATGTCGTCGTTGATGATATGTGAAATGACCGGGGAGTGTCAGCTACTCCCCGGTTTATTCATTGCCTGCGTAATTGCGTCGGTAATTTCGCCGACGTTACACCGTGACtcta                                                               <  2:408903/143‑1 (MQ=255)
                                                                           ttGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCGCCGGTTGATTGATTGCATGGGTAATTGCGTAGGTAATTTCGGGGAAGTTACAACGTGAAtctatagg                                                           >  2:584482/1‑140 (MQ=255)
                                                                            tGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGTTTGCCTGCGAAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatctgc                                                          >  1:288089/1‑141 (MQ=255)
                                                                               cggcgACCACGCACGCGCCGATGATGTCTGCGGTGATGAGATGTGAGATGGTCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGc                                                       <  1:1021230/143‑1 (MQ=255)
                                                                                ggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTATAAGCGGAGATGATTATTGGCAAACAAATGGCGGCGGTAATTCCGCGGCAGTTACCCCGTGACTCTCTCTACCGcc                                                      >  2:1034923/1‑143 (MQ=255)
                                                                                  cgACCGCGCACGCGCGGATTATGTGGACGTTGCTGATATGTGAAATGACGGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAg                                                    <  2:1056781/143‑1 (MQ=255)
                                                                                  cgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGCCATTTCGCGGACGTTCCACCGTGACTCCATCTACCGccgc                                                    >  2:996183/1‑141 (MQ=255)
                                                                                    aCCACGCAGGCGCCGATTATGTCGACGTTGAGTGTATGTGAGCTGTCCGCTGAGTATCAGCTACTCCCCGGCTGATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGca                                                  <  2:916337/143‑1 (MQ=255)
                                                                                      cacgcacgCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACAGCGTGACTCTATCTACCGCCAGcaca                                                >  1:595128/1‑143 (MQ=255)
                                                                                      cacgcacgCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGcaca                                                <  1:709071/143‑1 (MQ=255)
                                                                                         gcacgcGCCGATTATGTCGACGTTGATAATATGTGAAATAACCGGGGAGTATCAGCTACTCCCCGGGTTAGTGATTGCCTGCGTAATTGCGTCGGTAAGTTCGCGGACGTTACACCGTGACTCTATCTACTGCCAgaaacacg                                             >  2:391741/1‑136 (MQ=255)
                                                                                           ccgcgCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgcg                                           <  1:707531/142‑1 (MQ=255)
                                                                                              cgcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATGGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag                                        <  1:377742/143‑1 (MQ=255)
                                                                                                ccGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTGCACCGTGACTCTATCTACCGCCAGCACACTGCgcagca                                      >  1:660919/1‑143 (MQ=255)
                                                                                                    ttATGTCGACGTTGATGATATGTTAAAGGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGc                                  <  2:28340/143‑1 (MQ=255)
                                                                                                        gTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAGTTGCGCGGCCGTTACACCGTGACTCTATCCACCGCCAGCGCACTTGTCAGCACAGCtggc                              >  2:852211/1‑140 (MQ=255)
                                                                                                        gTCGACCTTGATGATACGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAg                              <  1:67534/143‑1 (MQ=255)
                                                                                                         tCGACGGTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGc                             <  1:314691/143‑1 (MQ=255)
                                                                                                              gTTGATGATCTGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGCCGTTACACCGTGACTCTATCTACCGCCAGCACACTGCCCAGAATACCCAAGCGGAta                        >  2:1021318/1‑142 (MQ=255)
                                                                                                              gTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCGACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATg                        >  1:442609/1‑143 (MQ=255)
                                                                                                                 gatgatATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTAt                     >  1:613642/1‑143 (MQ=255)
                                                                                                                        tgtgCAATGACCGGGGAGTATCAGCTACTCCACGGTTTATTGATTGCCTGCGTACTTGCGTCGGTACTTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTTCGCAGCATAGCTAAGCGGATGTATTGCCCCa              <  2:862959/143‑1 (MQ=255)
                                                                                                                              ctgcccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCAc        <  2:779761/139‑1 (MQ=255)
                                                                                                                               tGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACg       <  1:777237/143‑1 (MQ=255)
                                                                                                                                  ccGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTg    <  1:230935/143‑1 (MQ=255)
                                                                                                                                   cGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTACTTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCCTTGCTACGCGGATGTATTGCCCCACTTCACGCTGt   >  2:242828/1‑143 (MQ=255)
                                                                                                                                   cGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACg       >  1:41208/1‑139 (MQ=255)
                                                                                                                                   cGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACg       <  2:41208/139‑1 (MQ=255)
                                                                                                                                   cGGGGAGTAGCTGATACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGt   >  1:846738/1‑143 (MQ=255)
                                                                                                                                    ggggAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCAAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCCCGCTGtt  >  1:55979/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                     
CGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGGTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTT  >  NC_000913/1666326‑1666600

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