Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 354,548 T→A 8.9% P208P (CCA→CCT cbcW ← choline/betaine/carnitine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947354,5480TA 8.9% 122.1 / 3.1 45P208P (CCA→CCTcbcWcholine/betaine/carnitine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (20/21);  new base A (2/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

CCAGGCCCGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCGCAGCAATGGCGAACACCAC  >  NC_002947/354400‑354691
                                                                                                                                                       |                                                                                                                                               
ccAGGCCCGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCtt                                                                                                                                                 <  2:272164/149‑1 (MQ=255)
ccAGGCCCGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCtt                                                                                                                                                 <  7:266740/149‑1 (MQ=255)
ccAGGCCCGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCtt                                                                                                                                                 <  4:195795/149‑1 (MQ=255)
ccAGGCCCGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCtt                                                                                                                                                 <  5:24844/149‑1 (MQ=255)
ccAGGCCCGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCtt                                                                                                                                                 <  4:5465/149‑1 (MQ=255)
 cAGGCCCGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTg                                                                                                                                                <  1:74618/149‑1 (MQ=255)
 cAGGCCCGCCTCGCACCCCACGGACATATCGGCGGTGTTCCGTGCCTTGACCACAGGTGTGCCCAGGCCATCAGCGCCAACCATGGCGGCAATGCCCACCCTCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTg                                                                                                                                                <  2:119222/149‑1 (MQ=255)
      ccGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATc                                                                                                                                           <  1:120806/149‑1 (MQ=255)
      ccGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATc                                                                                                                                           >  3:181499/1‑149 (MQ=255)
             gAAGCCCAGGGGGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGc                                                                                                                                     <  5:27480/148‑1 (MQ=255)
                   cAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCTGCGATGCTTGGCATCGCGTGCGGCAGtt                                                                                                                              <  4:181499/149‑1 (MQ=255)
                             tCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACgg                                                                                                                      >  1:198308/1‑147 (MQ=255)
                             tCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACgg                                                                                                                      <  2:198308/147‑1 (MQ=255)
                                gcggTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTTGGTTCGCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGt                                                                                                                    >  1:110598/1‑148 (MQ=255)
                                      ttCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCt                                                                                                            >  1:36411/1‑148 (MQ=255)
                                        cAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTgac                                                                                                         >  1:295516/1‑149 (MQ=255)
                                                     acaGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGagcagc                                                                                             >  2:300382/1‑148 (MQ=255)
                                                         gTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCCAAGCATcgccg                                                                                                                                       >  1:2112/1‑100 (MQ=255)
                                                         gTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCCAAGCATcgccg                                                                                                                                       <  2:2112/102‑3 (MQ=255)
                                                          tttGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAgg                                                                                       >  3:282134/1‑149 (MQ=255)
                                                              cccAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGAGGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCtt                                                                                   <  2:309018/149‑1 (MQ=255)
                                                                cAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGc                                                                                 <  8:61966/149‑1 (MQ=255)
                                                                   gCCATGAGCGCCAACCAGGGCGTCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGc                                                                                      <  2:197871/141‑1 (MQ=255)
                                                                   gCCATCAGCGCCAACCAGGGCGGCAATGACGCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGCCATCGCGTCCGGCAGTTCGATACGGGTAAGCAGCTCACGACCCGAGCAGCCAAAGGc                                                                                      >  1:197871/1‑142 (MQ=255)
                                                                      aTCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGt                                                                           <  1:300382/149‑1 (MQ=255)
                                                                               aaCCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCa                                                                        <  2:297754/143‑1 (MQ=255)
                                                                               aaCCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCa                                                                        >  1:297754/1‑143 (MQ=255)
                                                                                   aGGGCGGCAATGACCACCTTCGACTGCGACAGCATGATGCTTTGAGTGACACCGGCGGCGATGCTTGGCATGGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTGTTg                                                              <  4:282134/149‑1 (MQ=255)
                                                                                    gggCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGc                                                             >  6:79821/1‑149 (MQ=255)
                                                                                          cAATAACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATTCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACTACGCGAGCAGCCAAAGGCCTTGCCGGc                                                                             <  2:218717/127‑1 (MQ=255)
                                                                                                        cgacagcgacagCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGcc                                            >  2:85781/1‑147 (MQ=255)
                                                                                                          acagcgacagcATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCagg                                        >  1:14544/1‑149 (MQ=255)
                                                                                                            agcgacagcATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCaggta                                      <  4:150196/149‑1 (MQ=255)
                                                                                                                    cATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGc                              <  3:210053/149‑1 (MQ=255)
                                                                                                                    cATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGc                              >  8:247298/1‑149 (MQ=255)
                                                                                                                        atgCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGaa                           <  2:295516/148‑1 (MQ=255)
                                                                                                                           cATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAgg                                                                                       <  6:41009/85‑1 (MQ=255)
                                                                                                                           cATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAgg                                                                                       >  5:41009/1‑85 (MQ=255)
                                                                                                                           cATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATggg                       <  2:175006/149‑1 (MQ=255)
                                                                                                                           cATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATggg                       <  2:14544/149‑1 (MQ=255)
                                                                                                                           cATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATgg                        <  1:85781/148‑1 (MQ=255)
                                                                                                                                      cacCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCGCAGCAATg             >  3:42314/1‑149 (MQ=255)
                                                                                                                                        ccGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCGCAGCAATGGc           <  5:79821/149‑1 (MQ=255)
                                                                                                                                          cggcggCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCGCAGCAATGGc           >  2:190679/1‑148 (MQ=255)
                                                                                                                                           ggcggcGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAAt                          >  8:286217/1‑132 (MQ=255)
                                                                                                                                            gcggcgATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGTAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCGCAGCAATGGCGaa        >  1:184461/1‑149 (MQ=255)
                                                                                                                                            gcggcgATGCTTGGCATCGCGTGCGGCAGTTCGAAACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAAt                          >  3:126190/1‑131 (MQ=255)
                                                                                                                                              cgcgATGCCAAGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGc                                                                                      >  1:276713/2‑69 (MQ=255)
                                                                                                                                              cgcgATGCCAAGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGc                                                                                      <  2:276713/68‑1 (MQ=255)
                                                                                                                                                  aTGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCGCAGCAATGGCGAAcaccac  >  4:257653/1‑149 (MQ=255)
                                                                                                                                                       |                                                                                                                                               
CCAGGCCCGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCGCAGCAATGGCGAACACCAC  >  NC_002947/354400‑354691

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