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.C95.8% 151.2 / ‑0.1 48intergenic (‑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 (24/21);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTC  >  NC_002947/1932500‑1932773
                                                                                                                                                |                                                                                                                                        
ctccGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAgg                                                                                                                                >  7:249921/4‑149 (MQ=255)
 ttggAAGAGGCTCTGGACCTTGGACAGCGGGCTTGCAGGTCTGACGAAAATAGCAGTAAACGCGGGGGCCCGGTTGCG‑‑‑GCCCATCGCGGGGAAGCCCGCGGCGACCAAGGTTATGCGGCCCCGGGAGGCGACGGCGTGCCCGCGGAggg                                                                                                                               >  6:82955/4‑149 (MQ=255)
  tCGAAGAGGATCTGGACATTCGACAGTGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTACTACAAAGGATATGCGTCCCCTGTACGAGCGGGCTTGACCGCGAAGGGc                                                                                                                              >  1:204525/1‑149 (MQ=38)
  tCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGc                                                                                                                              >  2:260431/1‑149 (MQ=255)
   cgactcggATCTGCACATTCGACAGCGGCATTCCCGGGCTCACCAAAATAGCCGTAAACGCGGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAggg                                                                                                                               <  8:249921/141‑1 (MQ=25)
    gAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTg                                                                                                                            >  7:250901/1‑149 (MQ=255)
       gagGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGGTATGCGTGCCCAGTAGGAGCGGGCTTGCCCGCGAAGGGCTGc                                                                                                                           >  3:229006/1‑147 (MQ=38)
          gATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAgg                                                                                                                       >  2:47682/1‑148 (MQ=255)
          gATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCGCCTACCAAGGATTTGCGTGCCCTGGAGGAGCGGGCGTGCCCGCGAAGGGCTGCAAGGc                                                                                                                      >  6:224572/1‑149 (MQ=38)
                  cATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCAGTAGGAGCGGGCTTGCACGCGAAGGGCTGCAAGGCAGCCCCa                                                                                                               >  2:279562/1‑148 (MQ=255)
                   aTTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCaaa                                                                                                             >  7:70118/1‑149 (MQ=255)
                   aTTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCaaa                                                                                                             >  2:98301/1‑149 (MQ=255)
                         cAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCg                                                                                                       <  4:119998/149‑1 (MQ=255)
                                   gaaGGTCTCACGAAAAAATCACTAATAACTGGGGCCGCTTTGCCGCATGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTAc                                                                                                  <  1:126482/146‑1 (MQ=21)
                                   gCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGa                                                                                             >  8:24960/1‑149 (MQ=255)
                                       gTCTCACTAAAATACCAGTAACCGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAACGTATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAg                                                                                         <  2:279413/149‑1 (MQ=37)
                                       gTCTCACTAAAATACCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGTATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAg                                                                                         <  7:28333/149‑1 (MQ=37)
                                         ctcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTc                                                                                                        <  5:89976/132‑1 (MQ=255)
                                         ctcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCGGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTc                                                                                                        >  6:89976/1‑132 (MQ=255)
                                          tcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTACTACAAAGGATATTCGTGCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGCCCCAAATTGTGGACTACGCGCAtacacg                                                                                       >  8:164116/1‑143 (MQ=25)
                                            aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATc                                                                                    >  8:156199/1‑149 (MQ=255)
                                            aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATc                                                                                    >  6:156311/1‑149 (MQ=255)
                                            aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGcccgcga                                                                                                                                  <  5:83823/102‑6 (MQ=255)
                                            aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGcccgcga                                                                                                                                  >  6:83823/1‑97 (MQ=255)
                                              gAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCg                                                                                  >  5:274482/1‑149 (MQ=255)
                                                     gCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGt                                                                           <  8:250901/149‑1 (MQ=255)
                                                          aaaCGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATg                                                                       <  8:70118/148‑1 (MQ=255)
                                                           aaCGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTCCAAAGCATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATg                                                                       <  7:164116/147‑1 (MQ=255)
                                                             ggcgggggCGGTTTTGCG‑‑‑GCCCTGCGCGGCCACGCCCGCTCCTAAAAAGTATCTGCGTCCCGTGTAGGCGGGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCa                                                                   <  7:127780/145‑1 (MQ=255)
                                                              gCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAg                                                                  >  7:223466/1‑149 (MQ=255)
                                                               cTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAACCCGCAACGAAAAAAGATATGCGTGCCCGGTAAGAGCGGGCTTTGCCGCCCAGGGCAGCAAAGGAGCCCCAAAATTTAGACAAAGCGGGAGCGCGAAACCGCAGGGTGATGttgat                                                                  >  2:334621/1‑145 (MQ=255)
                                                                    gCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTc                                                               >  5:88910/1‑146 (MQ=255)
                                                                           tGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCg                                                       <  1:260431/147‑1 (MQ=255)
                                                                           tGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGa                                                      <  1:98301/148‑1 (MQ=255)
                                                                           tGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGa                                                      <  6:274482/148‑1 (MQ=255)
                                                                            gCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGag                                                     >  6:83772/1‑148 (MQ=255)
                                                                              cG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagc                                                    <  5:83772/148‑1 (MQ=255)
                                                                                    cccTTCGCGGGCAAGCCCGCGCCTACAAAGGATATGCGTCCCCGGTAGGAGCGGGCGTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCt                                                >  4:269839/1‑149 (MQ=39)
                                                                                                  gCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCa                                  <  1:47682/149‑1 (MQ=255)
                                                                                                   cccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAg                                 >  6:228913/1‑149 (MQ=255)
                                                                                                   cccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAg                                 >  4:299921/1‑149 (MQ=255)
                                                                                                     cGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCt                               <  3:299921/149‑1 (MQ=255)
                                                                                                          cTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGata                          >  3:113700/1‑149 (MQ=255)
                                                                                                                    tatGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatcat                <  4:113700/149‑1 (MQ=255)
                                                                                                                     atGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATc                                                        >  2:350513/1‑108 (MQ=255)
                                                                                                                     atGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATc                                                        <  1:350513/108‑1 (MQ=255)
                                                                                                                      tGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACatcatc               <  5:228913/148‑1 (MQ=255)
                                                                                                                              cTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt      >  8:41361/1‑149 (MQ=255)
                                                                                                                                   ggAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTc  <  2:240834/148‑1 (MQ=255)
                                                                                                                                   ggAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTc  <  5:156311/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                        
CTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCG‑‑‑GCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTC  >  NC_002947/1932500‑1932773

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: