Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,079,235 C→A 11.4% noncoding (42/77 nt) PP_t56 ← tRNA‑Arg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,079,2350CA11.4% 92.9 / 3.3 35noncoding (42/77 nt)PP_t56tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base C (15/16);  new base A (2/2);  total (17/18)
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.97e-01

CTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCTTTTCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGC  >  NC_002947/5079090‑5079377
                                                                                                                                                 |                                                                                                                                              
cTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCAGCCTCAGCAGGATTCGAACCTGCGACCGCTCGGtt                                                                                                                                            >  1:114907/1‑150 (MQ=255)
         acacCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGa                                                                                                                                   <  2:1003826/150‑1 (MQ=255)
                     ttACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGACGCCGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCa                                                                                                                       >  1:900202/1‑150 (MQ=255)
                        cAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                                                                                                    <  1:296869/150‑1 (MQ=255)
                                  tAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAg                                                                                                          <  1:842089/150‑1 (MQ=255)
                                                  aGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGttt                                                                                          <  1:426404/150‑1 (MQ=255)
                                                                   aaacaaaAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACt                                                                         >  2:558992/1‑150 (MQ=255)
                                                                   aaacaaaAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACt                                                                         >  2:559068/1‑150 (MQ=255)
                                                                     acaaaAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTgg                                                                       <  2:865055/150‑1 (MQ=255)
                                                                     acaaaAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTgg                                                                       <  1:440589/150‑1 (MQ=255)
                                                                       aaaaaGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGTGGTTTTTAGCCAAGTTACCGCTGGCt                                                                     >  1:421198/1‑150 (MQ=18)
                                                                            gCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                                                                                             >  2:122662/1‑105 (MQ=255)
                                                                            gCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                                                                                             <  1:122662/105‑1 (MQ=255)
                                                                             cccGCGCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATGGTAGTTTTAGAACAGGGTACCACTGGCTTACTGa                                                               >  1:826292/1‑150 (MQ=255)
                                                                              ccGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGaa                                                              >  2:644027/1‑150 (MQ=255)
                                                                              ccGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGGTTTAGACCAGGTTAACACTGGCTTACCGaa                                                              >  2:657817/1‑150 (MQ=255)
                                                                                              gCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTAGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGc                                                         <  2:1059603/139‑1 (MQ=255)
                                                                                              gCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTAGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGc                                                         <  2:1059404/139‑1 (MQ=255)
                                                                                              gCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTAGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGc                                                         >  1:1059603/1‑139 (MQ=255)
                                                                                              gCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTAGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGc                                                         >  1:1059404/1‑139 (MQ=255)
                                                                                                tttttGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                                                                                             >  2:1091149/1‑85 (MQ=255)
                                                                                                tttttGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                                                                                             <  1:1091149/85‑1 (MQ=255)
                                                                                                tttttGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCTTTTCAAGCAGaa                                            >  1:920346/1‑150 (MQ=255)
                                                                                                tttttGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCTTTTCAAGCAGaa                                            >  1:629315/1‑150 (MQ=255)
                                                                                                        ggggTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCTTTTCAAGCAGAACCACTTCa                                    <  1:559068/150‑1 (MQ=255)
                                                                                                        ggggTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCTTTTCAAGCAGAACCACTTCa                                    <  1:558992/150‑1 (MQ=255)
                                                                                                         gggTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCTTTTCAAGCAGAACCACTTCaa                                   <  2:920346/150‑1 (MQ=255)
                                                                                                           gTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCtttt                                                     <  1:732325/130‑1 (MQ=255)
                                                                                                           gTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCtttt                                                     >  2:732325/1‑130 (MQ=255)
                                                                                                                gCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCTTTTCAAGCAGAACCACTTCAAGAAAGTg                            <  2:629315/150‑1 (MQ=255)
                                                                                                                  ggACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCTTTTCAAGCAGAACCACTTCAAGAAAGTGGc                          <  2:900202/150‑1 (MQ=255)
                                                                                                                        agcagGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCTTTTCAAGCAGAACCACTTCAAGAAAGTGGCGGACTc                    >  2:693359/1‑150 (MQ=255)
                                                                                                                            ggATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCTTTTCAAGCAGAACCACTTCAAGAAAGTGGCGGACTcagca                <  1:644027/150‑1 (MQ=255)
                                                                                                                                  gAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCTTTTCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTc          >  2:968917/1‑150 (MQ=255)
                                                                                                                                          gACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCTTTTCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGc  <  1:693359/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAATGGCTTTTCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGC  >  NC_002947/5079090‑5079377

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