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.C94.1% 104.8 / ‑0.0 34intergenic (‑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 (15/17);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 4.85e-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: Variant not supported by required number of reads on each strand.

AGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTT  >  NC_002947/1932506‑1932781
                                                                                                                                       |                                                                                                                                               
agagGATCTGGACATTCGACAGCGCCATTGCAGGTCTCACTAAAAAAGCAGTAAACGCTGGGGCCGCTTTGGGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCTAAGGGCt                                                                                                                                     <  4:300156/146‑1 (MQ=38)
   ggATCTGGACATTCGACCGCGCAATTCCAGGTCTCACTAAAATACCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCACTCCTACAAAGGATATGTGTCCCCTGTAGTAGCGGGCTTGCCCGCGAAGGGCTGCAAgg                                                                                                                               <  8:142034/149‑1 (MQ=37)
     agccggccATTCGCCACCGCCATAGCAGGCCGCCCCAAAATCCCACTAAACGCCGGGGCCGCATTGCGGCCCGGCGCGGCCACGCCCGCTCCTCCACAGGAAACGCGTCCCCTGTAGCAGCGGGCTGGCCCGCGAAGGGCTGCAAGGCa                                                                                                                             <  8:244270/142‑1 (MQ=255)
             aTTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCATACGCGGCCAAGCCCGATCCAAAAAAGGAATTGCGTCCCCTGTAGGACCGGGGTTGCCCCCGAAGGGCGGAAAGGCAGCCCCaaa                                                                                                                     >  1:217592/1‑149 (MQ=255)
                  aCAGCGGTATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGTATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTc                                                                                                                <  2:70085/149‑1 (MQ=255)
                    aGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAATCCCGCTACAACAAAGGATATGCGTACCCTGTAGGAGCGGGCTTGCCCGCCAAGGGCTGCAACGCAGCCCCAAATTGtag                                                                                                               >  4:239650/1‑146 (MQ=37)
                               aGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGAt                                                                                                    >  2:67441/1‑148 (MQ=255)
                                 gTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTGCGACAAAGGATATGCGGCCCCGGGAGGAGCGGGGGTGGCCGCGCAGGGCTGCAAGGGCGCCCCCAAATTTGGCCTACGCGGGTCAg                                                                                                 >  7:315630/1‑149 (MQ=25)
                                    tcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGGTTGGCCGCGAAGGGCTGCAAGGCAGCCCCAAAATGTc                                                                                                                >  1:132043/1‑131 (MQ=255)
                                           aaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGCATATGCGGCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCa                                                                                       <  3:239650/149‑1 (MQ=255)
                                               gccctAACCGCTGGGGCCGCTTTGCGGCCCTCCGCGGCCAACCCCGCTCCAACAAAGGAAATGCGTCCCCTGTAGGAGCGGCCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGt                                                                                   <  8:342926/145‑1 (MQ=38)
                                                cAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTTTAGACTACGCCGATCAGCGATCCGGCAGGGt                                                                                   >  6:109313/1‑148 (MQ=37)
                                                     aaCGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGtt                                                                             >  3:309110/1‑149 (MQ=255)
                                                          tGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGCCTACGCGGATCAGCGCTCCGGCAGGGTGATGTTCAGCt                                                                        >  3:115943/1‑149 (MQ=255)
                                                          tGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAg                                                                          >  1:54477/1‑147 (MQ=255)
                                                          tGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAg                                                                          <  2:54477/147‑1 (MQ=255)
                                                           ggggCCGCTTGCGGGCCGTGCGCGGCCAACCCGGCTCATACAATGTATATGCGTCCCGTTTCGGAGGGGGTTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTc                                                                       <  3:358435/149‑1 (MQ=255)
                                                               ccGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAg                                                                                      <  7:346584/130‑1 (MQ=255)
                                                               ccGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAAGTGTCGACTACGCGGATCAGCGATCCGGCAg                                                                                      >  8:346584/1‑130 (MQ=255)
                                                                cgtgttGCGGCCCTGCGCGGGCACGCCCGCTCCAACAACGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGaa                                                                  <  8:5339/145‑1 (MQ=38)
                                                                                gcgGGCAAGCCCGCTCCTACAACGGATTTGCGTCCCCTGTCGGAGCGGGCTGGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCt                                                   <  4:115943/148‑1 (MQ=39)
                                                                                gcgGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGTGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTg                                                  <  6:127188/149‑1 (MQ=255)
                                                                                gcgGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCt                                                   <  2:217592/148‑1 (MQ=255)
                                                                                       cgCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTcc                                           <  5:109313/148‑1 (MQ=255)
                                                                                           cGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCt                                       >  2:15815/1‑149 (MQ=255)
                                                                                           cGCTCCTAAAAAGTATATGCGGCCCCTGTAGCAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCt                                       <  8:315630/149‑1 (MQ=39)
                                                                                              tCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGa                                    >  5:164986/1‑149 (MQ=255)
                                                                                                             gCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGt                     >  6:352701/1‑149 (MQ=255)
                                                                                                                tCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTTCAAGGGAGCCCCAAATTTTCGCCGAAGCGGATCAGCGGTGCGGCAGGGGGATGTTGAGGTGCCGAATCGAGGCGCGGCCCTGGGTTTGCACCTCGGTATGCACAAGATCGTTgg                  >  4:8331/1‑148 (MQ=255)
                                                                                                                tCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGc                  >  3:152052/1‑149 (MQ=255)
                                                                                                                    cTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt              >  1:24555/1‑149 (MQ=255)
                                                                                                                         ggAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTc          <  4:121880/148‑1 (MQ=255)
                                                                                                                            gCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACat      <  4:152052/149‑1 (MQ=255)
                                                                                                                                 gCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATAttt  >  3:177023/1‑149 (MQ=255)
                                                                                                                                       |                                                                                                                                               
AGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTT  >  NC_002947/1932506‑1932781

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: