Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,010,385 C→T 11.7% intergenic (‑1291/‑377) PP_2628 ← / → PP_2631 ABC transporter ATP‑binding protein/cellulose biosynthesis protein BcsF/YhjT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,010,3850CT11.7% 83.6 / 3.0 34intergenic (‑1291/‑377)PP_2628/PP_2631ABC transporter ATP‑binding protein/cellulose biosynthesis protein BcsF/YhjT
Reads supporting (aligned to +/- strand):  ref base C (15/15);  new base T (2/2);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

TACGTCTGGCGTCGACACCTGCCAACCGACTTCGACGCCTTGCTGGCACGCCTGCGGCCGCTGGCCATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCT  >  NC_002947/3010256‑3010532
                                                                                                                                 |                                                                                                                                                   
ctcGTCTGGCGTCGACACCTGCCAAGCGACTTCGACGCCTTGCTGGCACGCCTGCGGCCGCTGGCCATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTcc                                                                                                                                   <  5:249900/146‑1 (MQ=255)
  cGTCTGGCGTCGACACCTGCCAACCGACTTCGACGCCTTGCTGGCACGCCTGCGGCCGCTGGCCATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTccac                                                                                                                                 <  7:166671/148‑1 (MQ=255)
   gTCTGGCGTCGACACCTGCCAACCGACTTCGACGCCTTGCTGGCACGCCTGCGGCCGCTGGCCATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGACGAAATCGTccacca                                                                                                                               <  8:141680/149‑1 (MQ=255)
   gTCTGGCGTCGACACCTGCCAACCGACTTCGACGCCTTGCTGGCACGCCTGCGGCCGCTGGCCATCTGCGGGTTGGTAGCGCCACGGGCGTTAGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTcc                                                                                                                                   >  6:249900/1‑145 (MQ=255)
     cTGGCGTCGACACCTGCCAACCGACTTCGACGCCTTGCTGGCACGCCTGCGGCCGCTGGCCATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGt                                                                                                                             >  5:273638/1‑149 (MQ=255)
        gCGTCGACACCTGCCAACCGACTTGGACGCATTGCTGGCAGGCGTGCGGCGGCGGGCATTCTGGGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCATAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGc                                                                                                                          <  1:259827/149‑1 (MQ=255)
        gCGTCGACACCTGCCAACCGACTTCGACGCCTTGCTGGCACGCCTGCGGCCGCTGGCCATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGtt                                                                                                                            <  7:332843/147‑1 (MQ=255)
        gCGTCGACACCTGCCAACCGACTTCGACGCCTTGCTGGCACGCCTGCGGCCGCTGGCCATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGc                                                                                                                          <  6:273638/149‑1 (MQ=255)
          atCGACACCTGCCAACCGACTTCGACGCCTTGCTGGCACGCCTGCGGCCGCTGGCCATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTg                                                                                                                        >  7:67892/2‑149 (MQ=255)
                  cTGCCAACCGACTTCGACGCCTTGCTGGCACGCCTGCGGCCGCTGGCCATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAAACGTCCACCAGTTGCTGGTGTTGCg                                                                                                                >  6:129102/1‑149 (MQ=255)
                   tGCCAACCGACTTCGACGCCTTGCTGGCACGCCTGCGGCCGCTGGCCATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGt                                                                                                               >  8:173098/1‑149 (MQ=255)
                                             gCACGCCTGCGGCCGCTGGCCATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTg                                                                                       >  4:10422/1‑147 (MQ=255)
                                             gCACGCCTGCGGCCGCTGGCCATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGATGCGTCCGGCGCCTGGCCTGACCCCGCTg                                                                                       >  7:160210/1‑147 (MQ=255)
                                                  ccTGCGGCCGCTGGCCATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCtt                                                                                <  7:173098/149‑1 (MQ=255)
                                                                 cATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGTGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCgcctggc                                                                                                   >  1:299565/1‑115 (MQ=255)
                                                                 cATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGTGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCgcctggc                                                                                                   <  6:122685/115‑1 (MQ=255)
                                                                 cATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGTGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCgcctggc                                                                                                   >  5:122685/1‑115 (MQ=255)
                                                                 cATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGTGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCgcctggc                                                                                                   <  2:299565/115‑1 (MQ=255)
                                                                          gTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTg                                                        >  7:331902/1‑149 (MQ=255)
                                                                                         ggCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGCTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCt                                         <  8:160210/149‑1 (MQ=255)
                                                                                          gCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCTCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACAtt                                             <  8:48653/144‑1 (MQ=255)
                                                                                          gCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCTCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACAtt                                             >  7:48653/1‑144 (MQ=255)
                                                                                           cGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGc                                       >  6:280237/1‑149 (MQ=255)
                                                                                             ttCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGa                                                       <  4:51381/131‑1 (MQ=255)
                                                                                             ttCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGa                                                       >  3:51381/1‑131 (MQ=255)
                                                                                               cGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAAACGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTg                                    <  5:129102/148‑1 (MQ=255)
                                                                                                        cGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCg                             >  1:86279/1‑146 (MQ=255)
                                                                                                        cGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCg                             <  2:86279/146‑1 (MQ=255)
                                                                                                         gTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGt                         >  1:332232/1‑149 (MQ=255)
                                                                                                          tGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTg                        >  5:328614/1‑149 (MQ=255)
                                                                                                           gCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGc                       >  5:185121/1‑149 (MQ=255)
                                                                                                             aaaGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGctg                     <  5:280237/149‑1 (MQ=255)
                                                                                                                          cGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGctgttcctgttcctgt        <  6:185121/149‑1 (MQ=255)
                                                                                                                               tGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGcc   <  4:70955/149‑1 (MQ=255)
                                                                                                                                 cgcAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCt  <  2:332232/148‑1 (MQ=255)
                                                                                                                                 cgcAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCt  <  8:67892/148‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                                   
TACGTCTGGCGTCGACACCTGCCAACCGACTTCGACGCCTTGCTGGCACGCCTGCGGCCGCTGGCCATCTGCGGGTTGGTAGCGCCACGGGCGTTCGCCGATGCCGTGCAAAGCATGTGGCACGGGGTGCGCAACGGCGAAATCGTCCACCAGTTGCTGGTGTTGCGTCCGGCGCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCT  >  NC_002947/3010256‑3010532

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: