Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,932,574 T→A 7.7% intergenic (‑41/+109) 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,5740TA 7.7% 144.0 / 3.2 52intergenic (‑41/+109)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 T (27/21);  new base A (2/2);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

CCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATC  >  NC_002947/1932428‑1932719
                                                                                                                                                   |                                                                                                                                                     
ccTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTg                                                                                                                                                     >  4:271935/1‑149 (MQ=255)
ccTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCACACGAAAATAGCAGTAAACGCTGGGGCAGCTTTg                                                                                                                                                     >  4:285476/1‑149 (MQ=255)
 cTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGc                                                                                                                                                    <  5:108877/149‑1 (MQ=255)
    gttgtCCTCGGCACGGCCAGGTACCGACAGCAGCAAGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCgg                                                                                                                                                  <  7:39514/148‑1 (MQ=255)
          cTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTcg                                                                                                                                           >  1:8222/1‑149 (MQ=255)
             ggCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTcgcg                                                                                                                                         <  7:40802/148‑1 (MQ=255)
              gCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCggg                                                                                                                                       <  2:6061/149‑1 (MQ=255)
                      aGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCg                                                                                                                               >  5:242838/1‑149 (MQ=255)
                        gTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCt                                                                                                                             <  7:71902/149‑1 (MQ=255)
                          aCCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTc                                                                                                                            >  2:98724/1‑148 (MQ=255)
                           ccGACAGCAGCAGGGTGGGTTTGTTGCTCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGcgaa                                                                                                                           <  4:217493/148‑4 (MQ=39)
                             gACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTa                                                                                                                        >  8:271558/1‑148 (MQ=255)
                              aCAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCAAAGCGGccccagcg                                                                                                                                          >  5:73137/1‑125 (MQ=255)
                              aCAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCAAAGCGGccccagcg                                                                                                                                          <  6:73137/130‑6 (MQ=255)
                                agcagcagGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACaaa                                                                                                                    <  6:242838/149‑1 (MQ=255)
                                agcagcagGGTAGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACaaa                                                                                                                    <  6:291161/149‑1 (MQ=255)
                                   agcagGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGa                                                                                                                 >  4:104301/1‑149 (MQ=255)
                                    gcagGGTCGGTTTTTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGat                                                                                                                <  3:90923/149‑1 (MQ=255)
                                      agGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGatat                                                                                                              <  3:104301/149‑1 (MQ=255)
                                       gggTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATg                                                                                                             <  2:8222/149‑1 (MQ=255)
                                                 gttgATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTg                                                                                                    <  1:130862/148‑1 (MQ=255)
                                                                    gcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGcccgcg                                                                             >  6:169425/1‑145 (MQ=255)
                                                                         ttCGAAGAGGATCTGGACATTCAACAGCGGCATTGCAGGTCTCAAGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAACCCCCACCAAAAAAAGAAATGCGTGCCCGGTAGGAACGGGCTTGACCGCCAAggg                                                                         >  7:194149/1‑149 (MQ=255)
                                                                           cGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTcgcg                                                                                                                                         >  8:174358/1‑86 (MQ=255)
                                                                           cGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTcgcg                                                                                                                                         <  7:174358/86‑1 (MQ=255)
                                                                                  gATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAgg                                                                 >  2:27508/1‑148 (MQ=255)
                                                                                    tCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAg                                                              >  7:31086/1‑149 (MQ=255)
                                                                                      tGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTATGAGCGGGCTTGCCCGCGAACGGCTGCAAGGCAGcc                                                            >  4:10415/1‑149 (MQ=255)
                                                                                         aCATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTGCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCCTTGCCCGCGAAGGGCTGCAAGGCAGcccc                                                          <  1:27508/149‑1 (MQ=255)
                                                                                           aTTCGACAGCGGCATTGCAGGTCTCACGAAAAGTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGat                                                                                                                >  1:36583/1‑95 (MQ=255)
                                                                                           aTTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGat                                                                                                                <  2:36583/94‑1 (MQ=255)
                                                                                             tCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAACCCCCCAACAACAACAGCTATTCGTCCCCAGTAGGAGGCGGGGTGGCAGCCCAGGGCGGACAGGCGGCCCCAAAtt                                                     >  8:248482/1‑149 (MQ=255)
                                                                                                    cGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCggg                                                                                         <  8:247461/107‑1 (MQ=255)
                                                                                                    cGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCggg                                                                                         >  7:247461/1‑107 (MQ=255)
                                                                                                            cAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCAACAAAGGAAATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGAt                                      >  1:94763/1‑149 (MQ=255)
                                                                                                                tctcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGc                                  >  1:63850/1‑149 (MQ=255)
                                                                                                                   cACGAAAATACCAGTAAACGCTGGGGCCGCTTTGCGGCCGTTCGCGGGCAAGCCCGCTCCTAAAATGGATATGCGTCCCCTGTAGTAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGa                                <  7:271558/148‑1 (MQ=38)
                                                                                                                    aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATc                              >  5:295220/1‑149 (MQ=255)
                                                                                                                            tAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAg                                   >  8:59538/1‑137 (MQ=255)
                                                                                                                            tAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAg                                   <  7:59538/137‑1 (MQ=255)
                                                                                                                                 ttacccGCTGGGCCCGCTTGGCGGCCCTGCGCGGGAACCCCCGCTCATACAACGTATAGGCGTCCCCTGTTGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGAt                  <  5:169425/144‑1 (MQ=25)
                                                                                                                                 gTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCCCTCCTACAAAGGATATTCGTCCCCGGTAGGAGCCGGCGTGCCCGCGAACGGCTGCAAGGCCGCCCCAAATTTTCGACAACCCGGATAAGCCAACCCGCAGGGTGAt                  >  4:192784/1‑149 (MQ=25)
                                                                                                                                 gTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCCCTCCGACAAAGGATTTGTGTACCCAGTAGGAGAGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCCAATTTTCGACGAAGCGGAAAAGCCAAACGGCAGGGTGAt                  >  2:269327/1‑149 (MQ=25)
                                                                                                                                 atAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCGCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCCGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGa                   >  1:67184/2‑149 (MQ=38)
                                                                                                                                    aaCGCTGGGGCCGCTTTGCGGCCATTCGCGGGCAAGCCCGCTACAAAAAAGCATATTCCTCCACAGTAGGATACGGCTTGCCCGCGAAGTGCTTCAAGGCACGCCCAAAATGTAGAATACGCCGATCAGCGAAACGGCACGGTGAttt                >  8:25735/1‑146 (MQ=255)
                                                                                                                                     aCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGCCTACGCGGATCAg                                   >  1:203903/1‑128 (MQ=255)
                                                                                                                                         aagggCCGCAAAGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTc              >  7:70397/3‑145 (MQ=38)
                                                                                                                                         aagggCCGCAAAGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTc              <  8:70397/143‑1 (MQ=38)
                                                                                                                                           gggCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCt          <  2:67184/147‑1 (MQ=255)
                                                                                                                                             gCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTc         >  8:177107/1‑146 (MQ=255)
                                                                                                                                              ccGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAg      >  4:222766/1‑148 (MQ=255)
                                                                                                                                                 cTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATc  <  8:31086/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATC  >  NC_002947/1932428‑1932719

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: