Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,112,642 T→A 9.1% intergenic (‑231/+35) rluB ← / ← PP_4497 23S rRNA pseudouridylate synthase B/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,112,6420TA 9.1% 110.0 / 3.4 44intergenic (‑231/+35)rluB/PP_449723S rRNA pseudouridylate synthase B/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (17/23);  new base A (2/2);  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 = 2.29e-01

GGGGAACGCACTGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATCCAGCTCCACCAACAGGC  >  NC_002947/5112495‑5112789
                                                                                                                                                   |                                                                                                                                                   
ggggAACGCACTGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTg                                                                                                                                                    >  5:343033/1‑149 (MQ=255)
       gCACTGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCCACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCcgcg                                                                                                                                              >  6:154811/1‑148 (MQ=255)
       gCACTGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCAAGCCcgctc                                                                                                                                             <  7:98839/149‑3 (MQ=255)
             gAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAagag                                                                                                                                         <  8:201607/147‑1 (MQ=255)
             gAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                       >  2:5243/1‑149 (MQ=255)
                 gTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGc                                                                                                                                   <  3:225191/149‑1 (MQ=255)
                    gagGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGca                                                                                                                                  >  3:215512/1‑147 (MQ=255)
                     agGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCACCCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGcaca                                                                                                                                <  5:154811/148‑1 (MQ=255)
                      gggTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTAGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCTCGCGAAGAGGCCAGCACAgg                                                                                                                              >  6:62577/1‑149 (MQ=255)
                       ggTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCGTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGgc                                                                                                                             <  6:290637/149‑1 (MQ=255)
                                      ggCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAgg                                                                                                              >  8:6848/1‑149 (MQ=255)
                                         aTGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCt                                                                                                           <  7:66669/149‑1 (MQ=255)
                                             ttttCTGCCGAGTGGGCTTTGTGCTGGGCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTTATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGc                                                                                                       <  5:313852/149‑1 (MQ=255)
                                                  tGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCGACATGGTCGCGCAACCTTCGGGGCCTGTGATACCCCAGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGCCCAGGCAGCAGCTCAAGCCGGGCTGAGCCTTg                                                                                                   >  7:252426/1‑148 (MQ=38)
                                                   gCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGAt                                                                                                 >  1:254072/1‑149 (MQ=255)
                                                    cccAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGAt                                                                                                 <  4:215512/148‑1 (MQ=255)
                                                        gTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATctcct                                                                                            >  7:303628/1‑149 (MQ=255)
                                                            gCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATctcctcctc                                                                                        <  2:267244/149‑1 (MQ=255)
                                                                    cTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCTCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGa                                                                                 <  5:62577/148‑1 (MQ=255)
                                                                       gCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTg                                                                                                   >  4:117012/1‑127 (MQ=255)
                                                                       gCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTg                                                                                                   <  3:117012/127‑1 (MQ=255)
                                                                       gCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGaggct                                                                             >  3:137338/1‑149 (MQ=255)
                                                                       gCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGaggct                                                                             >  4:46923/1‑149 (MQ=255)
                                                                        ccTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAATAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGaggctc                                                                            <  8:303628/149‑1 (MQ=255)
                                                                              cgcgGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGaggctcaggc                                                                        <  7:6848/147‑1 (MQ=255)
                                                                                   gCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTGGATCTCCTCCTCCACCGGAGGCTCAGGCTCCtctc                                                                 >  7:191034/1‑149 (MQ=255)
                                                                                         gcGCTCCTCGATGGTCGCGCAACCTGCGGGGCCTGTGATACCCCTGGGGGCGCGGGCTTGCCCGCGAAGAGGACAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCa                                                           <  8:252426/148‑1 (MQ=255)
                                                                                            cTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGat                                                        <  1:147601/149‑1 (MQ=255)
                                                                                                  cATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGAtcatca                                                   <  2:254072/148‑1 (MQ=255)
                                                                                                      ttCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGACTCCTCTCCCCGCAGATCATCAAAAt                                               <  1:5243/148‑1 (MQ=255)
                                                                                                                cTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCaaa                                    >  3:150792/1‑149 (MQ=255)
                                                                                                                 ttCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAAc                                   >  4:55007/1‑149 (MQ=255)
                                                                                                                     gggCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCgg                                            <  7:7167/136‑1 (MQ=255)
                                                                                                                     gggCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCgg                                            >  8:7167/1‑136 (MQ=255)
                                                                                                                            tgATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCa                          <  4:137338/147‑1 (MQ=255)
                                                                                                                            tgATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATg                        <  8:191034/149‑1 (MQ=255)
                                                                                                                              aTACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATgg                       <  3:46923/148‑1 (MQ=255)
                                                                                                                                  cccTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATcc                  <  4:150792/149‑1 (MQ=255)
                                                                                                                                  cccTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATc                   >  1:48122/1‑148 (MQ=255)
                                                                                                                                       tgGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATCCAGc               <  8:217095/147‑1 (MQ=255)
                                                                                                                                         ttcgCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGtt                                          <  7:264704/115‑1 (MQ=255)
                                                                                                                                         ttcgCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGtt                                          >  8:264704/4‑118 (MQ=255)
                                                                                                                                         ttcgCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATCCAGCTcca           >  1:43805/4‑149 (MQ=255)
                                                                                                                                                 cTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATCCAGCTCCACCAACAgg   <  2:48122/149‑1 (MQ=255)
                                                                                                                                                  ttGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATCCAGCTCCACCAACAGGc  <  6:343033/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGGGAACGCACTGGAGGGTAGAGGGTTTTGTGGGGCTGGGCATGATTTTCTGCCCAGTGGGCTTTGTGCTGGCCTCTTCGCGGGCACGCCCGCTCCTACATGTTCGCGCAACCTTCGGGGCCTGTGATACCCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCAGCACAGGCAGCAGATCAATCAGGGCTGAGCCTTGATCTCCTCCTCCACCGGAGGCTCAGGCTCCTCTCCCCGCAGATCATCAAAATCGGTTTTCAAACCGTCTTCCATGGCATCCAGCTCCACCAACAGGC  >  NC_002947/5112495‑5112789

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