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.C93.1% 81.7 / 0.5 29intergenic (‑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 (17/10);  total (19/10)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

CTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATT  >  NC_002947/1932500‑1932780
                                                                                                                                            |                                                                                                                                              
cTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGGCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAg                                                                                                                                        <  3:68291/148‑1 (MQ=255)
 ttCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAgg                                                                                                                                       >  4:71876/1‑148 (MQ=255)
     aaGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCCACCCCCCTCCGACAAAGGCAATTGGTCCCCCGGAGGAGCGGGCGTGCCCGCGACCGGCTg                                                                                                                                   >  1:156515/1‑148 (MQ=255)
     aaGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTATGAGCGGGCTTGCCCGCGAAGGGCt                                                                                                                                    >  8:46040/1‑147 (MQ=255)
       gagGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAACCCGCTCCGACAAAGGGTATGGGTCCCCGGTAGGAGCGGGCGTGCCCGCGAAGGGCTGCa                                                                                                                                 >  5:135310/1‑148 (MQ=38)
                   aTTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAACCCCCACCGACAAAGGATATGCGGGCCCGGTAGGAGCGGGCGTGCCCGCGCAGGGCTGCAACGCccccccc                                                                                                                      >  2:76326/1‑146 (MQ=255)
                    ttCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCaaa                                                                                                                    >  4:22938/1‑148 (MQ=255)
                       gACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTATGCGGCCCTTCGCGGGCAAGCCCGCAGCGACAAACGATATGCCTCCCCTGTAGGAGCGGGCATGCCCCCGCACGGCAGCAAGGACGCCCCAAAATg                                                                                                                 >  7:162994/1‑148 (MQ=255)
                         ccgCGGCTTTGCAGGTCTCACTAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCAACAAAGTATATGGGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTc                                                                                                               <  4:89118/146‑1 (MQ=37)
                         cAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGCATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTc                                                                                                               <  8:6014/148‑1 (MQ=255)
                                            aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCGGTAGGAGCGGGCTTGGCCGCGCAGGGCTGGCAGGCAGCCCCAAATTGTCGACTACGCGGGTCAGCGAt                                                                                            >  2:54786/1‑148 (MQ=38)
                                                    aGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTGCGCGGGCAAGCCCGCTCCTACAAAGTATATGCGTCCCCTTTTGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAgg                                                                                    <  2:37095/148‑1 (MQ=38)
                                                    aGCAGAACACGCGGGGGCCGCTTTGCGGCCCTGCGCGGGCAAGCCCGCTCCTACAACGTATAGGCGTCCCCTGTAGGAGCGGGCTGGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAgg                                                                                    <  3:19899/148‑1 (MQ=37)
                                                       aGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCGGTAGGCGCGGGCTTGCCCGCGCAGGGCTGCAAGGCAGCCCCAAATTGTCGCCTACGCGGCTCAGCGATCCGGCCGGGTg                                                                                 >  7:121726/1‑148 (MQ=37)
                                                        gTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCCCGCCGACAAAGGATATGCGGCCCCGGTAGGAGCGGGCTTGCCCGCGCAGGGCAGGAAGGCAGCCCCCAATTTTAGCCTACGCCGATAAGCGATCCCGCAGGGGGa                                                                                >  8:157923/1‑148 (MQ=21)
                                                         tAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCCAACCCCCAACAACAAACGCTATGCGGGCCCAGTAGGAGCGGGCTTGGCCCCCAACGGCTGCAAGGCCGCCCCCAATTGTCCACAACGCGGAAAACCGAACCGGCaccgggaa                                                                               >  6:509/1‑141 (MQ=255)
                                                          aaccgcgGGGGCCGCTTGGCGGCCCTTCGCGGGCACGCCCGCGCCTACAACGGATAGGCGGCCCCTGTCGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATg                                                                              <  5:78559/145‑1 (MQ=37)
                                                           aaCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCCGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGt                                                                             >  8:23840/1‑148 (MQ=255)
                                                            aCGCGGGGGCGGCTTTGCGGCCTTTCGGGGGAAAGCCGGATCATACAAAGTATATGGGTCCCTTTTAGTTGCGGGGTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGtt                                                                            <  2:46234/148‑1 (MQ=255)
                                                                                              gCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTcc                                          <  3:71876/148‑1 (MQ=255)
                                                                                                 cGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGc                                       >  5:7365/1‑148 (MQ=255)
                                                                                                  gCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCt                                      >  3:67671/1‑148 (MQ=255)
                                                                                                       tACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGata                                 >  2:199780/1‑148 (MQ=255)
                                                                                                                tatGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatca                        >  3:51435/1‑148 (MQ=255)
                                                                                                                     gTCCCCTGTTGGTGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGtt                   <  1:9218/148‑1 (MQ=255)
                                                                                                                       ccccTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGc                 >  5:150017/1‑148 (MQ=255)
                                                                                                                          cTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGa              >  4:30793/1‑148 (MQ=255)
                                                                                                                               gttgggggCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATTCACATCATCGTTGCCGATGtt          <  6:145880/142‑1 (MQ=255)
                                                                                                                                      gCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATAtt  >  8:52421/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                              
CTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATT  >  NC_002947/1932500‑1932780

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: