Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,759,521 A→G 10.6% intergenic (+130/+26) 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,5210AG10.6% 126.5 / 5.5 48intergenic (+130/+26)trmL/secBtRNA (cytosine(34)‑2'‑O)‑methyltransferase TrmL/protein‑export protein SecB
Reads supporting (aligned to +/- strand):  ref base A (28/14);  new base G (3/2);  total (32/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.03e-01

GTTTGCCGGGAGCTGATGTGTTGCCTGTGCGGGCCTCTTCGCGGGTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCGGGTCACCAGGCTG  >  NC_002947/5759376‑5759667
                                                                                                                                                   |                                                                                                                                                   
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                                                                                                                                         >  5:230122/1‑80 (MQ=21)
                                                                              gcAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCTCTCGGGGCGCTTTTTTg                                                                                                                                         <  6:230122/80‑1 (MQ=21)
                                                                                aaTCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTa                                                                  >  1:46139/1‑148 (MQ=255)
                                                                                    ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAg                                                               <  1:244286/149‑1 (MQ=255)
                                                                                     cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAggg                                                             >  2:367041/1‑149 (MQ=255)
                                                                                     cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAgg                                                              <  5:247012/148‑1 (MQ=255)
                                                                                     cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAgg                                                              >  6:247012/1‑148 (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                                          >  1:268535/1‑148 (MQ=255)
                                                                                                         ccGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGa                                          >  2:121760/1‑148 (MQ=255)
                                                                                                         ccGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGAc                                         >  5:106622/1‑149 (MQ=255)
                                                                                                          cgcgAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGt                                                    <  1:98105/137‑1 (MQ=255)
                                                                                                          cgcgAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGt                                                    >  2:98105/1‑137 (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)
                                                                                                                                                 agagGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCGGGTCACCAGGCTg  <  5:280634/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTTTGCCGGGAGCTGATGTGTTGCCTGTGCGGGCCTCTTCGCGGGTAAACCCGCTCCTACGCAGGCCGAGGGAACGGCGCAATCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCGCCGGCAACAAAAAAGCGCCCCGAGAGGCGCTTTCTTGAAACCGGCAGTAATTACTGCACAGTAGGCGCTTCACCCGCTTCCTGCATACGCTGCATTTCTTGCGCGTACAGGGCGTCGAAGTTGACCGGCGACAGCATCAGCGCCGGGAACGAACCGCGGGTCACCAGGCTG  >  NC_002947/5759376‑5759667

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: