Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,594,487 G→C 9.8% R312R (CGC→CGG kgtP ← alpha‑ketoglutarate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,594,4870GC 9.8% 178.8 / 10.3 61R312R (CGC→CGGkgtPalpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base G (26/29);  new base C (4/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-02

GTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTGGTGTAGGTGTA  >  NC_002947/1594340‑1594633
                                                                                                                                                   |                                                                                                                                                  
gTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGAcggc                                                                                                                                                   <  6:285983/149‑1 (MQ=255)
    aTCGAGGTATAGCCGCTTACGATGATCAGTGGGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGcc                                                                                                                                               >  1:232774/1‑149 (MQ=255)
     tCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGcc                                                                                                                                               >  8:211624/1‑148 (MQ=255)
       gAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGa                                                                                                                                             <  2:272629/148‑1 (MQ=255)
        aGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGCCGCCGTCCGATc                                                                                                                                           <  3:135054/149‑1 (MQ=255)
        aGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATc                                                                                                                                           >  4:73539/1‑149 (MQ=255)
        aGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATc                                                                                                                                           >  7:281362/1‑149 (MQ=255)
            ataGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGt                                                                                                                                       >  3:341476/1‑149 (MQ=255)
             taGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTc                                                                                                                                      <  1:357756/149‑1 (MQ=255)
               gCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCg                                                                                                                                     >  4:198359/1‑148 (MQ=255)
                ccGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGa                                                                                                                                   >  2:146094/1‑149 (MQ=255)
                ccGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGa                                                                                                                                   >  7:243920/1‑149 (MQ=255)
                      aCGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGac                                                                                                                              >  4:344735/1‑148 (MQ=255)
                      aCGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGac                                                                                                                              <  3:73539/148‑1 (MQ=255)
                      aCGATGATCAGTGCGGCAATGATCAGGAAGAACGCACCCCACCAGGTCTTGATGGTGTGCAAGGTGCTGAGGATGGGGAAGGTGAACAGTGAGCCCAGCACACCGAAGGAGATCAGGATCGGACGGATGCCCATCTTGTCGGACAGac                                                                                                                              >  4:241467/1‑148 (MQ=255)
                                         tGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAgg                                                                                                          <  2:232774/149‑1 (MQ=255)
                                           aTCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGca                                                                                                        <  5:298572/149‑1 (MQ=255)
                                            tCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGcac                                                                                                       <  1:146094/149‑1 (MQ=255)
                                            tCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGcac                                                                                                       <  1:138929/149‑1 (MQ=255)
                                             cAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGGCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGcaca                                                                                                      >  7:178649/1‑149 (MQ=255)
                                              aGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACAt                                                                                                     >  4:96560/1‑149 (MQ=255)
                                              aGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACAt                                                                                                     <  4:341476/149‑1 (MQ=255)
                                                gaagaaCGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATg                                                                                                    >  7:236206/1‑148 (MQ=255)
                                                 aagaaCGCACCCCACCAGGTCTGGATGGTGCGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATgaa                                                                                                  >  3:53796/1‑149 (MQ=255)
                                                           cccACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGc                                                                                        >  6:203569/1‑149 (MQ=255)
                                                           cccACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGc                                                                                        >  4:43927/1‑149 (MQ=255)
                                                             caccaGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGt                                                                                      >  6:120539/1‑149 (MQ=255)
                                                             caccaGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGt                                                                                      <  8:29689/149‑1 (MQ=255)
                                                               ccaGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTgg                                                                                    <  4:361398/149‑1 (MQ=255)
                                                                 aGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTTCAGGCACATGAACAGGAACAGCGTggcg                                                                                  <  8:275947/149‑1 (MQ=255)
                                                                    tCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGcc                                                                               <  8:281362/149‑1 (MQ=255)
                                                                                  aGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGa                                                                 <  7:211624/149‑1 (MQ=255)
                                                                                    ggTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGt                                                               <  3:96560/149‑1 (MQ=255)
                                                                                        cTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTggcg                                                                                  >  1:81270/1‑126 (MQ=255)
                                                                                        cTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTggcg                                                                                  <  2:81270/126‑1 (MQ=255)
                                                                                            ggATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATg                                                       <  8:243920/149‑1 (MQ=255)
                                                                                             gATGGGGACGGTGTACAGGGTGCCGTGCACACCGAAGGCGAGCAGGATCGTCCGGCGGCCGTTCTTGTGGGAGTGACCACCGATCATCGGATGCAGGCACATGAACAGGAACCGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGc                                                      <  3:241467/149‑1 (MQ=255)
                                                                                                 gggACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCg                                                                              <  4:353609/121‑1 (MQ=255)
                                                                                                 gggACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCg                                                                              >  3:353609/1‑121 (MQ=255)
                                                                                                 gggACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCAt                                                  >  1:204084/1‑149 (MQ=255)
                                                                                                      ggTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCa                                             >  8:252807/1‑149 (MQ=255)
                                                                                                        tGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACg                                           <  3:344735/149‑1 (MQ=255)
                                                                                                         gAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTggcggc                                                                                <  4:269759/111‑1 (MQ=255)
                                                                                                           aCAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGgtg                                        <  3:43927/149‑1 (MQ=255)
                                                                                                            gatggTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGTCAGACGACCGATCATCGGCTGCAGGCACATGAACATGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGgtgt                                       <  5:120539/146‑1 (MQ=255)
                                                                                                                   ccGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAg                                <  4:53796/149‑1 (MQ=255)
                                                                                                                   ccGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAg                                <  7:221002/149‑1 (MQ=255)
                                                                                                                     gAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAgg                               >  5:296835/1‑148 (MQ=255)
                                                                                                                       gCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTAt                            <  2:204084/149‑1 (MQ=255)
                                                                                                                           acCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCt                        >  8:153107/1‑149 (MQ=255)
                                                                                                                                 ggCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGt                  <  2:311990/149‑1 (MQ=255)
                                                                                                                                    gATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCGTGCCCACGGTGTTGACCAGGTATTTCTGCATGTAgg               <  5:203569/149‑1 (MQ=255)
                                                                                                                                     gacaagATCGGCCGCCGTCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAggt              >  3:346111/6‑149 (MQ=255)
                                                                                                                                      tCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAggtg             <  6:296835/149‑1 (MQ=255)
                                                                                                                                      acaagATCGGCCGCCGTCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCAt                                                  >  7:264441/5‑112 (MQ=255)
                                                                                                                                      acaagATCGGCCGCCGTCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCAt                                                  <  8:264441/108‑1 (MQ=255)
                                                                                                                                       caagATCGGCCGCCGTCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGccc                                              >  6:244970/4‑115 (MQ=255)
                                                                                                                                          gATCGGCCGCCGTCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGgtggtg          >  2:380392/1‑148 (MQ=255)
                                                                                                                                          gATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGgtggtg          >  5:353308/1‑148 (MQ=255)
                                                                                                                                            tCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTggtgtag       <  7:252807/149‑1 (MQ=255)
                                                                                                                                              ggACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTggtgtaggt     >  7:312691/1‑149 (MQ=255)
                                                                                                                                               gACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTggtgtaggtg    >  3:276109/1‑149 (MQ=255)
                                                                                                                                                  ggcggcCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAg                <  6:310415/134‑1 (MQ=255)
                                                                                                                                                  ggcggcCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAg                >  5:310415/1‑134 (MQ=255)
                                                                                                                                                  ggcggaCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTggtgtaggtgta  <  8:178649/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTGGTGTAGGTGTA  >  NC_002947/1594340‑1594633

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: