Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,932,638:1 +C 100% intergenic (‑105/+45) 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.C90.0% 118.6 / 8.0 40intergenic (‑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 . (3/1);  new base C (17/19);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTT  >  NC_002947/1932491‑1932781
                                                                                                                                                        |                                                                                                                                                 
tGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAATCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGc                                                                                                                                                  >  5:184978/1‑149 (MQ=255)
  gcggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGccc                                                                                                                                                <  8:311204/149‑2 (MQ=255)
             gAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTg                                                                                                                                     >  1:284700/1‑149 (MQ=255)
                gagctacTGCACATGCGACCGCGGCATTGCAGGTCGCACGAAAATCCCAGAAAACGCTGGGGCCGCTTTGCGGCCCTGCGCGGGCACGCCCCCTCCTAAAACGGATATGCGCCCCCGGTCGGAGCGGGCTTGCCCGCGAAGGGCTGCaa                                                                                                                                  <  4:207056/143‑1 (MQ=255)
                gagGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAACCCCGCTACTAAAAAACAAATTCCTCCCCTATAAGAACGGTCTTTCCCGCGAAAGGATgga                                                                                                                                   >  8:260162/1‑146 (MQ=255)
                        ggACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTGCTACAACGGGTATGCGTCCCCTGTTGGAGCGGGGGTGCCCGCGCAGGGCTGGCAGGCAGccc                                                                                                                          >  8:187974/1‑149 (MQ=38)
                            aTTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCCCGCCGACAAAGGATATGCGTCCCCGGTAGGAGCGGGCTTTGCCGCGCAGGGCTTCAAGGCCGCCCCaaa                                                                                                                      >  2:325345/1‑149 (MQ=38)
                                 aCAGCGGCATTGCAGGTCTCACGAAAATAGCAGTACACGCGGGGGCCGCTTTGCGGCCCTTCGCGGGCACGCCCGCTCCTACAAAGGATATGCGGCCCCTGTCGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTc                                                                                                                <  7:70480/149‑1 (MQ=37)
                                            gCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGa                                                                                                     >  5:6449/1‑149 (MQ=255)
                                               ggTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGAt                                                                                                    >  7:175940/1‑147 (MQ=255)
                                                     aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCCAGCCCCCAGCGACAAAGGGTATGCGTGCCCGGTTGGGGCGGGGGTGGCCGCGAAGGGCGGGCAGGGCGCCCCCAAATGGGGCCTACGCGGGTCacggatg                                                                                            >  2:273491/1‑143 (MQ=255)
                                                     aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATc                                                                                            >  7:17674/1‑149 (MQ=255)
                                                       gAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCg                                                                                          >  5:316772/1‑149 (MQ=255)
                                                        aaatTAGCATTAACGGTGGGGGCGTCTTTGGGGCCGTTCGCGGGAACCCCCTCCCCTACACAGTATTTGCGCCCCCTTTAGTCGGGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCgg                                                                                         <  4:150614/149‑1 (MQ=255)
                                                          aaTAGCAGTAAACGGTGGGGCCGCTTTGCGGCCGTGCGGGGGAAAGCCCGCTCCTACAATGTATTTGCGCCCCCTGTTGGTGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCa                                                                                       <  6:6449/149‑1 (MQ=37)
                                                          aaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCgg                                                                                         >  1:288352/1‑147 (MQ=255)
                                                                     aCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCCCTCCAACAAAGGATATGCGTCCCCTGTAGGAGCGCGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCTGATCAGCGATGCGGCAGGGAGATGtt                                                                             >  8:263694/1‑148 (MQ=37)
                                                                            ggCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCa                                                                     <  8:17674/149‑1 (MQ=255)
                                                                                     tgggccccTGCGGGGCCACGCCCGCTCAAACAACGTATATCCGCCCCTTGTAGCCGCGGGCTTCCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGa                                                              <  7:187974/143‑1 (MQ=255)
                                                                                     tGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGATCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATc                                                                <  1:89119/147‑1 (MQ=255)
                                                                                      gCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagc                                                            >  6:66754/1‑149 (MQ=255)
                                                                                       cGGCCCTTCGGGGGAACGCCCGCTCCTACAAAGCATAGGCGCCCCCTGTAGGTGGGGGGTGGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagca                                                           <  1:186422/149‑1 (MQ=255)
                                                                                       cGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagc                                                            <  2:288352/148‑1 (MQ=255)
                                                                                                      aaGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc                                            >  3:112526/1‑149 (MQ=255)
                                                                                                           cGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCa                                                                           <  3:114967/113‑1 (MQ=255)
                                                                                                           cGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCa                                                                           >  4:114967/1‑113 (MQ=255)
                                                                                                                cTACAAAGTATAGGCGTCCCCTGTAGGCGCGGGGTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGa                                    <  1:325345/147‑1 (MQ=255)
                                                                                                                   cAACGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATg                                <  1:273491/148‑1 (MQ=255)
                                                                                                                      aGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCATGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGtt                                                                             >  1:250818/1‑101 (MQ=38)
                                                                                                                      aGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGtt                                                                             <  2:250818/100‑1 (MQ=255)
                                                                                                                            tGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCg                      <  4:112526/149‑1 (MQ=255)
                                                                                                                              cGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGtt                    >  4:320664/1‑149 (MQ=255)
                                                                                                                               gTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTg                   >  8:303164/1‑149 (MQ=255)
                                                                                                                                tCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGc                  <  6:184978/149‑1 (MQ=255)
                                                                                                                                   ccTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGat                                   <  7:153346/129‑1 (MQ=255)
                                                                                                                                   ccTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGat                                   >  8:153346/1‑129 (MQ=255)
                                                                                                                                        aGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGtt           <  5:66754/148‑1 (MQ=255)
                                                                                                                                            gCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTcaca       <  6:316772/148‑1 (MQ=255)
                                                                                                                                             cGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagca                                                           >  2:210850/1‑95 (MQ=255)
                                                                                                                                             cGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagca                                                           <  1:210850/95‑1 (MQ=255)
                                                                                                                                                ggCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATAtt   <  2:284700/149‑1 (MQ=255)
                                                                                                                                                 gCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATAttt  >  6:141291/1‑149 (MQ=255)
                                                                                                                                                        |                                                                                                                                                 
TGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTT  >  NC_002947/1932491‑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: