Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,632,641 G→C 18.3% noncoding (35/76 nt) PP_t59 ← tRNA‑Gly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,632,6410GC18.3% 63.1 / 6.0 27noncoding (35/76 nt)PP_t59tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base G (12/10);  new base C (2/3);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAG  >  NC_002947/4632502‑4632786
                                                                                                                                           |                                                                                                                                                 
cGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGCCAAGGTCggg                                                                                                                                         <  1:434774/150‑3 (MQ=255)
 ggggTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGc                                                                                                                                        >  1:393280/1‑150 (MQ=255)
    gTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGctct                                                                                                                                     <  2:393280/150‑1 (MQ=255)
          aGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAAc                                                                                                                               >  1:42076/1‑150 (MQ=255)
                aaTGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCt                                                                                                                         >  2:99777/1‑150 (MQ=255)
                       tggaacgaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGGCTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCcg                                                                                                                  <  2:281731/148‑1 (MQ=255)
                        cgaacgaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCTGGAAACTAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCcgc                                                                                                                 <  2:273021/150‑1 (MQ=255)
                         gaacgaaAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCcgcg                                                                                                                >  2:619929/1‑150 (MQ=255)
                                       cccGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCGTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGcc                                                                                                  >  1:393064/1‑150 (MQ=255)
                                       cccGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGAACTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGcc                                                                                                  >  1:478785/1‑150 (MQ=255)
                                         cGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCAt                                                                                                <  1:99777/150‑1 (MQ=255)
                                                                  cAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCaa                                                                                                                                <  1:287003/93‑1 (MQ=255)
                                                                  cAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCaa                                                                                                                                >  2:287003/1‑93 (MQ=255)
                                                                       ttCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAAcccc                                                                  >  2:303724/1‑150 (MQ=255)
                                                                        tCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCt                                                                 >  2:68482/1‑150 (MQ=255)
                                                                             tGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGAt                                                                                       >  1:306034/1‑123 (MQ=255)
                                                                             tGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGAt                                                                                       <  2:306034/123‑1 (MQ=255)
                                                                                    tGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCaaaaaa                                                     <  2:393064/150‑1 (MQ=255)
                                                                                                              ggaaaCTTGATCTCGTGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTcgc                           <  2:42076/146‑1 (MQ=255)
                                                                                                                            tagagcacGACCTTGCCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggc                <  2:144022/140‑1 (MQ=255)
                                                                                                                            tagagcacGACCTTGCCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggc                >  1:143988/8‑147 (MQ=255)
                                                                                                                            tagagcacGACCTTGCCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggc                <  2:143988/140‑1 (MQ=255)
                                                                                                                            tagagcacGACCTTGCCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggc                >  1:144022/8‑147 (MQ=255)
                                                                                                                                 cccGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggccaggcggc        <  1:619929/150‑1 (MQ=255)
                                                                                                                                    gACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTgcggccaggcggccag     >  2:563495/1‑150 (MQ=255)
                                                                                                                                       cTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCTAGCCTTCGCGCAGTTCCGGCCTGGCGTCcagcag  >  2:243554/1‑150 (MQ=39)
                                                                                                                                       cTTGGCAAGGTCGTGCTCTACCAACTCAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCcagcag  <  1:303724/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
CGGGGTATTGAGGCAGAATGCATTCGAACGAAAAAAGGCCCCGAAGGGCCTTTTTTCAACGACTTACAGACTTCGACTGAAACCTGTAGATCTTTATTTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCGTCGTGATGGGTGCCATTCTAGCGATTTCTAAAACAGCGTCAACCCCTTGATTCAAAAAACTTTTATTTTTTTTCCGAGCCTTCGCGCAGGTGCGGCCAGGCGGCCAGCAG  >  NC_002947/4632502‑4632786

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