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.C97.7% 145.7 / ‑3.1 44intergenic (‑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 (25/18);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

CGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATT  >  NC_002947/1932494‑1932780
                                                                                                                                                   |                                                                                                                                                 
cggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATGGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCACTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGcccg                                                                                                                                               >  1:357547/1‑147 (MQ=255)
       ttCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCCACCCCCCTACTAACAAGGATAATCGTCCCCCATATGAGACGGCGTGCCCGCGAAggg                                                                                                                                         >  6:330270/1‑149 (MQ=255)
                gATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGAATTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGAAAGCCGGCTCCTAAAAAGGATTTGCGTCCCGTGTAGGTGGGGGGTTGCCCGCGAAGGGCTGCAAGGc                                                                                                                                <  8:38456/149‑1 (MQ=255)
                 aTCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAgg                                                                                                                                 >  1:300907/1‑147 (MQ=255)
                             gACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGt                                                                                                                   >  6:172934/1‑149 (MQ=255)
                                aGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGa                                                                                                                <  7:62752/149‑1 (MQ=255)
                                   ggCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCCCTCCAACAAAGGATATGCGTACCCTGTAGGAGCCGGCTTTCCCGCGAAGGGCTGCAAGGCAGCCCCAAAATGTCGAcaa                                                                                                             >  2:328111/1‑147 (MQ=37)
                                         gCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCgg                                                                                                        >  7:220938/1‑148 (MQ=255)
                                             ggctcACGAAAATCGCAGTAACCGCTGGGGCCGCTTTGCGGCCCTGCGCGGGCAAGCCCGCTCCTACAACGCATATGCGTCCCCTGTCGGAGCGGGGTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTc                                                                                                                  <  4:130986/132‑1 (MQ=25)
                                             gTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTc                                                                                                                  >  3:130986/1‑134 (MQ=255)
                                                tcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGa                                                                                                >  1:268302/1‑149 (MQ=255)
                                                  ccGAAAAAAGCCTTAACCGCTGGGGCCGCTTTGCGGCCCTGCGCGGGCAAGCCCGCTCCTACAAAGCATATGCGTCCCCTGTCGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGa                                                                                                <  5:330270/146‑1 (MQ=38)
                                                  aCGAAAATCGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCgg                                                                                                        <  2:286501/139‑1 (MQ=255)
                                                  aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCgg                                                                                                        >  1:286501/1‑140 (MQ=255)
                                                  aCGAAAATACCAGTAAACGCGGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATACCCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTc                                                                                                                  <  7:100004/129‑1 (MQ=37)
                                                                          gCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTc                                                                      >  6:358732/1‑146 (MQ=255)
                                                                            cGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGaa                                                                 <  2:93658/149‑1 (MQ=255)
                                                                             gCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCATGTAGGAGCGGGATTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCAGGCAGGGTGATGTTCAGCTCCAGAAt                                                                >  6:182329/1‑149 (MQ=38)
                                                                              cTTTGCGGCCTTTCGCGGGCAAGCCCGCTCATACAAAGTATATGCGTCCCTTTTTGGTGGGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATc                                                               <  2:300907/149‑1 (MQ=255)
                                                                                    agCCCTTCGCGGGCAAGCCCGCTCCTACAGGGGACGCATATCCTTTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagcag                                                         >  6:182475/2‑149 (MQ=18)
                                                                                               gggCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGTGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGa                                                             >  2:331251/1‑135 (MQ=255)
                                                                                               gggCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGTGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGa                                                             <  1:331251/135‑1 (MQ=255)
                                                                                                     gCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc                                           >  6:301536/1‑147 (MQ=255)
                                                                                                      cccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAAtt                                                                                                                     <  6:380983/75‑1 (MQ=255)
                                                                                                      cccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAAtt                                                                                                                     >  5:380983/1‑75 (MQ=255)
                                                                                                      cccACTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAg                                        <  2:357547/149‑1 (MQ=255)
                                                                                                       ccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCg                                                              <  1:349484/126‑1 (MQ=255)
                                                                                                        cGCTCCTACAATGGATATGGGTCCCCTGTAGTTGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCt                                      <  1:328111/149‑1 (MQ=255)
                                                                                                           tCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGTAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGa                                   >  1:278721/1‑149 (MQ=255)
                                                                                                             cTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGata                                 >  7:306967/1‑149 (MQ=255)
                                                                                                                 aaaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                                >  4:391636/1‑146 (MQ=255)
                                                                                                                 aaaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                                <  3:391636/146‑1 (MQ=255)
                                                                                                                 aaaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGca                             >  1:6673/1‑149 (MQ=255)
                                                                                                                        atGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACatcatc                      >  3:317895/1‑149 (MQ=255)
                                                                                                                         tGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACatcatc                      <  5:358732/148‑1 (MQ=255)
                                                                                                                           cGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGtt                   >  3:221983/1‑149 (MQ=255)
                                                                                                                             tcctttGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGc                 <  5:182475/144‑1 (MQ=255)
                                                                                                                             tCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGc                 <  5:172934/149‑1 (MQ=255)
                                                                                                                               cccTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGcc                >  8:254284/1‑148 (MQ=255)
                                                                                                                                      ggAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTca        >  5:343417/1‑149 (MQ=255)
                                                                                                                                      ggAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTca        >  4:296149/1‑149 (MQ=255)
                                                                                                                                          cGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACata    >  3:262505/1‑149 (MQ=255)
                                                                                                                                          cGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACata    <  2:268302/149‑1 (MQ=255)
                                                                                                                                           gggCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGAATGTAGGAACGAACAGAGAGGAGCAGCAGCCCTGGTTTTACAGCTCGATATTCACATCATAGTTGCAGATGTTcaca      >  1:321568/1‑149 (MQ=255)
                                                                                                                                            ggCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATAtt  >  2:115446/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATT  >  NC_002947/1932494‑1932780

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