Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 536,283 A→G 13.9% noncoding (178/190 nt) PP_mr06 ← P27 ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947536,2830AG13.9% 82.2 / 6.1 36noncoding (178/190 nt)PP_mr06P27 ncRNA
Reads supporting (aligned to +/- strand):  ref base A (18/13);  new base G (2/3);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01

TGAAGCCCTGATCGCTGACCTGAAGCGTATCAAGCCAGCTTCCTCGAAAGGTATCTACGTTAAGCGCGTTACCCTGAGCACCACCATGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAGCCCACGCAGATGGTGCTCCCGATTCGTTACCGAATC  >  NC_002947/536135‑536401
                                                                                                                                                    |                                                                                                                      
tGAAGCCCTGATCGCTGACCTGAAGCGTATCAAGCCAGCTTCCTCGAAAGGTATCTACGTTAAGCGCGTTACCCTGAGCACCACCATGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAg                                                                                                                       >  1:1170080/1‑150 (MQ=255)
       cTGATCGCTGACCTGAAGCGTATCAAGCCAGCTTCCTCGAAAGGTATCTACGTTAAGCGCGTTACCCTGAGCACCACCATGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACTGGTCGcgcc                                                                                                                <  2:1294143/150‑1 (MQ=255)
       cTGATCGCTGACCTGAAGCGTATCAAGCCAGCTTCCTCGAAAGGTATCTACGTTAAGCGCGTTACCCTGAGCACCACCATGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGcgcc                                                                                                                <  1:648934/150‑1 (MQ=255)
        tGATCGCTGACCTGAAGCGTATCAAGCCAGCTTCCTCGAAAGGTATCTACGTTAAGCGCGTTACCCTGAGCACCACCATGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGcgccg                                                                                                               <  2:196944/150‑1 (MQ=255)
                        gCGTATCAAGCCAGCTTCCTCGAAAGGTATCTACGTTAAGCGCGTTACCCTGAGCACCACCATGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTgggg                                                                                               >  2:1630372/1‑150 (MQ=255)
                                               aaGGTATCTACGTTAAGCGCGTTACCCTGAGCACCACCATGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATcc                                                                        >  1:867839/1‑150 (MQ=255)
                                                          gTTAAGCGCGTTACCCTGAGCACCACCATGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAgg                                                             <  1:324142/150‑1 (MQ=255)
                                                                 gcgTTACCCTGAGCACCACCATGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCcct                                                                                          >  1:1313022/1‑114 (MQ=255)
                                                                 gcgTTACCCTGAGCACCACCATGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCcct                                                                                          <  2:1313022/114‑1 (MQ=255)
                                                                   gTTACCCTGAGCACCACCATGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAg                                                    >  2:864356/1‑150 (MQ=255)
                                                                                   ccaTGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGcgccgc                                                                                                              <  1:1167826/76‑1 (MQ=255)
                                                                                   ccaTGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGcgccgc                                                                                                              >  2:1167826/1‑76 (MQ=255)
                                                                                   ccaTGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAAtt                                                                                                   >  2:1073220/1‑87 (MQ=255)
                                                                                   ccaTGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAAtt                                                                                                   <  1:1073220/87‑1 (MQ=255)
                                                                                   ccaTGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAGcc                                    >  1:1449349/1‑150 (MQ=255)
                                                                                         gCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAGCCCACGCa                              >  2:1303945/1‑150 (MQ=255)
                                                                                          cccAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAg                                      <  2:886385/141‑1 (MQ=255)
                                                                                          cccAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAg                                      >  1:886385/1‑141 (MQ=255)
                                                                                             aGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAGCCCACGCAGATg                          >  1:1506317/1‑150 (MQ=255)
                                                                                               gTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGCCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAGCCCACGCAGATGGt                        >  2:313991/1‑150 (MQ=255)
                                                                                                       aTCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTggg                                                                                                <  2:535997/70‑1 (MQ=255)
                                                                                                       aTCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTggg                                                                                                >  1:535997/1‑70 (MQ=255)
                                                                                                       aTCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAGCCCACGCAGATGGTGCTCCCGa                >  2:1233898/1‑150 (MQ=255)
                                                                                                            tCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAGCCCACGCAGATGGTGCTCCCGATTCGt           <  2:1506317/150‑1 (MQ=255)
                                                                                                                    cGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAGCCCACGCAGATGGTGCTCCCGATTCGTTACCGAAt   <  1:1630372/150‑1 (MQ=255)
                                                                                                                     gCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAggcg                                                           <  2:1375672/93‑1 (MQ=255)
                                                                                                                     gCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAggcg                                                           >  1:1375672/1‑93 (MQ=255)
                                                                                                                     gCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAGCCCACGCAGATGGTGCTCCCGATTCGTTACCGAATc  >  1:380114/1‑150 (MQ=255)
                                                                                                                             tgtAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAGCCCACGCAGATg                          >  1:1526832/1‑118 (MQ=255)
                                                                                                                             tgtAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAGCCCACGCAGATg                          <  2:1526832/118‑1 (MQ=255)
                                                                                                                                  gAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCggggg                                                                              <  2:842521/61‑1 (MQ=255)
                                                                                                                                  gAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCggggg                                                                              >  1:842521/1‑61 (MQ=255)
                                                                                                                                         gcggcgCGACTGGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAg                                      >  1:354034/1‑94 (MQ=255)
                                                                                                                                         gcggcgCGACTGGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAg                                      <  2:354034/94‑1 (MQ=255)
                                                                                                                                           ggcgCGACTGGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCa                                       >  1:341138/1‑91 (MQ=255)
                                                                                                                                           ggcgCGACTGGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCa                                       <  2:341138/91‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                      
TGAAGCCCTGATCGCTGACCTGAAGCGTATCAAGCCAGCTTCCTCGAAAGGTATCTACGTTAAGCGCGTTACCCTGAGCACCACCATGGGCCCAGGTCTGGTCATCGATCAGAGCTCGCTGAACGTGTAAGAACATGGCGGCGCGACCAGTCGCGCCGCCAGCAAAAATTGGGGTCCCTGCCTGGCGGGGGCTATCCAAGACCGTAGGCGGCGCAAGCCTTAAAACATCAGCCCACGCAGATGGTGCTCCCGATTCGTTACCGAATC  >  NC_002947/536135‑536401

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