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.2% 167.5 / 0.3 52intergenic (‑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 (23/27);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

GCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACAT  >  NC_002947/1932493‑1932777
                                                                                                                                                     |                                                                                                                                            
gcggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGccc                                                                                                                                           <  8:23697/149‑2 (MQ=255)
        ttCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAggg                                                                                                                                   >  5:108497/1‑149 (MQ=255)
            aaGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTg                                                                                                                                >  3:159303/1‑148 (MQ=255)
                       gACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTg                                                                                                                                >  1:78905/1‑137 (MQ=255)
                       gACATTCCACAGCGGCATTGCAGGTCTCACTAAAATAGCATTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTg                                                                                                                                <  2:78905/137‑1 (MQ=37)
                              gaccggggCATTGCTGGTCTCACGAAATTAGCAGTAAAGGGTGGGGCCGCTTTGCGGCCCTTCGGGGGCAAGCCGGCTCCTACAATGGATTTGCGTCCCCTGTAGGTGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTg                                                                                                              <  2:248048/142‑1 (MQ=37)
                              gACTGCGGCATTGCAGGTCTCACTAAAATACCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAgg                                                                                                                                    <  1:326371/126‑1 (MQ=37)
                                  gCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGAc                                                                                                         <  6:118396/149‑1 (MQ=255)
                                     gCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTa                                                                                                       >  4:119967/1‑148 (MQ=255)
                                     gCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATCCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTa                                                                                                       <  3:119967/148‑1 (MQ=255)
                                      cATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTAcg                                                                                                     >  5:170914/1‑149 (MQ=255)
                                           cAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGAt                                                                                                >  1:97286/1‑149 (MQ=255)
                                               tctcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGc                                                                                            >  6:84256/1‑149 (MQ=255)
                                                 tcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGa                                                                                          >  5:327443/1‑149 (MQ=255)
                                                   actaaaaTCGCAGTAACCGCGGGGGCCGTTTTGCGGCCGTGGGCGGGCAAGCCCGCTCATAAAACGGATATGCGCCCCTTTTTGGAGGGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATc                                                                                        <  4:159303/146‑1 (MQ=255)
                                                   aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTACTACAAAGGATATGAGTGCCCGGTAGGAGCGGGCTTGGCCGCGAAGGGCTTGAAGGCAGCCCCAAATTTTCGACAACCCGGATACGCGAta                                                                                        >  8:175089/1‑148 (MQ=37)
                                                            gctttAAACGGTGTGGCCGCTTTGCGGCCGTTGGCTGGCAAGCCGGCGACTACAAAGGTTATTCGTCCCCTGTAGTAGCGGGATTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGt                                                                               <  5:132192/145‑1 (MQ=255)
                                                            gaaGTATACGCTGTGGGCGCTTTGCGGCCCTTCTCGGGCAAGCCCGCTCATACAAATGATATGCGTCCCCTGTAGGAGCTGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGt                                                                               <  5:138182/147‑1 (MQ=38)
                                                              aGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGa                                                                             >  3:332182/1‑149 (MQ=255)
                                                               gTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGCTCAGCGATCCGGCAGGGTGAt                                                                            >  7:330025/1‑149 (MQ=255)
                                                                 aaaCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATg                                                                           <  6:327443/148‑1 (MQ=255)
                                                                   aCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTc                                                                        <  5:84256/149‑1 (MQ=255)
                                                                   aCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGCATATGCGGCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTc                                                                        <  3:170661/149‑1 (MQ=255)
                                                                       tGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCGACAAAGGATATGCGTCCCCGGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCCAGGCCGCCCCAAATTGTCGACTACGCGGATCAGCGAtg                                                                                        >  3:173367/1‑128 (MQ=37)
                                                                       tGGGGCCGCTTGGCGGCCCTGCGCGGCCACCCCCGCTCCTACAACGCATTTGCGCCCCCGGTAGCCGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATc                                                                                        <  4:173367/129‑1 (MQ=25)
                                                                           gCCGCTTTGCGGCCCTTCGCGGGCAAGCCCCCAACTACAAAGGATATTCGTCCCCGGTAGGAGCGGGCGTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTTTGGACAAAGCCGATCAGCGATCCGGGAGGGTGATGTTCAGCtg                                                                   >  2:164353/1‑145 (MQ=25)
                                                                             cGCTTTGCGGCCCTTCGCGGGCCAACCCCCAACAAAAAAAGGATATGCGTCCCATTTAGAACCGGGATTGCCCCCCAAGGCCTGAAAGCCACCCCCAAAACGACGAAAACCCGGAACAGCGAACCGCAAGGGTGATTTTAAgaacgga                                                                >  2:347593/1‑142 (MQ=255)
                                                                             cGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCATACAATGTTTATGCGTCCCCTTTCGGTGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATCGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGtt                                                                         <  1:347593/138‑1 (MQ=37)
                                                                                    cGGCCCTTGGCGGGAAAGCCGGATCATACAATGTATTTGCGTCCCGTGTAGTTGGGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagca                                                       <  2:97286/149‑1 (MQ=255)
                                                                                     ggCCCTTCGCGGGCAAGCCCGCTCCTACAGAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGa                                                                             >  1:247105/1‑127 (MQ=255)
                                                                                     ggCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGa                                                                             <  2:247105/126‑1 (MQ=255)
                                                                                                 gCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGttt                                          >  6:115619/1‑149 (MQ=255)
                                                                                                 gCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGttt                                          >  5:3239/1‑149 (MQ=255)
                                                                                                    acgCCCGCTCCTCCAATGGATTTGCGTCCCAGGTCGTAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc                                        <  8:330025/147‑1 (MQ=38)
                                                                                                    aaGCCCGCTCCTACAAAGTATATGCGCCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc                                        <  3:157939/149‑1 (MQ=255)
                                                                                                              cTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGata                              <  4:332182/149‑1 (MQ=255)
                                                                                                              cTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGa                                <  6:108497/147‑1 (MQ=255)
                                                                                                                aCAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATg                            >  6:62545/1‑149 (MQ=255)
                                                                                                                  aaaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATc                                                            >  4:237938/1‑115 (MQ=255)
                                                                                                                  aaaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATc                                                            <  3:237938/115‑1 (MQ=255)
                                                                                                                  aaaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGca                          <  7:175089/149‑1 (MQ=255)
                                                                                                                    aaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagc                                                        <  4:264838/118‑1 (MQ=255)
                                                                                                                    aaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagc                                                        >  3:264838/1‑118 (MQ=255)
                                                                                                                      ggATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcat                       >  5:398976/1‑149 (MQ=255)
                                                                                                                           tGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGag                                                         <  4:278998/110‑1 (MQ=255)
                                                                                                                           tGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGag                                                         >  3:278998/1‑110 (MQ=255)
                                                                                                                           ggCGGCCCCTGTAGGAGCGGGCTGGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACatcatc                   <  7:156347/147‑1 (MQ=255)
                                                                                                                            gCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGt                 >  8:262082/1‑149 (MQ=255)
                                                                                                                               tCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGc              <  5:62545/149‑1 (MQ=255)
                                                                                                                                        ggAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTc      <  7:244234/148‑1 (MQ=255)
                                                                                                                                        ggAGCGGGCTTGCCAGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCCGCAGGGTGATGTTCATCTCCAGAATCGAGCAGCTCCCCTCGTTTACCAGCTCGATATGCACATCATCGTTGCCGATGTTca     >  6:410548/1‑149 (MQ=255)
                                                                                                                                         gAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTca     >  4:411125/1‑148 (MQ=255)
                                                                                                                                           gCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACat  <  6:3239/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACAT  >  NC_002947/1932493‑1932777

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: