Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 819,241:1 +C 5.8% coding (395/627 nt) PP_0705 → DNA‑3‑methyladenine glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947819,2411.C 5.8% 277.3 / 8.6 69coding (395/627 nt)PP_0705DNA‑3‑methyladenine glycosylase
Reads supporting (aligned to +/- strand):  ref base . (29/36);  new base C (2/2);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

GTTGCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGA  >  NC_002947/819098‑819384
                                                                                                                                                |                                                                                                                                               
gTTGCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATc                                                                                                                                             >  4:309487/1‑148 (MQ=255)
    cTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAAc                                                                                                                                         <  3:65686/148‑1 (MQ=255)
       gCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGc                                                                                                                                       <  6:166471/147‑1 (MQ=255)
       gCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCCTTAACCATGGCTGACGAGGCGTTGATCGAACGCCt                                                                                                                                     >  4:321841/1‑149 (MQ=255)
        cGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATCAACGCCtc                                                                                                                                    <  5:297254/149‑2 (MQ=255)
          gTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTCGATCAACGCCtcgt                                                                                                                                  <  1:27572/149‑1 (MQ=255)
              ccccGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGgcg                                                                                                                              >  3:23569/1‑149 (MQ=255)
               cccGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGgcgc                                                                                                                             >  7:40005/1‑149 (MQ=255)
               cccGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGgcgc                                                                                                                             >  7:174762/1‑149 (MQ=255)
               cccGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGgcgc                                                                                                                             <  5:380693/149‑1 (MQ=255)
                  ggAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGc                                                                                                                          <  5:377834/149‑1 (MQ=255)
                      aCCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTgg                                                                                                                      <  1:303320/149‑1 (MQ=255)
                        ccTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTgggg                                                                                                                    >  6:8489/1‑149 (MQ=255)
                         cTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGt                                                                                                                        >  8:7626/1‑144 (MQ=255)
                         cTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGt                                                                                                                        <  7:7626/144‑1 (MQ=255)
                               gCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcg                                                                                                             >  8:91782/1‑149 (MQ=255)
                                 cTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggt                                                                                                           <  7:233368/149‑1 (MQ=255)
                                     ggCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggac                                                                                                       <  3:309487/149‑1 (MQ=255)
                                        ttttCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTg                                                                                                                       >  5:392330/1‑130 (MQ=255)
                                        ttttCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTg                                                                                                                       <  6:392330/130‑1 (MQ=255)
                                          ttCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtg                                                                                                   <  3:321841/148‑1 (MQ=255)
                                          ttCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtg                                                                                                   <  4:377768/148‑1 (MQ=255)
                                          ttCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtg                                                                                                   <  6:248748/148‑1 (MQ=255)
                                            cGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtgga                                                                                                 <  6:113576/148‑1 (MQ=255)
                                                  gCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGAtgctg                                                                                          <  2:201905/149‑1 (MQ=255)
                                                   cAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGAtgctgc                                                                                         >  6:361256/1‑149 (MQ=255)
                                                   cAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGAtgctgc                                                                                         <  2:307178/149‑1 (MQ=255)
                                                       cTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGAt                                                                                     <  4:23569/149‑1 (MQ=255)
                                                           tGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCtat                                                                                 <  7:91782/149‑1 (MQ=255)
                                                                 tCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTg                                                                           <  5:8489/149‑1 (MQ=255)
                                                                 tCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTg                                                                           >  2:321994/1‑149 (MQ=255)
                                                                   cAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtg                                                                                                   <  3:324596/123‑1 (MQ=255)
                                                                   cAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtg                                                                                                   >  4:324596/1‑123 (MQ=255)
                                                                        cATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGaa                                                                                                                                          >  8:211273/1‑79 (MQ=255)
                                                                        cATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGtt                                                                    <  1:322250/149‑1 (MQ=255)
                                                                        cATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGtt                                                                    <  7:55955/149‑1 (MQ=255)
                                                                             cTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGGGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGa                                                                >  8:86277/1‑148 (MQ=255)
                                                                              tCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACa                                                              >  7:57944/1‑149 (MQ=255)
                                                                                 aGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGGGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTc                                                           <  7:97449/149‑1 (MQ=255)
                                                                                   ccACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTg                                                         <  7:306/149‑1 (MQ=255)
                                                                                    cACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGc                                                        <  8:174762/149‑1 (MQ=255)
                                                                                    cACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGc                                                        >  7:24231/1‑149 (MQ=255)
                                                                                     aCCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGcc                                                       >  3:384665/1‑149 (MQ=255)
                                                                                       ccTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGACCGCCTGGTGGCGCTTCGTGGCGTAGGGCGGGGGACGGTGGAGATGCTGCTGATCTATAGCCTGCGGCGGTATGACATTCTGCAgg                                                     >  8:28807/1‑149 (MQ=255)
                                                                                       ccTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCTCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCgg                                                     <  8:268877/149‑1 (MQ=255)
                                                                                        cTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCCTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGt                                                    <  1:321994/149‑1 (MQ=255)
                                                                                         tGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTg                                                   >  7:365305/1‑149 (MQ=255)
                                                                                          ggAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTg                                                   <  8:365305/148‑1 (MQ=255)
                                                                                               gCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGAt                                             >  3:316484/1‑149 (MQ=255)
                                                                                               gCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGGGGACGGTGGAGATGCTGCTGAGCTATAGCCTGGGGCGTACTGACATGCTGCCGGTGGACCAt                                             >  3:350571/1‑149 (MQ=255)
                                                                                                  tGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGTGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGAtttt                                          >  4:314678/1‑149 (MQ=255)
                                                                                                     tACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGCCGATTTTGGc                                       >  4:381750/1‑149 (MQ=255)
                                                                                                     tACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGc                                       <  5:361256/149‑1 (MQ=255)
                                                                                                     tACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGc                                       >  2:374181/1‑149 (MQ=255)
                                                                                                        cGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTg                                    >  6:27842/1‑149 (MQ=255)
                                                                                                           cgcgGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTgcgc                                 <  8:57944/149‑1 (MQ=255)
                                                                                                            gcgGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGCACGCCTGGTGGCGCTGCGTGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTgcgc                                 >  4:285724/1‑147 (MQ=255)
                                                                                                             cgGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGa                               <  4:344877/149‑1 (MQ=255)
                                                                                                               ggAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAgg                             <  5:27842/149‑1 (MQ=255)
                                                                                                                                tGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAg            <  8:40005/149‑1 (MQ=255)
                                                                                                                                  gcggACGAGGGGTTGATCGTACGGCTGGTGGCGCCGCGGTGGGTACGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAggg          <  4:350571/146‑1 (MQ=255)
                                                                                                                                  gCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAggg          <  4:316484/149‑1 (MQ=255)
                                                                                                                                      accaGGCGTTCGATCAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATAGGCGGCTCAAGGGGATg      >  8:351753/4‑149 (MQ=255)
                                                                                                                                      aCGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGc        <  1:374181/147‑1 (MQ=255)
                                                                                                                                      aCGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTg      >  5:127678/1‑149 (MQ=255)
                                                                                                                                      aCGAGGCGTTGATCGAACGCCTGGGGGCGCTGCGGGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATCGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTg      <  7:28807/149‑1 (MQ=255)
                                                                                                                                       ccaGGCGTTCGATCAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTg      >  6:279560/3‑148 (MQ=255)
                                                                                                                                          ggCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGgaga  >  7:160102/1‑149 (MQ=255)
                                                                                                                                          ggCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGgaga  >  5:117327/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GTTGCTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGA  >  NC_002947/819098‑819384

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: