Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 819,246 Δ1 bp 5.7% coding (400/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,2460G. 5.7% 228.3 / 6.8 70coding (400/627 nt)PP_0705DNA‑3‑methyladenine glycosylase
Reads supporting (aligned to +/- strand):  ref base G (29/37);  new base . (2/2);  total (31/39)
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.12e-01

CTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACGGCGCC  >  NC_002947/819102‑819391
                                                                                                                                                 |                                                                                                                                                 
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                                                                                                                                    <  5:380693/149‑1 (MQ=255)
           cccGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGgcgc                                                                                                                                    >  7:40005/1‑149 (MQ=255)
           cccGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGgcgc                                                                                                                                    >  7:174762/1‑149 (MQ=255)
              ggAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGc                                                                                                                                 <  5:377834/149‑1 (MQ=255)
                  aCCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTgg                                                                                                                             <  1:303320/149‑1 (MQ=255)
                    ccTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTgggg                                                                                                                           >  6:8489/1‑149 (MQ=255)
                     cTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGt                                                                                                                               <  7:7626/144‑1 (MQ=255)
                     cTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGt                                                                                                                               >  8:7626/1‑144 (MQ=255)
                           gCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcg                                                                                                                    >  8:91782/1‑149 (MQ=255)
                             cTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggt                                                                                                                  <  7:233368/149‑1 (MQ=255)
                                 ggCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggac                                                                                                              <  3:309487/149‑1 (MQ=255)
                                    ttttCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTg                                                                                                                              <  6:392330/130‑1 (MQ=255)
                                    ttttCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTg                                                                                                                              >  5:392330/1‑130 (MQ=255)
                                      ttCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtg                                                                                                          <  3:321841/148‑1 (MQ=255)
                                      ttCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtg                                                                                                          <  6:248748/148‑1 (MQ=255)
                                      ttCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtg                                                                                                          <  4:377768/148‑1 (MQ=255)
                                        cGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGcggtggacggtgga                                                                                                        <  6:113576/148‑1 (MQ=255)
                                              gCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGAtgctg                                                                                                 <  2:201905/149‑1 (MQ=255)
                                               cAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGAtgctgc                                                                                                <  2:307178/149‑1 (MQ=255)
                                               cAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGAtgctgc                                                                                                >  6:361256/1‑149 (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                                                               >  7:24231/1‑149 (MQ=255)
                                                                                cACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGc                                                               <  8:174762/149‑1 (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                                              >  2:374181/1‑149 (MQ=255)
                                                                                                 tACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGc                                              <  5:361256/149‑1 (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)
                                                                                                                                              aTCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAgg                                    <  6:234153/115‑1 (MQ=255)
                                                                                                                                              aTCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACGGCGcc  >  8:94255/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CTGGCGGTATCCCCGGAAACCCTGCGTGCCTGCGGCTTTTCGGCGAGCAAGCTGGTGACGGTCCAGGGCATTGCTCAAGCCACCCTGGAGGGCCTGGTACCGACGCGGGAGCAGGCGTTAACCATGGCTGACGAGGCGTTGATCGAACGCCTGGTGGCGCTGCGTGGGGTAGGGCGGTGGACGGTGGAGATGCTGCTGATCTATAGCCTGGGGCGTTCTGACATTCTGCCGGTGGACGATTTTGGCGTGCGCGAGGGGTATCGGCGGCTCAAGGGGCTGGAGACGGCGCC  >  NC_002947/819102‑819391

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: