Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,438,757 T→C 8.9% Q200R (CAG→CGG)  rluE ← 23S rRNA pseudouridine(2457) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,438,7570TC 8.9% 122.3 / 2.5 45Q200R (CAG→CGG) rluE23S rRNA pseudouridine(2457) synthase
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base C (2/2);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

TTCTTCGTGCTGGGCCTGGGGGTGCTGGCCTTCATCATCAAAAGCGTGGTGGTGGGGTTTCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCTTCGCGAATCACCAGTTCCAGCCAGCTGGTCGGCACGCTCTTGC  >  NC_002947/5438613‑5438904
                                                                                                                                                |                                                                                                                                                   
ttcttcGTGCTGGGCCTGGGGGTGCTGGCCTTCATCATCAAAAGCGTGGTGGTGGGGTTTCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTc                                                                                                                                                 >  7:154149/1‑149 (MQ=255)
         cTGGGCCTGGGGGTGCTTGCCTTCATCATCAAAAGCGTGGTGGCGGGGTTGCTGGAAGGCGTATAAAAGATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCTCCAGCAACCCTGCTCCAGCCCATc                                                                                                                                        <  3:198400/149‑1 (MQ=255)
          tGGGCCTGGGGGTGCTGGCCTTCATCATCAAAAGCGTGGTGGTGGGGTTTCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATc                                                                                                                                        >  8:59846/1‑148 (MQ=255)
                     gTGCTGGCCTTCATCATCAAAAGCGTGGTGGTGGGGTTTCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTc                                                                                                                            <  8:154149/149‑1 (MQ=255)
                         tGGCCTTCATCATCAAAAGCGTGGTGGTGGGGTTTCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCAcc                                                                                                                         >  5:94846/1‑148 (MQ=255)
                           gCCTTCATCATCAAAAGCGTGGTGGTGGGGTTTCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAAt                                                                                                                      <  2:58833/149‑1 (MQ=255)
                             cTTCATCATCAAAAGCGTGGTGGTGGGGTTTCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATgc                                                                                                                    >  5:59288/1‑149 (MQ=255)
                                catcatcaAAAGCGTGGTGGTGGGGTTTCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATgcgc                                                                                                                  <  2:224752/148‑1 (MQ=255)
                                     tcaAAAGCGTGGTGGTGGGGTTTCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATgcgcacgcgc                                                                                                            >  7:14986/1‑149 (MQ=255)
                                       aaaaGCGTGGTGGTGGGGTTTCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATgcgcacgcgcac                                                                                                          <  3:35541/149‑1 (MQ=255)
                                        aaaGCGTGGTGGTGGGGTTTCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATgcgcacgcgcac                                                                                                          >  6:25509/1‑148 (MQ=255)
                                                        gTTTCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCa                                                                                                                               >  4:53008/1‑111 (MQ=255)
                                                         tttCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGc                                                                                          >  2:54038/1‑147 (MQ=255)
                                                         tttCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCaa                                                                                        <  7:59846/149‑1 (MQ=255)
                                                         tttCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAACGCGCACGCGCACCAGGCGCAAGGTCGGCaa                                                                                        >  6:203293/1‑149 (MQ=255)
                                                             tGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGcc                                                                                     >  6:115796/1‑148 (MQ=255)
                                                                          aaaCATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCa                                                                       <  6:59288/149‑1 (MQ=255)
                                                                             cATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAg                                                                                                >  2:118940/1‑121 (MQ=255)
                                                                             cATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAg                                                                                                <  1:118940/121‑1 (MQ=255)
                                                                               tttCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGGCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATccgc                                                                  <  5:25509/149‑1 (MQ=255)
                                                                                     ggCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCAcc                                                            <  5:115796/149‑1 (MQ=255)
                                                                                       ccATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTg                                                          >  5:90321/1‑149 (MQ=255)
                                                                                          ttCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCg                                                       >  6:273522/1‑149 (MQ=255)
                                                                                                ggCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACg                                                 <  1:54038/149‑1 (MQ=255)
                                                                                                 gCCCTGACCTGATCTAAAGCCTGGCAGGGACCTCCAGCCAGCAACCCCGCTCAAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATGCGCCGCACCTGGCGGTTaaga                                                >  2:57165/1‑146 (MQ=255)
                                                                                                 gCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCCGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGc                                                >  2:60937/1‑149 (MQ=255)
                                                                                                  cccTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGc                                                                                          >  8:36220/1‑106 (MQ=255)
                                                                                                  cccTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGc                                                                                          <  7:36220/106‑1 (MQ=255)
                                                                                                    ctgacctgaACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAACGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCt                                             <  5:203293/149‑1 (MQ=255)
                                                                                                               tAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGc                                                        >  7:43517/1‑127 (MQ=255)
                                                                                                               tAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGc                                                        <  8:43517/127‑1 (MQ=255)
                                                                                                                 aaGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGTCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCTTCGCGAATCAc                                  <  8:14986/147‑1 (MQ=255)
                                                                                                                 aaGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCTTCGCGAATCAc                                  <  6:94846/147‑1 (MQ=255)
                                                                                                                 aaGCCTGGCAGGGACCTCCCGCCAGCAACCCCGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCTTCGCGAATCACCa                                <  1:60937/149‑1 (MQ=255)
                                                                                                                 aaGCCTGGCAGGGACCTCCCGCCAGCAACCCCGCTCCAGCACATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCTTCGCGAATCACCa                                <  1:57165/149‑1 (MQ=255)
                                                                                                                  aGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCTTCGCGAATCACCAg                               <  6:90321/149‑1 (MQ=255)
                                                                                                                     cTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCTTCGCGAATCACCAGTTc                            >  5:240120/1‑149 (MQ=255)
                                                                                                                          aGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATccgc                                                                  >  7:52224/1‑106 (MQ=255)
                                                                                                                          aGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCTTCGCGAATCACCAGTTccagc                        >  8:247811/1‑148 (MQ=255)
                                                                                                                               ccTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCTTCGCGACTCACCAGTTCCAGCCAGCTg                  >  6:142286/1‑149 (MQ=255)
                                                                                                                                cTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCTTCGCGAATCACCAGTTCCAGCCAGCTgg                 >  8:193199/1‑149 (MQ=255)
                                                                                                                                           aaCCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCTTCGCGAATCACCAGTTccagcca                      <  3:22535/133‑1 (MQ=255)
                                                                                                                                           aaCCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCTTCGCGAATCACCAGTTccagcca                      >  4:22535/1‑133 (MQ=255)
                                                                                                                                             cccTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCTTCGCGAATCACCAGTTCCAGCCAGCTGGTCGGCACGCTCtt    <  7:247811/149‑1 (MQ=255)
                                                                                                                                               cTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCTTCGCGAATCACCAGTTCCAGCCAGCTGGTCGGCACGCTATTGc  >  6:14389/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TTCTTCGTGCTGGGCCTGGGGGTGCTGGCCTTCATCATCAAAAGCGTGGTGGTGGGGTTTCTGGAAGGCGTATAAAACATTTCAGGGCCATTCGCGGGCCCTGACCTGAACTAAAGCCTGGCAGGGACCTCCCGCCAGCAACCCTGCTCCAGCCCATCCAGCGCCCAGTCACCAATGCGCACGCGCACCAGGCGCAAGGTCGGCAAGCCCACTGCTGCCGTCATCCGCCGCACCTGGCGGTTACGCCCTTCGCGAATCACCAGTTCCAGCCAGCTGGTCGGCACGCTCTTGC  >  NC_002947/5438613‑5438904

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: ‑