Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,601,331 C→T 12.9% intergenic (‑148/‑39) bglX ← / → PP_1404 beta‑D‑glucoside glucohydrolase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,601,3310CT12.9% 73.8 / 3.3 31intergenic (‑148/‑39)bglX/PP_1404beta‑D‑glucoside glucohydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (13/14);  new base T (2/2);  total (15/16)
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.16e-01

GTTTTGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGCCATTGCCCGTGCCGAACGGCATACCGATGCCGAACTGGTGACGGTGCTTGCCCGTCGTGCCGACGA  >  NC_002947/1601191‑1601479
                                                                                                                                            |                                                                                                                                                    
gTTTTGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATaa                                                                                                                                               >  3:201563/1‑148 (MQ=255)
 ttttGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATaac                                                                                                                                              >  6:75052/1‑148 (MQ=255)
   ttGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATaacaa                                                                                                                                            <  4:71583/148‑1 (MQ=255)
      gggCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACta                                                                                                                                         <  3:51335/148‑1 (MQ=255)
              gCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGc                                                                                                                                >  8:213338/1‑149 (MQ=255)
                cGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAAc                                                                                                                                  <  6:258765/145‑1 (MQ=255)
                cGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAAc                                                                                                                                  >  5:258765/1‑145 (MQ=255)
                  tATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCaa                                                                                                                            <  4:154615/149‑1 (MQ=255)
                                       gACAGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGAccccc                                                                                                       >  7:202397/1‑149 (MQ=255)
                                         ccGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCtt                                                                                                     <  4:201563/149‑1 (MQ=255)
                                           gAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTcg                                                                                                   <  7:213338/149‑1 (MQ=255)
                                             aGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTcgac                                                                                                 <  5:75052/149‑1 (MQ=255)
                                                 ttCAGGCGGCTATTGTTGTACGATTTTCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCTGGGATaaca                                                                                                                                             >  8:119437/1‑101 (MQ=255)
                                                 ttCAGGCGGCTATTGTTGTACGATTTTCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCTGGGATaaca                                                                                                                                             <  7:123662/101‑1 (MQ=255)
                                                 ttCAGGCGGCTATTGTTGTACGATTTTCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCTGGGATaaca                                                                                                                                             <  7:119437/101‑1 (MQ=255)
                                                 ttCAGGCGGCTATTGTTGTACGATTTTCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCTGGGATaaca                                                                                                                                             >  8:123662/1‑101 (MQ=255)
                                                         gcTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAg                                                                                     >  3:24935/1‑149 (MQ=255)
                                                                   tACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTg                                                                           >  4:269301/1‑149 (MQ=255)
                                                                                    tcttcCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGCCATTGCCCGt                                                          >  3:153738/1‑149 (MQ=255)
                                                                                         cAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGCCATTGCCCGTGCCg                                                      >  4:72529/1‑148 (MQ=255)
                                                                                         cAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGCCATTGCCCGTGCCg                                                      >  8:148869/1‑148 (MQ=255)
                                                                                                       gCCTGTATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGCCATTGCCCGTGCCGAACGGCATACCGATg                                       >  5:277681/1‑149 (MQ=255)
                                                                                                        ccTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGTGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGCCATTGCCCGTGCCGAACGGCATACCGATGc                                      <  8:202397/149‑1 (MQ=255)
                                                                                                         cTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGCCATTGc                                                              <  2:108726/124‑1 (MQ=255)
                                                                                                              tatGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGCCATTGCCCGTGCCGAACGGCATACCGATGCCGAACt                                <  3:269301/149‑1 (MQ=255)
                                                                                                                    gCCTGCGGATGAAGCCCCGAATGCCGTGATAACAACTATAAGAAAGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGCCTTTGCCCGTGCCGAACGGCATACCGATGCCGAACTGGTGa                           <  3:206296/148‑1 (MQ=255)
                                                                                                                         cGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGCCATTGCCCGTGCCGAACGGCATACCGATGCCGAACTGGTGACGGt                       <  7:148869/147‑1 (MQ=255)
                                                                                                                              gAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGCCATTGCCCGTGCCGAACGGCATACCGATGCCGAACTGGTGACGGTGCTTGcc                >  2:306326/1‑149 (MQ=255)
                                                                                                                                    ccGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGCCATTGCCCGTGCCGAACGGCATACCGATGCCGAACTGGTGACGGTg                      >  2:236677/1‑137 (MQ=255)
                                                                                                                                    ccGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGCCATTGCCCGTGCCGAACGGCATACCGATGCCGAACTGGTGACGGTg                      <  1:236677/137‑1 (MQ=255)
                                                                                                                                       aaTGCCGGGATAAGAACTATAAGAACGCTCAAGAGACCGACGGCATGACCACCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGCCATTGCCCGTGCCGAACGGCATACGGATGCCGAACTGGTGACGGTGCTTGCCCGTCGTGcc       <  8:62395/149‑1 (MQ=255)
                                                                                                                                            cGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGGCATTGCCCGTGCCGAACGGCATACCGATGCCGAACTGGTGACGGTGCTTGCCCGTCGTGCcgacga  >  1:52655/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
GTTTTGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTGCCGAAGCCATTGCCCGTGCCGAACGGCATACCGATGCCGAACTGGTGACGGTGCTTGCCCGTCGTGCCGACGA  >  NC_002947/1601191‑1601479

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