Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 728,963 C→A 25.0% T150N (ACC→AAC)  rluD → 23S rRNA pseudouridine(1911/1915/1917) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947728,9630CA25.0% 34.1 / 6.5 16T150N (ACC→AAC) rluD23S rRNA pseudouridine(1911/1915/1917) synthase
Reads supporting (aligned to +/- strand):  ref base C (7/5);  new base A (2/2);  total (9/7)
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.00e+00

CAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCTGATGGTGGTGGCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCA  >  NC_002947/728844‑729094
                                                                                                                       |                                                                                                                                   
cAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGGGCCGCGCGCCGGTATCGTGCACCGCCGGGACAAGGACACCACCGGCCTGATGGTGGTGGCCAAGACCCTGCAGGCGCAGAcc                                                                                                                     <  1:247147/136‑1 (MQ=255)
                        gctgcACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCTGATGGTGGTGGCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGcc                                                                                             <  2:78418/136‑1 (MQ=255)
                        gcggcACCATGGGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCTGATGGTGGTGGCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGcc                                                                                             <  2:285323/133‑1 (MQ=255)
                                      cAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCTGATGGTGGTGGCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCg                                                                               >  1:275045/1‑136 (MQ=255)
                                                cAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCTGATGGTGGTGGCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAg                                                                     >  2:69161/1‑136 (MQ=255)
                                                       ccGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCTGATGGTGGTGGCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCg                                                              <  2:234143/136‑1 (MQ=255)
                                                         gcgcgcCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCTGATGGTGGTGGCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTg                                                            <  1:69161/136‑1 (MQ=255)
                                                                          aCCGCCTGGACAAGGACACCACCGGCCTGATGGTGGTGGCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGc                                           >  1:45220/1‑136 (MQ=255)
                                                                                 ggACAAGGACACCACCGGCCTGATGGTGGTGGCCAAGAACCTGCAGGCGCAGACCAAGCt                                                                                                                <  1:246135/60‑1 (MQ=255)
                                                                                 ggACAAGGACACCACCGGCCTGATGGTGGTGGCCAAGAACCTGCAGGCGCAGACCAAGCt                                                                                                                <  1:247586/60‑1 (MQ=255)
                                                                                 ggACAAGGACACCACCGGCCTGATGGTGGTGGCCAAGAACCTGCAGGCGCAGACCAAGCt                                                                                                                >  2:246135/1‑60 (MQ=255)
                                                                                 ggACAAGGACACCACCGGCCTGATGGTGGTGGCCAAGAACCTGCAGGCGCAGACCAAGCt                                                                                                                >  2:247586/1‑60 (MQ=255)
                                                                                                                gCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCAcggcg     >  1:125894/1‑136 (MQ=255)
                                                                                                                gCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCAcggcg     >  2:29147/1‑136 (MQ=255)
                                                                                                                gCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCAcggcg     >  2:8904/1‑136 (MQ=255)
                                                                                                                   aaGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCa  >  1:275343/1‑136 (MQ=255)
                                                                                                                       |                                                                                                                                   
CAGCGGCACCCTGCTCAACGCGCTGCTGCACCATGTGCCAGACATCGTCAACGTGCCGCGCGCCGGTATCGTGCACCGCCTGGACAAGGACACCACCGGCCTGATGGTGGTGGCCAAGACCCTGCAGGCGCAGACCAAGCTGGTCGAGCAGATGCAGGCCCGCAAGGTCAGCCGCATCTACGAGTGCATCGTGATTGGCGTGGTGACTGCCGGTGGCAAGATCGACGCGCCCATCGGCCGCCACGGCGGCA  >  NC_002947/728844‑729094

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