Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,981,977 A→T 8.7% noncoding (28/207 nt) PP_mr43 ← Cobalamin regulatory region

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,981,9770AT 8.7% 126.7 / 3.2 46noncoding (28/207 nt)PP_mr43Cobalamin regulatory region
Reads supporting (aligned to +/- strand):  ref base A (21/21);  new base T (2/2);  total (23/23)
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.06e-01

TGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGCCCCC  >  NC_002947/3981830‑3982122
                                                                                                                                                   |                                                                                                                                                 
tGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAg                                                                                                                                                  >  3:188019/1‑149 (MQ=255)
 ggCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAgg                                                                                                                                                 >  6:349800/1‑149 (MQ=255)
        aCCGTTACATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCgaag                                                                                                                                           >  6:17463/1‑148 (MQ=255)
          cGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAAc                                                                                                                                        >  1:305043/1‑149 (MQ=255)
            ttCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCt                                                                                                                                      >  6:310772/1‑149 (MQ=255)
                  gTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGt                                                                                                                                >  7:220119/1‑149 (MQ=255)
                     cccGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGc                                                                                                                             <  2:165917/149‑1 (MQ=255)
                       cGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCaa                                                                                                                           <  8:289146/149‑1 (MQ=255)
                         ttCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAgg                                                                                                                         <  4:188019/149‑1 (MQ=255)
                                  gCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGt                                                                                                                <  2:305043/149‑1 (MQ=255)
                                   cGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGtg                                                                                                               <  5:17463/149‑1 (MQ=255)
                                         gACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGtggctggc                                                                                                         >  4:183902/1‑149 (MQ=255)
                                           cGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGGCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCgg                                                                                                       <  5:310772/149‑1 (MQ=255)
                                            gCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCTGGCCgaagc                                                                                                                                          >  5:98465/1‑112 (MQ=255)
                                            gCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCTGGCCgaagc                                                                                                                                          <  6:98465/113‑2 (MQ=255)
                                              cTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGggtg                                                                                                    <  5:349800/149‑1 (MQ=255)
                                              cTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGggtg                                                                                                    <  7:130079/149‑1 (MQ=255)
                                              cTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGggtg                                                                                                    >  7:134456/1‑149 (MQ=255)
                                                  ccGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTc                                                                                                 >  3:130596/1‑148 (MQ=255)
                                                    gcgcgTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAAc                                                                                              >  4:137856/1‑149 (MQ=255)
                                                            gcgcAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATtgct                                                                                       <  8:10372/148‑1 (MQ=255)
                                                              gcAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCg                                                                                    <  3:137856/149‑1 (MQ=255)
                                                                  tGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGc                                                                                <  3:183902/149‑1 (MQ=255)
                                                                  tGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGc                                                                                <  8:220119/149‑1 (MQ=255)
                                                                     cTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCgg                                                                             <  8:84659/149‑1 (MQ=255)
                                                                      tgtgTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGcc                                                                                                                                               >  8:154943/1‑82 (MQ=255)
                                                                      tgtgTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGcc                                                                                                                                               <  7:154943/82‑1 (MQ=255)
                                                                            gATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACaa                                                                      <  2:260557/149‑1 (MQ=255)
                                                                            gATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACa                                                                       <  4:130596/148‑1 (MQ=255)
                                                                                cGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGGTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGGGGGTGGTCAACCATTGCTGCGGGGCCGGTCGCCAACGCt                                                                  >  2:317564/1‑149 (MQ=255)
                                                                                 gTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGAGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTg                                                                 >  7:113653/1‑149 (MQ=255)
                                                                                            ttACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGc                                                      >  2:142149/1‑149 (MQ=255)
                                                                                                    gCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGt                                                        <  4:211648/139‑1 (MQ=255)
                                                                                                    gCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGt                                                        >  3:211648/1‑139 (MQ=255)
                                                                                                       ggggCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACa                                                >  4:12669/1‑144 (MQ=255)
                                                                                                       ggggCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACa                                                <  3:12669/144‑1 (MQ=255)
                                                                                                         ggCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGgcgc                                         >  6:59040/1‑149 (MQ=255)
                                                                                                              ccGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTAt                                    >  1:33312/1‑149 (MQ=255)
                                                                                                                gcgGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATc                                   >  8:111795/1‑148 (MQ=255)
                                                                                                                     ttGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGt                             >  1:4858/1‑149 (MQ=255)
                                                                                                                       gTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGt                             <  7:111795/147‑1 (MQ=255)
                                                                                                                         aCCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCaa                          <  8:113653/148‑1 (MQ=255)
                                                                                                                                  ccGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCa                >  6:90133/1‑149 (MQ=255)
                                                                                                                                        gttcTTCGGCCTGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGa                       <  2:44482/133‑1 (MQ=255)
                                                                                                                                        gttcTTCGGCCTGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGa                       >  1:44482/4‑136 (MQ=255)
                                                                                                                                                gCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGccccc  <  5:59040/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGGCTTGCACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGCCCCC  >  NC_002947/3981830‑3982122

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