Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,768,939:1 +C 5.8% noncoding (35/209 nt) PP_mr34 ← Cobalamin regulatory region

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,768,9391.C 5.8% 280.9 / 10.9 69noncoding (35/209 nt)PP_mr34Cobalamin regulatory region
Reads supporting (aligned to +/- strand):  ref base . (35/30);  new base C (2/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

GGCTTGCGAGGGGCATGAAACCTTCACCTTCCCATGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTGCGGCGCCGAGACCAAATGTGGCATGAACATCA  >  NC_002947/2768798‑2769084
                                                                                                                                              |                                                                                                                                                  
ggCTTGCGAGGGGCATGAAACCTTCACCTTCCCATGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCgcggcg                                                                                                                                             >  7:52745/1‑149 (MQ=255)
         ggggCATGAACCCATCCCGTTCCCATGCTTGATCACAGTGGAATCTCGAACGGTTCGCTCGCTTACCGTTGCGGGTGCAGCACCGGACTTCTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCCGCCGCGCGCAggg                                                                                                                                    <  7:139306/148‑1 (MQ=255)
          gggCATGAAACCTTCACCTTCCCATGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCa                                                                                                                                  >  6:371723/1‑149 (MQ=255)
           ggCATGAAACCTTCACCTTCCCATGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCCGCCGCGCGCAgggtga                                                                                                                                 <  6:88005/149‑4 (MQ=255)
            gCATGAAACCTTCACCTTCCCATGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCAcc                                                                                                                                >  2:255074/1‑149 (MQ=255)
             cATGAAACCTTCACCTTCCCATGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCt                                                                                                                               <  7:203490/149‑1 (MQ=255)
              aTGAAACCTTCACCTTCCCATGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCt                                                                                                                               >  1:348699/1‑148 (MQ=255)
                    ccttcaccttcCCATGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTg                                                                                                                        <  7:113013/149‑1 (MQ=255)
                     cttcaccttcCCATGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTgg                                                                                                                       >  3:362168/1‑149 (MQ=255)
                          ccttcCCATGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGgcgcg                                                                                                                   >  8:180859/1‑148 (MQ=255)
                           cttcCCATGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCa                                                                                                                 <  2:159541/149‑1 (MQ=255)
                           cttcCCATGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCa                                                                                                                 >  4:218328/1‑149 (MQ=255)
                               ccATGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCt                                                                                                              >  2:353761/1‑148 (MQ=255)
                               ccATGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCt                                                                                                              <  3:218328/148‑1 (MQ=255)
                                 aTGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGAc                                                                                                           >  2:163490/1‑149 (MQ=255)
                                 aTGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGAc                                                                                                           >  5:362457/1‑149 (MQ=255)
                                 aTGCTTGAGCACAGTGGTAGCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGAc                                                                                                           <  7:53218/149‑1 (MQ=255)
                                  tGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGAcc                                                                                                          <  2:348699/149‑1 (MQ=255)
                                     ttGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACt                                                                                                       >  8:243496/1‑149 (MQ=255)
                                          cacaGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAAt                                                                                                  >  6:350849/1‑149 (MQ=255)
                                             aGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCATGGGCAACACGACACACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCt                                                                                               >  7:91399/1‑149 (MQ=255)
                                              gTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTg                                                                                              <  2:365789/149‑1 (MQ=255)
                                                   aTCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAggg                                                                                                                                    >  7:142543/1‑106 (MQ=255)
                                                   aTCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAggg                                                                                                                                    <  8:142543/106‑1 (MQ=255)
                                                    tctcGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGc                                                                                          >  2:97961/1‑147 (MQ=255)
                                                    tctcGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGc                                                                                          <  1:97961/147‑1 (MQ=255)
                                                    tctcGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGt                                                                                        >  7:244350/1‑149 (MQ=255)
                                                    tctcGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGt                                                                                        >  1:181643/1‑149 (MQ=255)
                                                       cGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCa                                                                                      >  6:36145/1‑148 (MQ=255)
                                                        gAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTc                                                                                       <  7:180859/146‑1 (MQ=255)
                                                             gTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTg                                                                               >  4:109845/1‑149 (MQ=255)
                                                                  ctcgctTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGtgtg                                                                          <  5:36145/149‑1 (MQ=255)
                                                                      ctTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCtgt                                                                      <  1:255074/149‑1 (MQ=255)
                                                                               gCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTg                                                              <  8:52745/148‑1 (MQ=255)
                                                                                  ggggCAGCACCGGACTTGTCGTGGGCAACCCGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCa                                                                                                                 >  8:18725/1‑94 (MQ=255)
                                                                                  ggggCAGCACCGGACTTGTCGTGGGCAACCCGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCa                                                                                                                 <  7:18725/94‑1 (MQ=255)
                                                                                    ggCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCa                                                                                                                 <  7:181023/92‑1 (MQ=255)
                                                                                    ggCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCa                                                                                                                 >  8:181023/1‑92 (MQ=255)
                                                                                       agcaCCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGaaaaaa                                                      >  6:354256/1‑148 (MQ=255)
                                                                                       agcaCCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAaac                                                     <  2:181643/149‑1 (MQ=255)
                                                                                       agcaCCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAaac                                                     <  5:350849/149‑1 (MQ=255)
                                                                                         caCCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAaacaa                                                   <  6:362457/149‑1 (MQ=255)
                                                                                          aCCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAaacaa                                                   >  1:247450/1‑148 (MQ=255)
                                                                                          aCCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAATGTTGGTGTGCTGTGAATGCTGAGAAAAaacaa                                                   >  1:205066/1‑147 (MQ=255)
                                                                                                 ttGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTtgt                                           >  7:68483/1‑149 (MQ=255)
                                                                                                   gTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGgtgt                                           <  4:362168/147‑1 (MQ=255)
                                                                                                    tCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGc                                         <  8:244350/148‑1 (MQ=255)
                                                                                                    tCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGc                                         <  1:353761/148‑1 (MQ=255)
                                                                                                    tCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCg                                        >  5:379265/1‑149 (MQ=255)
                                                                                                         ggCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGttttt                                   >  3:218800/1‑149 (MQ=255)
                                                                                                         ggCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGAGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGttttt                                   >  6:223744/1‑149 (MQ=255)
                                                                                                           cAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCTTCTGACCACTGAAATGATGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTgc                                 <  7:243496/149‑1 (MQ=255)
                                                                                                            aaCACGACGCACCGGTTTCCCGTTTCACCCTGAGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTgcg                                <  5:371723/149‑1 (MQ=255)
                                                                                                              cacGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTgcggc                              <  4:218800/149‑1 (MQ=255)
                                                                                                              cacGACGCACCGGTTTCCCGTTTCACCCTGCGCCCGGCCGCGCAGGGCACATGAAACTGGCGCGCATCTGACAACTTAAATGATGGGGCGTCAAATGTTTGTATGCTCTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTgcggc                              >  6:76565/1‑149 (MQ=255)
                                                                                                               acgacgCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTgcggcg                             >  5:313799/1‑149 (MQ=255)
                                                                                                                 gacgCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTGCGGCGcc                           >  6:315736/1‑149 (MQ=255)
                                                                                                                      aCCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTGCGGCGCCGAGAc                      <  3:109845/149‑1 (MQ=255)
                                                                                                                        cGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTGCGGCGCCGAGAcc                     >  1:161062/1‑148 (MQ=255)
                                                                                                                         ggTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTGCGGCGCCGAGACCa                    >  7:71904/1‑148 (MQ=255)
                                                                                                                          gTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTGCGGCGCCGAGACCaaa                  <  1:163490/149‑1 (MQ=255)
                                                                                                                               ccGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTGCGGCGCCGAGACCAAATGTgg             <  8:91399/149‑1 (MQ=255)
                                                                                                                                 gTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTGCGGCGCCGAGACCAAATGTGGCa           <  6:379265/149‑1 (MQ=255)
                                                                                                                                   ggtgcccTGCGCCGCCGCGCGCAGGGTACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGTTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTGCGGCGcc                           <  7:99397/127‑1 (MQ=255)
                                                                                                                                   ggtgcccTGCGCCGCCGCGCGCAGGGTACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGTTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTGCGGCGcc                           >  8:99397/5‑131 (MQ=255)
                                                                                                                                      gcccTGCGCCGCCGCGCGCAGGGTACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCt                                                               <  4:113789/91‑1 (MQ=38)
                                                                                                                                      gcccTGCGCCGCCGCGCGCAGGGTACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGGCAAATGTTGGTGTGCTGTGAATGCt                                                               >  3:113789/2‑92 (MQ=37)
                                                                                                                                      aCCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTGCGGCGCCGAGACCAAATGTGGCATGaa       >  1:58016/1‑148 (MQ=255)
                                                                                                                                      aCCCTGCGAGCGGCCGCGCAGGGAACCTGAAACTGGCGCGCATATGACCACTGAAATGCTGGGGCGTAAAATGTAGGTGTGCTCTGAATCCTGAGAAAAAACAACAGGTTGAGAGAGTTTGCGCCTCCGAGACCAAATGTGGCATGaa       >  4:77218/1‑148 (MQ=255)
                                                                                                                                            cgcgcgGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTGCGGCGCCGAGACCAAATGTGGCATGAACATCa  <  2:247450/147‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GGCTTGCGAGGGGCATGAAACCTTCACCTTCCCATGCTTGAGCACAGTGGTATCTCGAAGGGTTCGCTCGCTTACCGTTGCGGGGGCAGCACCGGACTTGTCGTGGGCAACACGACGCACCGGTTTCCCGTTTCACCCTGCGCGCGGCGGCGCAGGGCACCTGAAACTGGCGCGCATCTGACCACTGAAATGCTGGGGCGTCAAATGTTGGTGTGCTGTGAATGCTGAGAAAAAACAACAGGTTGTGCGTTTTTGCGGCGCCGAGACCAAATGTGGCATGAACATCA  >  NC_002947/2768798‑2769084

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: