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.0% 106.4 / ‑3.0 33intergenic (‑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 (17/15);  total (18/15)
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.83e-01

ATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTT  >  NC_002947/1932490‑1932781
                                                                                                                                                          |                                                                                                                                               
aTGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑g                                                                                                                                                  <  1:6708/149‑1 (MQ=255)
       cgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑Gcccgcg                                                                                                                                           >  1:274952/1‑144 (MQ=255)
             cGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACTAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAggg                                                                                                                                       <  6:131970/147‑1 (MQ=255)
               aaGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGTAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTg                                                                                                                                    >  1:273396/1‑148 (MQ=255)
                                cGACAGCGGCTTTGCAGGTCTCACTAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCGTTCGGGGGCATGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAAtt                                                                                                                   <  8:219461/148‑1 (MQ=38)
                                   cAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTc                                                                                                                <  2:273396/148‑1 (MQ=255)
                                         cATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGCTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTAc                                                                                                          <  1:102494/149‑1 (MQ=255)
                                              cAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGa                                                                                                     >  3:117657/1‑148 (MQ=255)
                                               aGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATc                                                                                                   <  2:13112/149‑1 (MQ=255)
                                                 gTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAg                                                                                                 <  2:274952/149‑1 (MQ=255)
                                                      aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTAcgc                                                                                                        >  6:120212/1‑137 (MQ=255)
                                                      aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTAcgc                                                                                                        <  5:120212/137‑1 (MQ=255)
                                                      aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATc                                                                                            >  5:122033/1‑149 (MQ=255)
                                                                   tAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATg                                                                               >  5:118699/1‑149 (MQ=255)
                                                                                   tttGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGTCGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGGGATGTTCAGCTCCAGAAt                                                                 >  1:118146/1‑149 (MQ=38)
                                                                                               cgcgGGCAAGCCCGCTCCTACACAGTATATGCGTCCCCTGTAGGCGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGccc                                                    <  2:118146/148‑1 (MQ=39)
                                                                                                 cgGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCCGCTCCAGAATCGAGCAGCTGCCCTgg                                                 >  7:367474/1‑149 (MQ=255)
                                                                                                  tgtcacgCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGTTGCGGGATT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGt                                                <  4:233037/143‑1 (MQ=255)
                                                                                                    gCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATc                                                                                            >  1:391408/1‑103 (MQ=255)
                                                                                                      aaGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc                                            <  4:117657/149‑1 (MQ=255)
                                                                                                      aaGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc                                            <  6:122033/149‑1 (MQ=255)
                                                                                                         cccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAg                                         >  7:272264/1‑149 (MQ=255)
                                                                                                                 tACAAAGGATATGCGTCCCCTGTAGGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGata                                  >  1:4374/1‑149 (MQ=255)
                                                                                                                 tACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                                 >  8:181649/1‑149 (MQ=255)
                                                                                                                  aCAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATc                                                                >  2:83330/1‑117 (MQ=255)
                                                                                                                  aCAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATc                                                                <  1:83330/117‑1 (MQ=255)
                                                                                                                     aaGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGcac                             <  7:181649/149‑1 (MQ=255)
                                                                                                                      aGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGcaca                            >  3:298979/1‑149 (MQ=255)
                                                                                                                             gCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGt                     >  6:81374/1‑149 (MQ=255)
                                                                                                                               gTCGCCTGTAGGCGGGGTCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTg                   <  3:273775/149‑1 (MQ=255)
                                                                                                                                    cTGTAGGAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt              >  3:118680/1‑149 (MQ=255)
                                                                                                                                           ggAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGc                  >  2:228078/1‑140 (MQ=255)
                                                                                                                                           ggAGCGGGCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGc                  <  1:228078/140‑1 (MQ=255)
                                                                                                                                                 ggCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAAACGAGCAGCTGCCCTGGTTTTCCAGCTCGATATTAAAATCATAGTTGACGATGTTGACATAtt   >  6:333134/1‑149 (MQ=255)
                                                                                                                                                  gCTT‑‑‑GCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATAttt  >  6:310900/1‑149 (MQ=255)
                                                                                                                                                          |                                                                                                                                               
ATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑‑‑GCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTT  >  NC_002947/1932490‑1932781

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: