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.C97.8% 147.5 / ‑3.0 45intergenic (‑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 (22/22);  total (23/22)
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.84e-01
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.

GCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTT  >  NC_002947/1932493‑1932781
                                                                                                                                                        |                                                                                                                                                 
gcggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACCAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGccc                                                                                                                                                >  5:289612/1‑148 (MQ=255)
  ggcgTCTTGGTAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGAAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGcccgc                                                                                                                                              <  2:328315/149‑4 (MQ=255)
     gTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCTGGGGGCGAAACAACAAAGGATATGCGTCCCCTGAAGCACCGGGCTTGccc                                                                                                                                                >  1:136412/1‑148 (MQ=255)
      tCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAg                                                                                                                                           <  2:87869/149‑1 (MQ=255)
           gAATAGGATGTGCACATCCGACAGGGGCATTGCAGGTCTCACTAAAATCGCAGTAACCGCTGGGGCGGCTTTGGGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAACGGATATGCGGCCCCTGTCGGAGCGGGCTTGCCCGCGAAGGGCTg                                                                                                                                      <  4:225976/149‑1 (MQ=25)
                ggATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAACC‑‑‑‑‑CCGCTCCTACAAAGCATATGCCTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAgg                                                                                                                                 <  1:247718/149‑1 (MQ=38)
                   tCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCACCTAAAAAAGATATTCGTACCCTGTATGAGCCGGCGTTCCCGCGAAGGGCTGCAATGCAg                                                                                                                              >  2:332530/1‑149 (MQ=255)
                      ggACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGccc                                                                                                                          >  5:344367/1‑149 (MQ=255)
                                 aGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTTCAAAGGATATGCGTCCCCAGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAAATGTCGa                                                                                                              >  2:90750/1‑149 (MQ=37)
                                           cAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGCCTACGCGGAt                                                                                                    >  2:275107/1‑149 (MQ=255)
                                               tctcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGc                                                                                                <  3:117735/149‑1 (MQ=255)
                                                 tcACTAAAATAGCAGAAAAGGGTGGGGCGGTTTGGGGGCCGTGGGGGGGCAAGC‑‑‑‑‑CGGCTCTTACAAAGTATTTGCGCCCCCTGTTGGTGGGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAg                                                                                                 <  4:186282/146‑1 (MQ=255)
                                                 tcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCGCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCCCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGGCTACGCGGCTCAg                                                                                                 >  3:186282/1‑146 (MQ=37)
                                                    ctaaaaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCGTTCGCGGGAAAGC‑‑‑‑‑CCGCTCATACAAAGTATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATcc                                                                                           <  6:344367/147‑1 (MQ=38)
                                                     gAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCaa                                                                                                                                   >  4:73288/1‑110 (MQ=255)
                                                     gAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCaa                                                                                                                                   <  3:73288/110‑1 (MQ=255)
                                                        aaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTTTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCa                                                                                       <  4:57896/149‑1 (MQ=255)
                                                             cAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGAAATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGt                                                                                   >  6:294002/1‑148 (MQ=255)
                                                               gTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGAt                                                                                >  3:185545/1‑149 (MQ=255)
                                                               gTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGccc                                                                                                                                                <  7:65463/86‑2 (MQ=255)
                                                               gTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGccc                                                                                                                                                >  8:65463/1‑85 (MQ=255)
                                                                tAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGcc                                                                                                                           >  7:35005/1‑106 (MQ=255)
                                                                tAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTAAAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGcc                                                                                                                           <  8:35005/106‑1 (MQ=255)
                                                                    cggtGGGGCCCGCTTGCGGCCCTTGGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCGGTAGGCGGGGGCTTGCCCGCGAATGGCTGCAAGGCAGCCCCAAATTGTCTGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTc                                                                            <  1:275107/146‑1 (MQ=255)
                                                                    ccgtGGGGCCGCTTTGCGGCCCTTGGCGGGCAAGC‑‑‑‑‑CCGCTCGTTCAAAGGTTATCCGTCGCCTGTAGGAGCGGGGTTGCCCGCGAAGGGCTGCAAGGCAGCCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTc                                                                            <  1:90750/146‑1 (MQ=255)
                                                                        ggggCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATACGGCAGGGTGATGTTCAGCTc                                                                       >  6:322331/1‑149 (MQ=255)
                                                                         gggCCGCTTTTCGGCCCTTCGCTGGCAACC‑‑‑‑‑CCGCTCCTTCAAAGCATCTCCGGCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCt                                                                        <  6:345410/147‑1 (MQ=37)
                                                                           gCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCGGCTACAAAGGATATGGGTCCCCTGTAGGAGCGGGCTTGCCCGCGACGGGCTGCAAGGCAGCCCCAAAATGTc                                                                                                                 >  7:102026/1‑105 (MQ=25)
                                                                                  tGCGGCCGTTCGCTGGAAAGC‑‑‑‑‑CCGTTCCTACAAAGGATATGCGTCCCCTGTAGGTGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGa                                                              <  2:136412/148‑1 (MQ=39)
                                                                                   gCGGCCAGTAGCTGGAAAGT‑‑‑‑‑CCTCTCCTACAAAGGATTTGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGag                                                             <  6:289612/148‑1 (MQ=37)
                                                                                          ttCGCGGGCAAGC‑‑‑‑‑CCGCTCCTAAAAACGAAATGCGTCCCCAGTAGGAGCCGGCTTGCCCGCCAAGGGCTGCAACGCAGCACCAAATTGACGACTACGCGGATCAGCGCTCCGGGAGGGTGATGTTGAGATCCAGAATCGagcagcgtca                                                     >  4:307847/1‑145 (MQ=255)
                                                                                           gagcgGGCTTGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAg                                                                                      >  4:328070/3‑115 (MQ=38)
                                                                                           gagcgGGCTTGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAg                                                                                      <  3:328070/113‑1 (MQ=38)
                                                                                               gggCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATc                                                                                                   <  5:403003/98‑1 (MQ=255)
                                                                                               gggCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATc                                                                                                   >  6:403003/1‑98 (MQ=255)
                                                                                                   aaGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc                                            <  5:294002/149‑1 (MQ=255)
                                                                                                                tCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGt                                                                              <  5:80851/107‑1 (MQ=255)
                                                                                                                tCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGt                                                                              >  6:80851/1‑107 (MQ=255)
                                                                                                                   tACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                                 >  5:235947/1‑149 (MQ=255)
                                                                                                                        aGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGa                                    <  1:145133/141‑1 (MQ=255)
                                                                                                                        aGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGa                                    >  2:145133/1‑141 (MQ=255)
                                                                                                                               gCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATcc                                                                                           <  2:352124/79‑1 (MQ=255)
                                                                                                                               gCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATcc                                                                                           >  1:352124/1‑79 (MQ=255)
                                                                                                                                 gTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTg                   >  4:68235/1‑149 (MQ=255)
                                                                                                                                 gTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTg                   >  5:6940/1‑149 (MQ=255)
                                                                                                                                         tAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCa                                                                           >  6:384505/1‑85 (MQ=255)
                                                                                                                                         tAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCa                                                                           <  5:384505/85‑1 (MQ=255)
                                                                                                                                           ggAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTca         <  4:163303/149‑1 (MQ=255)
                                                                                                                                            gAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTca         >  4:122430/1‑148 (MQ=255)
                                                                                                                                                  gCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATAttt  >  2:385475/1‑149 (MQ=255)
                                                                                                                                                        |                                                                                                                                                 
GCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTT  >  NC_002947/1932493‑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: