Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,981,977 A→T 6.9% 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 6.9% 162.4 / 2.9 58noncoding (28/207 nt)PP_mr43Cobalamin regulatory region
Reads supporting (aligned to +/- strand):  ref base A (26/28);  new base T (2/2);  total (28/30)
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

CACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGCCCCCGC  >  NC_002947/3981837‑3982124
                                                                                                                                            |                                                                                                                                                   
cACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCgaa                                                                                                                                              <  3:11616/148‑1 (MQ=255)
 aCCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCgaag                                                                                                                                             >  1:213989/1‑148 (MQ=255)
     ttCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCTGGCCgaagc                                                                                                                                            >  5:218020/1‑144 (MQ=255)
     ttCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCTGGCCgaagc                                                                                                                                            <  6:218020/145‑2 (MQ=255)
         aTGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGaa                                                                                                                                    <  3:195299/149‑1 (MQ=255)
             ccccGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTg                                                                                                                                 >  6:258363/1‑148 (MQ=255)
              cccGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTgg                                                                                                                                >  6:80611/1‑148 (MQ=255)
                     gggCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTa                                                                                                                        >  1:345393/1‑149 (MQ=255)
                             gggATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGtgg                                                                                                                >  4:234955/1‑149 (MQ=255)
                                  gACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGtggctggc                                                                                                           >  6:42761/1‑149 (MQ=255)
                                    cGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCgg                                                                                                         >  6:61793/1‑149 (MQ=255)
                                     gCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCggg                                                                                                        >  5:287848/1‑149 (MQ=255)
                                      ccTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGggt                                                                                                       <  3:234955/149‑1 (MQ=255)
                                           ccGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCaa                                                                                                                             >  1:62184/1‑122 (MQ=255)
                                           ccGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCaa                                                                                                                             <  2:62184/122‑1 (MQ=255)
                                             gcgcgTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCTGTGGCTGGCGGGTGGTCAAc                                                                                                <  2:345393/149‑1 (MQ=255)
                                                  ttGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATtg                                                                                           <  5:50682/149‑1 (MQ=255)
                                                  ttGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATtg                                                                                           <  5:42761/149‑1 (MQ=255)
                                                       gcAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCg                                                                                      >  3:347453/1‑149 (MQ=255)
                                                             gCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCg                                                                                >  1:283860/1‑149 (MQ=255)
                                                                 tgtCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTc                                                                             <  2:283860/148‑1 (MQ=255)
                                                                   tCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGa                                                                           >  5:149487/1‑148 (MQ=255)
                                                                     gATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCTGGTGGTCAACCATTGCTGCGGGGCCGGTCGACaa                                                                        <  5:258363/149‑1 (MQ=255)
                                                                         cGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCt                                                                    >  1:203501/1‑149 (MQ=255)
                                                                          gTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTg                                                                   >  2:196832/1‑149 (MQ=255)
                                                                                ggTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAg                                                                                                                            <  1:128189/86‑1 (MQ=255)
                                                                                ggTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAg                                                                                                                            >  2:128189/1‑86 (MQ=255)
                                                                                    cTTACAGTTGCGGTGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGt                                                          <  2:213989/148‑1 (MQ=255)
                                                                                     ttACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGc                                                        >  2:32048/1‑149 (MQ=255)
                                                                                     ttACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGAGCAAGTGc                                                        >  4:96120/1‑149 (MQ=255)
                                                                                       aCAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCGCGTCTTAGCTTCGGCGAGGCCGAAGAAGCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCg                                                       <  5:281497/148‑1 (MQ=255)
                                                                                       aCAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGc                                                        <  4:21357/147‑1 (MQ=255)
                                                                                       aCAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGc                                                        <  3:27411/147‑1 (MQ=255)
                                                                                       aCAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCg                                                       <  6:287848/148‑1 (MQ=255)
                                                                                        cAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGt                                                      <  2:203501/148‑1 (MQ=255)
                                                                                           ttGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACa                                                  <  2:121316/149‑1 (MQ=255)
                                                                                           ttGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACACAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGc                                                        >  6:44240/1‑143 (MQ=255)
                                                                                           ttGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACACAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGc                                                        <  5:44240/143‑1 (MQ=255)
                                                                                              cGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTg                                               >  6:316925/1‑149 (MQ=255)
                                                                                               gggggCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTgg                                              <  6:149487/149‑1 (MQ=255)
                                                                                                 gggCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTc                                                                                                   >  4:24999/1‑94 (MQ=255)
                                                                                                 gggCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTc                                                                                                   <  3:24999/94‑1 (MQ=255)
                                                                                                 gggCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGGCCACCATTGCGGCGGGGCCGGGCGCCACCCCGGGGGGCCAGGGCGTAAcccagcc                                             >  7:52330/1‑142 (MQ=255)
                                                                                                     aGCCGCTGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATtgctgc                                                                                       <  3:33214/102‑1 (MQ=255)
                                                                                                     aGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATtgctgc                                                                                       >  4:33214/1‑102 (MQ=255)
                                                                                                        cgcgGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATc                                     <  5:80611/149‑1 (MQ=255)
                                                                                                         gcgGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATcc                                    >  3:30669/1‑149 (MQ=255)
                                                                                                            gCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACg                                 <  5:316925/149‑1 (MQ=255)
                                                                                                                gTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGt                               <  5:61793/147‑1 (MQ=255)
                                                                                                                gTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTc                              >  3:246610/1‑148 (MQ=255)
                                                                                                                 tACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCaa                            <  4:30669/149‑1 (MQ=255)
                                                                                                                   ccGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAg                           <  1:196832/148‑1 (MQ=255)
                                                                                                                   ccGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAg                           <  7:213424/148‑1 (MQ=255)
                                                                                                                          cccGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCAc                    >  5:195876/1‑148 (MQ=255)
                                                                                                                             gTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGCCGCACCATg                >  7:242514/1‑149 (MQ=255)
                                                                                                                                     ttCGGCCTGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACgg                                >  2:14645/1‑125 (MQ=255)
                                                                                                                                     ttCGGCCTGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACgg                                <  1:14645/125‑1 (MQ=255)
                                                                                                                                          ccAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGCCCCCg   <  4:246610/149‑1 (MQ=255)
                                                                                                                                            aGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGCCCCCGc  >  2:52723/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
CACCGTTCCATGTCCCCGTTCGGGCTGGCGGGATGACGCCTTCCCGCGCGTTGGCGCAGTGGCTGTGTCGATCCGTTGTCGGTGCTTACAGTTGCGGGGGCAGCCGCGGCTTGTACCGCGTTCCCGTCTTAGCTTCGGCCAGGCCGAAGAACCTCGAAGTGGCAAGGCTACGCAGTGGCTGGCGGGTGGTCAACCATTGCTGCGGGGCCGGTCGACAACGCTGGGCGCAAGTGCGTTACAGTGGCGCCTTATCCACGGTCAAGGACGCACCATGCCCCTCGCCCCCGC  >  NC_002947/3981837‑3982124

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