Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,601,266 G→T 14.3% intergenic (‑83/‑104) 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,2660GT14.3% 72.9 / 4.8 28intergenic (‑83/‑104)bglX/PP_1404beta‑D‑glucoside glucohydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (12/12);  new base T (2/2);  total (14/14)
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.43e-01

GGGCGGCGAGGGCCGCCTGACTGGCAAGGCCCATGGCCAGGCCCAGCAAAGACAGTTTCATCATCAATTCCGTTTTGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTG  >  NC_002947/1601120‑1601406
                                                                                                                                                  |                                                                                                                                            
gggCGGCGAGGGCCGCCTGACTGGCAAGGCCCATGGCCAGGCCCAGCAAAGACAGTTTCATCATCAATTCCGTTTTGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCa                                                                                                                                            <  8:30696/149‑1 (MQ=255)
           gccgccTGACTGGCAAGGCCCATGGCCAGGCCCAGCAAAGACAGTTTCATCATCAATTCCGTTTTGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAAttcttc                                                                                                                                 <  4:156441/149‑1 (MQ=255)
                   aCTGGCAAGGCCCATGGCCAGGCCCAGCAAAGACAGTTTCATCATCAATTCCGTTTTGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAgcggcg                                                                                                                         <  1:55050/149‑1 (MQ=255)
                     tGGCAAGGCCCATGGCCAGGCCCAGCAAAGACAGTTTCATCATCAATTCCGTTTTGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAgg                                                                                                                                          <  1:295546/130‑1 (MQ=255)
                     tGGCAAGGCCCATGGCCAGGCCCAGCAAAGACAGTTTCATCATCAATTCCGTTTTGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAgg                                                                                                                                          >  2:295546/1‑130 (MQ=255)
                        cAAGGCCCATGGCCAGGCCCAGCAAAGACAGTTTCATCATCAATTCCGTTTTGGGGCCTCTCCCGTATCAACGGTGATACGAGCGCGACCGAACTTTACAGGCGGCTATTGGTGTACGATTTGCAGGCAATTGTACCAGCGGCGGATTa                                                                                                                    >  4:206296/1‑149 (MQ=255)
                                      aGGCCCAGCAAAGACAGTTTCATCATCAATTCCGTTTTGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTc                                                                                                      <  8:292597/149‑1 (MQ=255)
                                         cccAGCAAAGACAGTTTCATCATCAATTCCGTTTTGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTTGcc                                                                                                   >  7:62395/1‑149 (MQ=255)
                                                              atcaATTCCGTTTTGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGc                                                                              >  4:51335/1‑149 (MQ=255)
                                                              atcaATTCCGTTTTGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGc                                                                              >  3:71583/1‑149 (MQ=255)
                                                                       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                                                         >  5:258765/1‑145 (MQ=255)
                                                                                       cGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAAc                                                         <  6:258765/145‑1 (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                                                                    <  7:123662/101‑1 (MQ=255)
                                                                                                                        ttCAGGCGGCTATTGTTGTACGATTTTCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCTGGGATaaca                                                                    >  8:119437/1‑101 (MQ=255)
                                                                                                                        ttCAGGCGGCTATTGTTGTACGATTTTCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCTGGGATaaca                                                                    >  8:123662/1‑101 (MQ=255)
                                                                                                                        ttCAGGCGGCTATTGTTGTACGATTTTCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCTGGGATaaca                                                                    <  7:119437/101‑1 (MQ=255)
                                                                                                                                gcTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAg            >  3:24935/1‑149 (MQ=255)
                                                                                                                                          tACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTg  >  4:269301/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
GGGCGGCGAGGGCCGCCTGACTGGCAAGGCCCATGGCCAGGCCCAGCAAAGACAGTTTCATCATCAATTCCGTTTTGGGGCCTCTGCCGTATCCACGGTGATACGTGCGCGACCGAAGTTTTCAGGCGGCTATTGTTGTACGATTTGCAGGCAATTCTTCCAGCGGCGGATTATGCCTGAATATGTCGCCTGCGGATGAAGCCCCGAATGCCGGGATAACAACTATAAGAACGCTCAAGAGACCGACGGCATGACCCCCTTCGACGAATACCATCAGCGGCAGATTG  >  NC_002947/1601120‑1601406

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: