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.C96.8% 101.1 / ‑3.0 31intergenic (‑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 (18/12);  total (19/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATAT  >  NC_002947/1932494‑1932779
                                                                                                                                                  |                                                                                                                                             
cggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGccc                                                                                                                                            <  5:145350/148‑2 (MQ=255)
      cTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGcccgcgcag                                                                                                                                      >  7:107742/1‑141 (MQ=255)
       ttCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAgg                                                                                                                                      >  7:47213/1‑148 (MQ=255)
          gAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCCCTCCTACAAAGGATATGCGTCCCCTGTATGAGCGGGCTTTCCCGCGAAGGGCt                                                                                                                                   >  1:55307/1‑148 (MQ=38)
                       aCATTCCACAGCGGAATTGCAGGTCTCACTAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCACCCCCGCGCCTAAAACGTATATGCGTCCCCTGTCGGAGGGGGCTTGCCCGCGAAGGGCTGCAAGGCa                                                                                                                           <  2:123333/143‑1 (MQ=255)
                              aCAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGt                                                                                                               >  3:163439/1‑148 (MQ=255)
                               cAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGt                                                                                                               >  6:178022/1‑147 (MQ=255)
                                cgCGGCATTGAAGGTCTAACTAAAATAGCATTAAACGCGGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAATGTATATGCGGCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGt                                                                                                               <  6:176960/145‑1 (MQ=25)
                                             gTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCCAGCCCGCGCCTACAAAGGATATGCGTCCCCGGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCCGCCCCAAATTGTCGCCTACGCGGATCa                                                                                                >  2:41554/1‑148 (MQ=38)
                                              tctcCCGAAAATAGCTGTAAACGCTGGGGCCGGTTTGCGGCCCGTCGCGGGCACGCCCGCTCCTACAACGGATTTCCGTCCCGTGTCGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAg                                                                                               <  2:109586/148‑1 (MQ=37)
                                               ctcACTAAAATAGCTGTAAAGGTTGGGGCGGATTTGGGGCCCTTGGGGGGCATCCCCGCTCATACAAAGTATATGCTTCCCCTGTAGTAGGGGGGTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAg                                                                                               <  2:55307/147‑1 (MQ=255)
                                                tcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTAAGCGGATCAGc                                                                                              >  4:28713/1‑147 (MQ=255)
                                                          aGCATTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGGCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAgg                                                                                   <  6:178776/148‑1 (MQ=255)
                                                          aGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAgg                                                                                   <  3:58293/148‑1 (MQ=255)
                                                           gCAGTAAACGCTGGGGCCGCTTGGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATAGGCGGCCCCTGTCGGAGCGGGCTGGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCGCTGATCCgcgtagt                                                                                  <  7:17070/148‑7 (MQ=25)
                                                              gTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCGCCTACAAAGGCTATGCGTCCCCGGTAGGAGCGGGCGTGCCCGCGCAGGGCTGCCAGGCCGCCCCAAATTGTCGCCTCCGCGGCTCCGCGCTCCGGCCGGGGGa                                                                               >  5:193205/1‑148 (MQ=25)
                                                                  aCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAACGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGt                                                                            >  8:47580/1‑147 (MQ=255)
                                                                            cGCTTTGCGGCCCTTCGCGGGCAAGCCCCCTCCAAAAAACGATATTCGTACCCTGTATGAGCCGGCTTGCCCGCGAAGGGCAGCAAGGCAGCCCCAAAATGTCGACAACGCGGATCAGCGATCCCGCCCGGTGATGTTCAGCTCCAg                                                                  >  2:52550/1‑147 (MQ=255)
                                                                                        cTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGCGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCCAATTGTCGACTACGCGGCTGAGCGATCCGGCAGGGGGATGTTCAGCTCCAGAATCGAGCAGCTg                                                     >  6:98374/1‑148 (MQ=38)
                                                                                                    gCCCGCTCATACAAAGGATATGCTTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTcc                                         <  3:188526/148‑1 (MQ=255)
                                                                                                     cccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCa                                        <  5:178022/148‑1 (MQ=255)
                                                                                                      ccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAg                                       >  5:38649/1‑148 (MQ=255)
                                                                                                      ccGCTCCTACAAAGGATATGCGCCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAg                                       >  7:19352/1‑148 (MQ=255)
                                                                                                       cGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAAt                                                                                                                  >  6:163091/1‑72 (MQ=255)
                                                                                                        gCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAAtt                                                                                                                 <  5:163091/72‑1 (MQ=255)
                                                                                                               cAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGGTTTCCAGCTCGatat                               >  7:56853/1‑147 (MQ=255)
                                                                                                                aaaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGc                             >  1:32869/1‑148 (MQ=255)
                                                                                                                aaaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCAGGCAGGGTGATGTTCAGCTCCAGAATCGACCAGCTGCCCTGGTTTTCCAGCTCGATATGc                             >  8:152944/1‑148 (MQ=255)
                                                                                                                         gCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCg                    >  8:9051/1‑148 (MQ=255)
                                                                                                                          cGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGt                   >  3:99211/1‑148 (MQ=255)
                                                                                                                                   tAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGt          >  1:46955/1‑148 (MQ=255)
                                                                                                                                      gAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTc        <  4:99211/147‑1 (MQ=255)
                                                                                                                                           ggCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACatat  <  7:110870/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATAT  >  NC_002947/1932494‑1932779

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: