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.C95.7% 64.6 / ‑2.6 23intergenic (‑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 . (1/0);  new base C (10/12);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

TCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCA  >  NC_002947/1932499‑1932774
                                                                                                                                             |                                                                                                                                        
tCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGGCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGTGGGCAAGCCCGCTCCTACAGGGGACGCATATCCTTTGTAGGAGCGGGCTTGCCCGCGaa                                                                                                                                   <  3:78191/148‑1 (MQ=18)
 ctgctaaGAGGATCTGCACATTCGACAGCGGCATTTCAGGTCTAACGAAAATAGCAGTAAACGCTGGGGACGCTTTGCGGCCCTGCGCGGGCATGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGTAGCGGGCTTGCCCGCGAAg                                                                                                                                  <  6:173055/143‑1 (MQ=255)
             tCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCa                                                                                                                      >  7:132882/1‑148 (MQ=255)
                 gACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTACTAAAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAACGGCAGCAAGGccccc                                                                                                                   >  2:66879/1‑143 (MQ=38)
                     ttCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTg                                                                                                                                           <  1:6336/119‑1 (MQ=255)
                        caccgCGGAATTGCCGGTCTCACTAAAATAGCCGTAAACGCGGGGGCCGCTTGGCGGCCCTGCGCGGGCACGCCCGCTCCTACAACGGATACGCGGCCCCTGTCGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAAtt                                                                                                            <  8:3910/143‑1 (MQ=25)
                          ccgCGGCATTGCAGGGCTCACTAAAATAGCAGTAACCGGGGGGGCGGCTTTGCGGCCCGGCGCGGGAAAGCCCGCTCCAACAACGCATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTc                                                                                                         <  8:111802/146‑1 (MQ=25)
                          cAGCGGCATTGCAGGGCGCACGAAAATAGCAGTAACCGCGGGGGCCGCTTGGCGGCCCTTCGCGGGCACGCCCGCTCCTACAACGGATAGGCGGCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTc                                                                                                         <  1:14/148‑1 (MQ=37)
                          cAGCGGCAGTGCAGGGCCCACGAAAATCCCAGAAAACGCGGGGGCCGCTTTGCGGCCCGTCGCGGGCAACCCCGCTCCAACAAAGGATATGCGGCCCCTGAAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTc                                                                                                         <  3:166530/148‑1 (MQ=37)
                          cAGCGGAATTCCAGGTCTCACTAAAATAGCCGTAAACGCGGGGGCCGCTTTGCGGCCCTTCGCGGGCACGCCCGCTCATACAACGGATAGGGGCCCCGTGTCGGCGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTc                                                                                                         <  6:161411/148‑1 (MQ=25)
                                        gTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCa                                                                                           >  5:82634/1‑148 (MQ=255)
                                           tcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCGACAAAGGATATGCGTCCCCTGTAGGCGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTTTCGACTACGCGGATCAGCg                                                                                        >  8:80389/1‑148 (MQ=38)
                                             aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAACCCCCGCCTAAAAAGGATATGCGGCCCCAGTAGGAGCGGGCGTTCCCCCCAAAGGCTGCAAGGCAGCCCCCAAATTTGGACAAACCCGATCAGCgct                                                                                      >  2:68683/1‑146 (MQ=255)
                                                                          tttGCGGCCGTGCGGGGGCACGCCCGCTCCTACAAAGTATAGGCGGCCCCTGTCGGCGGGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATc                                                         <  6:113772/148‑1 (MQ=255)
                                                                                                 ccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAg                                  >  1:79904/1‑148 (MQ=255)
                                                                                                 ccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGATCCAGAATAGAGCAGCTGCCCTGGTTTTCAAg                                  >  3:57402/1‑148 (MQ=39)
                                                                                                         aCAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                          >  1:175626/1‑148 (MQ=255)
                                                                                                             aGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGcac                      >  4:152216/1‑148 (MQ=255)
                                                                                                              ggATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGcaca                     >  4:61849/1‑148 (MQ=255)
                                                                                                              ggATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGcaca                     <  2:175626/148‑1 (MQ=255)
                                                                                                                                 gAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTca  <  3:19326/148‑1 (MQ=255)
                                                                                                                                 gAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTc   <  2:28602/147‑1 (MQ=255)
                                                                                                                                 gAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTc   <  6:82634/147‑1 (MQ=255)
                                                                                                                                 gAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATAAGCGATCCGGCAGGCAGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTACAGCTCGATATGCACATCATCGTTGCCGATGtt  >  2:114421/1‑146 (MQ=39)
                                                                                                                                             |                                                                                                                                        
TCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCA  >  NC_002947/1932499‑1932774

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: