Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,759,516 C→T 10.8% intergenic (+125/+31) trmL → / ← secB tRNA (cytosine(34)‑2'‑O)‑methyltransferase TrmL/protein‑export protein SecB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,759,5160CT10.8% 123.3 / 5.3 46intergenic (+125/+31)trmL/secBtRNA (cytosine(34)‑2'‑O)‑methyltransferase TrmL/protein‑export protein SecB
Reads supporting (aligned to +/- strand):  ref base C (28/13);  new base T (3/2);  total (31/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

GTTTGCCGGGAGCTGATGTGTTGCCTGTGCGGGCCTCTTCGCGGGTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCGGGTCAC  >  NC_002947/5759376‑5759660
                                                                                                                                              |                                                                                                                                                 
gTTTGCCGGGAGCTGATGTGTTGCCTGTGCGGGCCTCTTCGCGGGTAAACCCGCTCCTACGCAGGCCGAGGGCACGGCGCCATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGgc                                                                                                                                           >  2:244286/1‑149 (MQ=255)
     ccGGGAGCTGATGTGTTGCCTGTGCGGGCCTCTTCGCGGGTAAACCCGCTACTAAGCAGGCCGAGGGAACGGCGCAATACCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCtt                                                                                                                                       >  2:30962/1‑148 (MQ=255)
             tgtggtgtTGCCTGGGCGGCCCTCTGCGCGGGAAACCCCGCTCCTCCGCCGCCCGCGGCACCGGCGAAACCCCCTGTCGGAGGGGTTTAACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAAc                                                                                                                              <  6:162164/145‑1 (MQ=255)
             tGATGTGTTGCCTGTGCGGGCCTCTTCGCGGGTAAAACCCCGGCGGCGCAGGGCCAGGGGACCGCGCGATACCCGGGTGGGGGGGGGTTTACCCCGCGGGGGCCCGCGCCCGGCAAAAAAAAACCCCCCGCGACGCGGGTTGATGAAAc                                                                                                                              >  4:82625/1‑149 (MQ=255)
                    ttGCCTGTGCGGGCCTCTTCGCGGGTTAACCCGCGCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTTCCCGCGAAGAGGCCCCCGCCGGGAACAAAAAAGCGCCCCGAGCGGCGCGTTGTTGAAACCGGCAg                                                                                                                        >  6:280634/1‑148 (MQ=255)
                                  ctctTCGCGGGTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCCGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAg                                                                                                          >  6:279236/1‑148 (MQ=255)
                                          gggTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCtt                                                                                                 >  7:269456/1‑149 (MQ=255)
                                                  ccGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGAGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGAGCTTCACCCGCt                                                                                        >  7:113503/1‑149 (MQ=255)
                                                     cTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTcc                                                                                     <  5:279236/149‑1 (MQ=255)
                                                      tCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCt                                                                                    >  6:240781/1‑149 (MQ=255)
                                                            gCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCa                                                                                                                         >  1:169847/1‑107 (MQ=255)
                                                            gCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCa                                                                                                                         <  2:169847/107‑1 (MQ=255)
                                                                gCCGAGGGAACGGCGCAATCCCCTGAAGGAGCGCGTAAACCCGCGAAGAGGCCCCCGCCCGCAACAAAAAAGCGCCCCGACAGGCGCTTACAACAAACCCGAAGAAATTACTTCACAGTAGGCGCTTCACCAGCTTCCTGCATACGCTg                                                                          >  4:125940/1‑149 (MQ=255)
                                                                        aaCGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCtt                                                                  >  3:282271/1‑149 (MQ=255)
                                                                              gcAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCTCTCGGGGCGCTTTTTTg                                                                                                                                  <  6:230122/80‑1 (MQ=21)
                                                                              gcAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCTCTCGGGGCGCTTTTTTg                                                                                                                                  >  5:230122/1‑80 (MQ=21)
                                                                                aaTCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTa                                                           >  1:46139/1‑148 (MQ=255)
                                                                                    ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAg                                                        <  1:244286/149‑1 (MQ=255)
                                                                                     cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAggg                                                      >  2:367041/1‑149 (MQ=255)
                                                                                     cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAgg                                                       >  6:247012/1‑148 (MQ=255)
                                                                                     cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAgg                                                       <  5:247012/148‑1 (MQ=255)
                                                                                         gTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTc                                                  <  4:118071/149‑1 (MQ=255)
                                                                                           aGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCg                                                 >  4:8539/1‑148 (MQ=255)
                                                                                              aGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGt                                             >  2:228719/1‑149 (MQ=255)
                                                                                               gCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGtt                                            >  6:12500/1‑149 (MQ=255)
                                                                                                   gTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGAcc                                        >  3:120135/1‑149 (MQ=255)
                                                                                                         ccGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGa                                   >  2:121760/1‑148 (MQ=255)
                                                                                                         ccGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGa                                   >  1:268535/1‑148 (MQ=255)
                                                                                                         ccGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGAc                                  >  5:106622/1‑149 (MQ=255)
                                                                                                          cgcgAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGt                                             >  2:98105/1‑137 (MQ=255)
                                                                                                          cgcgAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGt                                             <  1:98105/137‑1 (MQ=255)
                                                                                                                       gcgcCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGTCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCgg                    <  3:350449/149‑1 (MQ=255)
                                                                                                                           ccggtttcaaGAAAGCGCCTCTCGGGGCGCTTTTTTGTTGCCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaa                 >  6:222249/8‑149 (MQ=37)
                                                                                                                               caacaaAAAAGCGCCCCGAGAGGCGCTTTCTAGAAACCGCAAGTAATTACTGCACAGTAGGCGCTTAACCCGCTACCTCCATACGCTGCATTTCTTGCGCGGACAGGGCGTAGAAGTTGACCGGCGACAGCATCAGCGCCGGgaacgaa             >  8:117986/1‑149 (MQ=255)
                                                                                                                                aacaaAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaacgaac            >  6:59729/1‑149 (MQ=255)
                                                                                                                                aacaaAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaacgaac            >  4:353261/1‑149 (MQ=255)
                                                                                                                                aacaaAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaacgaac            <  4:120135/149‑1 (MQ=255)
                                                                                                                                aacaaAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaacgaac            <  3:353261/149‑1 (MQ=255)
                                                                                                                                aacaaAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaacgaac            <  2:268535/149‑1 (MQ=255)
                                                                                                                                  caaAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaacgaac            <  6:236025/147‑1 (MQ=255)
                                                                                                                                  caaAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaacgaac            >  5:236025/1‑147 (MQ=255)
                                                                                                                                  caaAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAAGTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACcg          >  6:337387/1‑149 (MQ=255)
                                                                                                                                     gaaaGCGCCTCTCGGGGCGCTTTTTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCAc                                                                                             >  1:126504/2‑63 (MQ=21)
                                                                                                                                     gaaaGCGCCTCTCGGGGCGCTTTTTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCAc                                                                                             <  2:126504/62‑1 (MQ=21)
                                                                                                                                     aaaaGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCgg       <  1:228719/149‑1 (MQ=255)
                                                                                                                                           gcCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGCACCGCGGGTCAc  >  4:254437/1‑148 (MQ=255)
                                                                                                                                             cccGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaa                 <  7:305217/131‑1 (MQ=255)
                                                                                                                                             cccGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGgaa                 >  8:305217/1‑131 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GTTTGCCGGGAGCTGATGTGTTGCCTGTGCGGGCCTCTTCGCGGGTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCGGGTCAC  >  NC_002947/5759376‑5759660

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: