Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,960,051 A→T 10.2% noncoding (74/103 nt) PP_mr63 ← rnk_pseudo ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,960,0510AT10.2% 116.6 / 2.9 39noncoding (74/103 nt)PP_mr63rnk_pseudo ncRNA
Reads supporting (aligned to +/- strand):  ref base A (21/14);  new base T (2/2);  total (23/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 = 9.94e-01

CGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGT  >  NC_002947/5959916‑5960195
                                                                                                                                       |                                                                                                                                                
cgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTga                                                                                                                                    <  2:527350/150‑1 (MQ=255)
cgAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTga                                                                                                                                    <  1:407871/150‑1 (MQ=255)
     cTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAgg                                                                                                                               <  2:156322/150‑1 (MQ=255)
          gATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTg                                                                                                                          <  1:572305/150‑1 (MQ=255)
            tCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATCAGAAATAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGc                                                                                                                        >  1:8521/1‑150 (MQ=255)
                                   tCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACACCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGacac                                                                                                 >  1:178561/1‑150 (MQ=255)
                                   tCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAAAAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCTCGCAAGGGCGGGGTGTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGacac                                                                                                 >  2:736330/1‑150 (MQ=255)
                                                 gATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGc                                                                                   >  2:332462/1‑150 (MQ=255)
                                                  aTGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGttt                                                                                                                                      >  1:508840/1‑98 (MQ=255)
                                                  aTGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGttt                                                                                                                                      <  2:508840/98‑1 (MQ=255)
                                                    gAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCTTGGGCggg                                                                                                                                          >  1:557859/1‑92 (MQ=255)
                                                    gAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCTTGGGCggg                                                                                                                                          <  2:557859/92‑1 (MQ=255)
                                                         aaGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGGTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTc                                                                           >  2:831784/1‑150 (MQ=255)
                                                                    cTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGc                                                                                                                        <  1:179958/94‑1 (MQ=255)
                                                                    cTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGc                                                                                                                        >  2:179958/1‑94 (MQ=255)
                                                                    cTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTgcg                                                                >  1:671833/1‑150 (MQ=255)
                                                                     tCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTgcgc                                                               >  2:69604/1‑150 (MQ=255)
                                                                                      tCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCt                                              >  1:588436/1‑150 (MQ=255)
                                                                                              aCAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGcc                                                     >  1:719303/1‑135 (MQ=255)
                                                                                              aCAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGcc                                                     <  2:719303/135‑1 (MQ=255)
                                                                                                aGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCaa                                    >  1:837704/1‑150 (MQ=255)
                                                                                                        ttCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTggg                            <  2:667919/150‑1 (MQ=255)
                                                                                                               tgCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTcgcg                     >  1:81606/1‑150 (MQ=255)
                                                                                                                gCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATCCACGATCCCCTGCGCTAAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCgg                    >  1:225857/1‑150 (MQ=255)
                                                                                                                 cAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTAc                                                                                              >  2:419803/1‑75 (MQ=255)
                                                                                                                 cAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTAc                                                                                              <  1:419803/75‑1 (MQ=255)
                                                                                                                    aTAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGa                >  2:436521/1‑150 (MQ=255)
                                                                                                                       aaGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCt             <  2:837704/150‑1 (MQ=255)
                                                                                                                          cAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGcc          >  2:676571/1‑150 (MQ=255)
                                                                                                                            aaaCCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCgg        >  2:20984/1‑150 (MQ=255)
                                                                                                                             aaCCCCGCCCTTGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCt             <  1:829995/144‑1 (MQ=255)
                                                                                                                             aaCCCCGCCCTTGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCt             >  2:829995/1‑144 (MQ=255)
                                                                                                                              acaccGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCgcga      <  2:178561/147‑4 (MQ=255)
                                                                                                                               ccccGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGc           <  2:586836/144‑1 (MQ=255)
                                                                                                                               ccccGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGc           >  1:586836/1‑144 (MQ=255)
                                                                                                                               ccccGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCgcgag     <  1:332462/150‑5 (MQ=255)
                                                                                                                               ccccGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGgcgc     <  2:588436/150‑1 (MQ=255)
                                                                                                                                cccGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGcc    >  1:113622/1‑150 (MQ=255)
                                                                                                                                  cGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGt  >  1:495265/1‑150 (MQ=255)
                                                                                                                                       |                                                                                                                                                
CGAGGCTGTCGATCAGACGCTCCAGACGCTGTACGTCCAATCGGGTGAGGATGAGGGAAGGCTTGGTGCTCATGATCCAGTCAGACTCCTTTGAACAGGCAGTTTTCAATGTGCAAATAAAGCAAAACCCCGCCCAAGGGCGGGGTTTGAGAAGGCTTTGGCGTCACCGACGCAGAACCCGACACTACCACAGCCGGGCAAATACTCAAGCCCCTGCGCTCAGTGTGCCTTAGCCTGCTCGCGCAAGACCTGGGCCTCGCGGCAGATCTGCCGGCGCCGT  >  NC_002947/5959916‑5960195

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: ‑