Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,144,974 C→T A12V (GCT→GTT)  rluC → 23S rRNA pseudouridine(955/2504/2580) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,144,9740CT95.5% 61.6 / NA 22A12V (GCT→GTT) rluC23S rRNA pseudouridine(955/2504/2580) synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (13/8);  minor base G (0/1);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GGCCTCATTATTCCACTGCAAGCCTTGTTATAGCAAGATGACTTTTACCATTTATCACCCGCTTACTCACAGTTTTTTCACTTCTTGCTGGTGATTGGTTTAATAACCATCAATAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGCTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAATATAAACTCGAAGCGGGTGATGAGGTGCGTATTCCACCGGTTCGCGTTGCTGAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGT  >  NC_000913/1144729‑1145220
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
ggCCTCATTATTCAACTGCAAGCCTTGTTATAGCAAGATGACTTTTACCATTTATCACCCGCTTACTCACAGTTTTTTCACTTCTTGCTGGTGATTGGTTTAATAACCATCAATAAAGTGATCACGAAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATc                                                                                                                                                                                                                                                    >  1:321660/1‑250 (MQ=255)
    tCATTATTCCACTGCAAGCCTTGTTATAGCAAGATGACTTTTACCATTTATCACCCGCTTACTCACAGTTTTTTCACTTCTTGCTGGTGATTGGTTTAATAACCATCAATAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGc                                                                                                                                                                                                                                               >  1:55297/1‑251 (MQ=255)
         attCCACTGCAAGCCTTGTTATAGCAAGATGACTTTTACCATTTATCACCCGCTTACTCACAGTTTTTTCACTTCTTGCTGGTGATTGGTTTAATAACCATCAATAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGAc                                                                                                                                                                                                                                           <  1:404167/250‑1 (MQ=255)
             cACTGCAAGCCTTGTTATAGCAAGATGACTTTTACCATTTATCACCCGCTTACTCACAGTTTTTTCACTTCTTGCTGGTGATTGGTTTAATAACCATCAATAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGc                                                                                                                                                                                                                                      >  1:360995/1‑251 (MQ=255)
                 gCAAGCCTTTGTATAGCAAGATGACTTTTACCAGTTATCCCCCGCTTACTCACAGTTTTTTCACTTCTTGCTGGTGATTGGTTTAATAACCATCAATAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATATAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCggg                                                                                                                                                                                                                                   <  2:428663/250‑1 (MQ=255)
                           ttATAGCAAGATGACTTTTACCATTTATCACCCGCTTACTCACAGTTTTTTCACTTCTTGCTGGTGATTGGTTTAATAACCATCAATAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTTCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCg                                                                                                                                                                                                                        >  1:398566/1‑251 (MQ=255)
                                       gACTTTTACCATTTATCACCCGCTTACTCACAGTTTTTTCACTTCTTGCTGGTGATTGGTTTAATAACCATCAATAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGc                                                                                                                                                                                                            >  1:192349/1‑251 (MQ=255)
                                           tttaccatttaccaCCCGCTTACTCACAGTTTTTTCACTTCTTGCTGGTGATTGGTTTAATAACCATCAATAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTa                                                                                                                                                                                                         <  2:166475/250‑1 (MQ=255)
                                           tttACCATTTATCACCCGCTTACTCACAGTTTTTTCACTTCTTGCTGGTGATTGGTTTAATAACCATCAATAAAGTGATCACGCAGTGAAGTAATAAGGCGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCACCCATGAAGACAGAGAATCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTAc                                                                                                                                                                                                        >  1:316612/1‑251 (MQ=255)
                                                                            ttCACTTCTTGCTGGTGATTGGTTTAATAACCATCAATAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGAttt                                                                                                                                                                        >  1:180588/1‑250 (MQ=255)
                                                                                        tGGTGATTGGTTTAATAACCATCAATAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGAGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAACAATTGTTGTTATCACCGCTGCCGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCg                                                                                                                                                           >  2:132964/1‑251 (MQ=255)
                                                                                                                    aGTGAACACGCAGTGAAGTAATAAGGGGCAAAAGTAACTATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACTATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGATGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTAACCAATTGAAAGGCGTACCAAAAAGTATGATTTACAGTATTTTGCGTAACGGAGAAGAGCGGGTGAACacaac                                                                                                                                >  2:109662/1‑246 (MQ=255)
                                                                                                                        aTCACGCAGTGAATTAATAAGGTGAAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTTGCTAATTAGAATCGCCAACCATGTAAACAGAGACTCTATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTa                                                                                                                           <  1:284534/251‑1 (MQ=255)
                                                                                                                                  gAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTg                                                                                                                  >  1:11511/1‑250 (MQ=255)
                                                                                                                                     gTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAat                                                                                                               >  1:33297/1‑250 (MQ=255)
                                                                                                                                                           ataAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGATAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGACTGAATATAAACTCGAAGCggggtatg                                                                                          >  1:371997/1‑244 (MQ=255)
                                                                                                                                                                       aaatgctgggCGCGAATCTGGTTGGTTAATTAGACTCGCCAACCATTAAAACAGAGACTCCATCCTTAAAAATTGGTGTTATCACCGCTTACGAAGCGGGGCAACGTATCTATAACTTTTGGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAATATAACCTCGCAGCTGGTGTTGAGGTGAGTATTCCa                                                                            <  2:381931/243‑1 (MQ=255)
                                                                                                                                                                                   cgAATCCTGTGGGATAATTAGAATCGCCAACCATGAAAACAGAGACTACATCCGTAAAAATTGTGGTTACCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAATATAAACTCGAAGCGGGTGATGAGGTGCGTATTCCACCGGTTCGCGtt                                                                <  2:192349/251‑1 (MQ=255)
                                                                                                                                                                                           gTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAATATAAACGCGAAGCGGGTGATGAGGTGCGTATTCCACCGGTTCGCGTTGCTGAGCg                                                        <  1:250287/251‑1 (MQ=255)
                                                                                                                                                                                              ggCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAATATAAACTCGAAGCGGGTGATGAGGTGCGTATTCCACCGGTTCGCGTTGCTGAGCggg                                                      <  1:216200/250‑1 (MQ=255)
                                                                                                                                                                                               gctaattacaatacccaccaaTGAAAACAGAGATCCCATACTTAAAATTTGTTTGTATTAACGATTACGAATCAGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCCGACCAAAAAGTATGATTTACCGTATTTTGCTTAAAGGCGCAGTTCGTGTGAATAAAAAACGTAGTAATCCTGAATAGTAACTCGAAGAGGGTTAAGAGGTGCGTATTCCACCTGTCCTCGTTGCTGAGCgttaa                                                    <  2:303898/231‑5 (MQ=255)
                                                                                                                                                                                                         aTCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGTTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAATATAAACTCGAAGCGGGTGATGAGGTGCGTATTCCACCGGTTCGCGTTGCTGAGCGGGAAGAAGAGGCgg                                          >  1:325154/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 gttgttATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAATATAAACTCGAAGCGGGTGATGAGGTGCGTATTCCACCGGTTCGCGTTGCTGAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGt  <  2:290624/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
GGCCTCATTATTCCACTGCAAGCCTTGTTATAGCAAGATGACTTTTACCATTTATCACCCGCTTACTCACAGTTTTTTCACTTCTTGCTGGTGATTGGTTTAATAACCATCAATAAAGTGATCACGCAGTGAAGTAATAAGGTGCAAAAGTAAATATAAGCATAGAAAAATGAGTGGCGCGAATCCTGTTGGCTAATTAGAATCGCCAACCATGAAAACAGAGACTCCATCCGTAAAAATTGTTGCTATCACCGCTGACGAAGCGGGGCAACGTATCGATAACTTTTTGCGTACCCAATTGAAAGGCGTACCAAAAAGTATGATTTACCGTATTTTGCGTAAAGGCGAAGTGCGGGTGAACAAAAAACGTATTAAGCCTGAATATAAACTCGAAGCGGGTGATGAGGTGCGTATTCCACCGGTTCGCGTTGCTGAGCGGGAAGAAGAGGCGGTTTCGCCACATCTGCAAAAGGTGGCGGCGCTGGCGGACGT  >  NC_000913/1144729‑1145220

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: