Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 1,660,340 C→T P257S (CCC→TCC)  clcB → H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,660,3400CT95.5% 131.1 / ‑4.4 44P257S (CCC→TCC) clcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/24);  minor base G (1/1);  total (19/25)
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.00e+00

CCACTGTTGTTAACGTTAATGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGCCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCTGCTTTCCCTGTTCACACCTGCAGTGT  >  CP009273/1660277‑1660403
                                                               |                                                               
ccACTGTTGTTAACGTTAATGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTg                                                             <  1:997219/68‑1 (MQ=255)
    tgttgtTAACGTTAATGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCaa                                                         >  1:2750759/1‑68 (MQ=255)
       tgttAACGTTAATGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTg                                                      >  1:1396344/1‑68 (MQ=255)
         ttAACGTTAATGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGc                                                    <  1:4488599/68‑1 (MQ=255)
         ttAACGTCAATGAACGCCGGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGGCCTGGCAACTGGc                                                    >  1:875787/1‑68 (MQ=255)
          tAACGTTGATGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCa                                                   >  1:869292/1‑68 (MQ=255)
          tAACGTTAATGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCa                                                   >  1:838055/1‑68 (MQ=255)
          tAACGTTAATGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCa                                                   <  1:4738832/68‑1 (MQ=255)
          tAACGTTAATGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCa                                                   >  1:3632022/1‑68 (MQ=255)
            aCGTTAATGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACt                                                 >  1:2344215/1‑68 (MQ=255)
            aCGTTAATGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACt                                                 >  1:2586228/1‑68 (MQ=255)
            aCGTTAATGAACGCCTGTCATCGTGGATTTGTGACTCTCAAACTTGCGCCGTCCTGGCAACTGGCACt                                                 >  1:2852199/1‑68 (MQ=255)
              gTTAATGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAg                                               >  1:4254424/1‑68 (MQ=255)
              gTTAATGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAg                                               <  1:2996150/68‑1 (MQ=255)
               ttACTGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAgg                                              <  1:3129118/68‑1 (MQ=255)
                 aaTGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCg                                            <  1:4518225/68‑1 (MQ=255)
                   tGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCggg                                          >  1:103918/1‑68 (MQ=255)
                   tGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCggg                                          >  1:10375/1‑68 (MQ=255)
                   tGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCggg                                          <  1:1211154/68‑1 (MQ=255)
                   tGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCggg                                          <  1:15154/68‑1 (MQ=255)
                    gAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGt                                         >  1:4876601/1‑68 (MQ=255)
                       cGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGa                                      <  1:2931583/68‑1 (MQ=255)
                       cGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGa                                      <  1:5394652/68‑1 (MQ=255)
                        gCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGAt                                     >  1:1298147/1‑68 (MQ=255)
                        gCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGAt                                     >  1:2649059/1‑68 (MQ=255)
                         ccTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGATc                                    <  1:4684979/68‑1 (MQ=255)
                            gTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGATCGTg                                 <  1:5154272/68‑1 (MQ=255)
                                tCGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTc                             >  1:1625235/1‑68 (MQ=255)
                                 cGTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCt                            <  1:2275585/68‑1 (MQ=255)
                                  gTGGATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCTg                           <  1:181248/68‑1 (MQ=255)
                                    ggATTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCGCGTTGAACGTGGGACTGca                         >  1:3026719/1‑67 (MQ=255)
                                      aTTTGTGAGTCTCAAACTTGCGCCGTCCTGGCAACAGGCACTAGGCGGGTTGATCGTGGGTCTGCttt                       >  1:261259/1‑68 (MQ=255)
                                         tgtgAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCTGCTTTccc                    >  1:5597977/1‑68 (MQ=255)
                                            gAGTCTCAAACTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCTGCTTTCCCTGt                 <  1:4090724/68‑1 (MQ=255)
                                            gAGTCTCAAACTTGCGCCGGCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCTGCTTTCCCTGt                 <  1:379148/68‑1 (MQ=255)
                                                   aaaCTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCTGCTTTCCCTGTTCACAcc          <  1:5099636/68‑1 (MQ=255)
                                                   aaaCTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCTGCTTTCCCTGTTCACAcc          <  1:4398513/68‑1 (MQ=255)
                                                   aaaCTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCTGCTTTCCCTGTTCACAcc          <  1:4713820/68‑1 (MQ=255)
                                                   aaaCTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCTGCTTTCCCTGTTCACAcc          <  1:1855567/68‑1 (MQ=255)
                                                    aaCTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCTGCTTTCCCTGTTCACACCt         <  1:4106830/68‑1 (MQ=255)
                                                    aaCTTGCGCCGTCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCTGCTTTCCCTGTTCACACCt         <  1:337026/68‑1 (MQ=255)
                                                         gcgcCGTCCTGGCAACTGGGACTAGGCGGGTTGATCGTGGGTCTGCTTTCCCTGTTCACACCTGCAgt    <  1:2412589/68‑1 (MQ=255)
                                                         gcgcCGTCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCTGCTTTCCCTGTTCACACCTGCAgt    <  1:3982153/68‑1 (MQ=255)
                                                         gcgcCGTCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGGCTGCTTTCCCTGTTCACACCTGCAgt    <  1:2877525/68‑1 (MQ=255)
                                                           gccgtcCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCTGCTTTCCCTGTTCACACCTGCAgtgt  <  1:1714596/68‑1 (MQ=255)
                                                           gccgtcCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCTGCTTTCCCTGTTCACACCTGCAgtgt  >  1:3618535/1‑68 (MQ=255)
                                                               |                                                               
CCACTGTTGTTAACGTTAATGAACGCCTGTCATCGTGGATTTGTGAGTCTCAAACTTGCGCCGCCCTGGCAACTGGCACTAGGCGGGTTGATCGTGGGTCTGCTTTCCCTGTTCACACCTGCAGTGT  >  CP009273/1660277‑1660403

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: