Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,932,640 +C intergenic (‑107/+43) rluA ← / ← minE bifunctional tRNA pseudouridine(32) synthase TrmQ/ribosomal large subunit pseudouridine synthase RluA/cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,932,6381.C91.7% 69.8 / 0.9 24intergenic (‑105/+45)rluA/minEbifunctional tRNA pseudouridine(32) synthase TrmQ/ribosomal large subunit pseudouridine synthase RluA/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (11/11);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

TCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTC  >  NC_002947/1932499‑1932773
                                                                                                                                              |                                                                                                                                       
ggttCGAAGAGTAGCTGGACATTCGACTGCGGCATTGAAGGTCTCACTAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCGTTGGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGaa                                                                                                                                  <  2:99733/146‑1 (MQ=38)
         cgcaTCTGGACATTCGACCGCGCCATTGCAGGCCTCACTAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCCCGGGCAAGCCCGCTCCTACAAAGCATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCaa                                                                                                                         <  6:23910/145‑1 (MQ=38)
              cTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTGCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCCAGGCAg                                                                                                                    >  8:89683/1‑148 (MQ=255)
               tGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGCTATGCGTCCCCGGTAGGAGCGGGCTTGCCCGCGCAGGGCTGCCAGGCAGc                                                                                                                   >  3:6651/1‑148 (MQ=38)
                ggACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGcc                                                                                                                  >  4:9999/1‑148 (MQ=255)
                    aTTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCaa                                                                                                              >  3:36837/1‑148 (MQ=255)
                             cGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAACGAAATGCCGCCCCTGTTGGGGCGGGCTTGGCCCCGCACGGGTGCCAGGGAGCACCAAAATggagg                                                                                                      >  4:262/1‑143 (MQ=25)
                                 aTTGCAGGTCTCACTAAAATAGCAGTAACCGCGGGGGCCGCTTTGCGGCCCTGCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGGGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTAcg                                                                                                 <  8:4100/148‑1 (MQ=37)
                                           tcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAACCCCGCGCCTAAAAAGGATATTCGGCCCCAGTTGGAGACGGCGTTCCCGCGCAGGGCAGCAAGGCAACCCCAAATTGTCGCAGACGCCGATCAGc                                                                                        >  6:218/1‑147 (MQ=255)
                                             aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCGGGAGGAGCGGGCTTGCCCGCGCAGGGCTGCAAGGCCGCCCCAAATTGTCGACTACGCGGCTCAGCGAt                                                                                     >  6:120069/1‑148 (MQ=37)
                                              cGAAAATAGCAGTAAACGCGGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGCATAGCCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATc                                                                                    <  4:107219/148‑1 (MQ=37)
                                                   aTAGCAGTAACCGCTGGGGCCGCTTGGCGGCCCTTCGCGGGCAACCCCGCTCCTACAAATGATAGGCGGCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCa                                                                               <  8:120979/148‑1 (MQ=38)
                                                    tAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCGCTGATCCgcgta                                                                              <  4:36837/148‑5 (MQ=37)
                                                                          tttGGGGCCGTTCGCGGGCAAGCCCGCGCCTACAACGTATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATc                                                        <  4:6651/148‑1 (MQ=38)
                                                                               ggCCGTGCGCGGGCAAGCCCGCTCCTACAACGGATATGCGGCCCCTGTAGGAGGGGGCTGGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagca                                                   <  5:38566/148‑1 (MQ=255)
                                                                                 cccTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagcagc                                                 >  6:31723/1‑148 (MQ=255)
                                                                                             aaGCCCGCTCCTACAAAGGATATGGGTCCCCTGTAGGAGCGGGCTGGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGtttt                                     <  3:127769/148‑1 (MQ=255)
                                                                                             aaGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGGTCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGtttt                                     >  6:25311/1‑148 (MQ=255)
                                                                                                           aaaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGc                       >  5:28407/1‑148 (MQ=255)
                                                                                                             cggatcggCGGCCCCGGTAGGCGGGGGGTGGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGcac                     <  3:3586/141‑1 (MQ=255)
                                                                                                                atatGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatc                  <  4:36870/148‑1 (MQ=255)
                                                                                                                      gTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGtt            >  2:78392/1‑148 (MQ=255)
                                                                                                                            tGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt      >  8:91512/1‑148 (MQ=255)
                                                                                                                                 gAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTc  >  2:10720/1‑147 (MQ=255)
                                                                                                                                              |                                                                                                                                       
TCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTC  >  NC_002947/1932499‑1932773

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