Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,478 G→C 19.6% intergenic (+290/‑641) RCAP_rcc01306 → / → RCAP_rcc01307 Cas2 family CRISPR‑associated protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,395,4780G→C19.6% 118.4 ‑1.1 ‑0.0240intergenic (+290/‑641)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (18/29):  ref base (90/103):  total (108/132)
Fisher's exact test strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.30e-01

Only 100 of 250 total aligned reads displayed.

CGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCC  >  NC_014034/1395378‑1395578
                                                                                                    |                                                                                                    
cgGGGGCGTGGATCGAAACTGCCCCCGCTGGCCCGGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTgg                                                                                                      <  1:1333816/101‑1 (MQ=255)
     gCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTc                                                                                                 >  1:261755/1‑101 (MQ=255)
      cGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCg                                                                                                >  1:2938501/1‑101 (MQ=255)
      cGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCg                                                                                                >  1:2438469/1‑101 (MQ=255)
       gTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGa                                                                                               >  1:3875017/1‑101 (MQ=255)
        tGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGAc                                                                                              >  2:2726699/1‑101 (MQ=255)
         ggATCGAAACTGCACCCGCTGGCCCTGCCTGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACa                                                                                             >  1:2624157/1‑101 (MQ=255)
         ggATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACa                                                                                             >  1:988216/1‑101 (MQ=255)
           aTCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAgg                                                                                           >  2:117086/1‑101 (MQ=255)
              gAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTg                                                                                        <  1:848362/101‑1 (MQ=255)
                 aCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTa                                                                                     <  1:3466905/101‑1 (MQ=255)
                   tGCCCCCGCTGCCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAc                                                                                   <  1:4159381/101‑1 (MQ=255)
                   tGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAc                                                                                   <  1:3913955/101‑1 (MQ=255)
                     cccccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACgt                                                                                 <  2:1625999/99‑1 (MQ=255)
                     cACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACgt                                                                                 >  1:1655428/1‑101 (MQ=255)
                      aCCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGGAACgtg                                                                                >  2:3057989/1‑101 (MQ=255)
                       cccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACgtgt                                                                               <  1:516525/101‑1 (MQ=255)
                        ccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGTGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTc                                                                              <  2:335861/101‑1 (MQ=255)
                         cGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCg                                                                             >  2:329753/1‑101 (MQ=255)
                         cGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCg                                                                             <  1:1212377/101‑1 (MQ=255)
                          gCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGGAACGTGTCGc                                                                            >  2:3001589/1‑101 (MQ=255)
                           cTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCt                                                                           <  2:1713654/101‑1 (MQ=255)
                             ggCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcc                                                                         >  1:789910/1‑101 (MQ=255)
                               cccTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccc                                                                       >  2:4022521/1‑101 (MQ=255)
                                ccTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTccccc                                                                      <  1:100034/101‑1 (MQ=255)
                                 cTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccccc                                                                     <  2:1587702/101‑1 (MQ=255)
                                 cTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccccc                                                                     <  2:52103/101‑1 (MQ=255)
                                         gatcggtgatcAGCGGTCGCTCCCCCCGCGTGGGCGTGGATCGAAACATCGCCGCTTTGGCAGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGTGc                                                             <  2:3481269/101‑1 (MQ=255)
                                         gatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGc                                                             >  2:566937/1‑101 (MQ=255)
                                         gatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGc                                                             <  2:971995/101‑1 (MQ=255)
                                          atcggtgatcAGCGGTCGCTCCCCCCGCGGGGGTTTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCg                                                            <  1:3630275/101‑1 (MQ=255)
                                           tcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGt                                                           <  2:4177425/101‑1 (MQ=255)
                                            cggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTg                                                          >  2:3937025/1‑101 (MQ=255)
                                            cggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTg                                                          <  2:1526036/101‑1 (MQ=255)
                                            cggtgatcAGCGGACGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTg                                                          <  2:3089563/101‑1 (MQ=255)
                                             ggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTgg                                                         <  1:3371570/101‑1 (MQ=255)
                                             ggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTgg                                                         <  2:3912547/101‑1 (MQ=255)
                                             ggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTgg                                                         >  1:1832167/1‑101 (MQ=255)
                                              gtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGa                                                        >  1:107331/1‑101 (MQ=255)
                                              gtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGa                                                        <  2:4127914/101‑1 (MQ=255)
                                              gtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGa                                                        >  2:3691073/1‑101 (MQ=255)
                                              gctggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGGGGAt                                                       >  1:2216980/10‑101 (MQ=2)
                                               tgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGAt                                                       <  1:7210/101‑1 (MQ=255)
                                               ctggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                      >  1:1539092/9‑101 (MQ=0)
                                                gatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATc                                                      <  1:3176755/101‑1 (MQ=255)
                                                 ggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                    >  2:1743447/7‑101 (MQ=2)
                                                 atcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                     >  2:366066/1‑101 (MQ=255)
                                                  gacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                   >  2:3206203/6‑101 (MQ=2)
                                                   cAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                   <  1:4228421/101‑1 (MQ=255)
                                                    cgtggTCGCTCCCCCCGCGGGGGGGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                 >  1:398216/4‑101 (MQ=2)
                                                    cgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                 >  2:599069/4‑101 (MQ=2)
                                                    aGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                  <  1:1283247/101‑1 (MQ=255)
                                                     gtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                <  2:2800312/99‑1 (MQ=2)
                                                     gtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                >  1:4108360/3‑101 (MQ=2)
                                                     gCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                 <  2:204807/101‑1 (MQ=255)
                                                     gCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                 <  2:2941449/101‑1 (MQ=255)
                                                      tggTCGCTCCCCCCGCGGGGGCTTGGATCGAAACTACACCGGGTCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACag                                               <  2:4226238/100‑2 (MQ=2)
                                                      tggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACag                                               <  2:1564782/100‑2 (MQ=2)
                                                       ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACagc                                              >  2:640582/1‑99 (MQ=2)
                                                       ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACagc                                              <  1:3625243/101‑3 (MQ=2)
                                                         tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtct                                             <  2:1977338/101‑1 (MQ=255)
                                                         tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtct                                             >  2:2483019/1‑101 (MQ=255)
                                                          cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtctc                                            <  1:2853850/101‑1 (MQ=255)
                                                          cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtctc                                            <  1:5830/101‑1 (MQ=255)
                                                           gCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttc                                          >  1:1979495/1‑98 (MQ=2)
                                                           gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg                                           >  1:1447941/1‑101 (MQ=255)
                                                            cTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGTCGTGGATCGAAACAgctttca                                         <  1:2399969/101‑5 (MQ=2)
                                                            cTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttca                                         <  2:1259552/101‑5 (MQ=2)
                                                             tCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGt                                         >  1:2359869/1‑101 (MQ=255)
                                                             tCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGt                                         >  1:2066742/1‑101 (MQ=255)
                                                             cccccccGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaa                                        <  2:2682562/100‑6 (MQ=0)
                                                              ccccccGCGGGGGTGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTg                                        <  2:1047259/101‑1 (MQ=255)
                                                              ccccccGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaat                                       <  2:2921057/101‑7 (MQ=0)
                                                              ccccccGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaat                                       <  1:2748894/101‑7 (MQ=0)
                                                               cccccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGc                                       <  1:2472404/101‑1 (MQ=255)
                                                               cccccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGc                                       <  1:1897925/101‑1 (MQ=255)
                                                                ccccGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaatgg                                     <  2:1083294/101‑9 (MQ=2)
                                                                 cccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCg                                     <  1:3261659/101‑1 (MQ=255)
                                                                   cgcgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTc                                   >  1:554238/1‑101 (MQ=255)
                                                                     cgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGc                                 <  1:158392/101‑1 (MQ=255)
                                                                          gCGTTGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCTCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTc                            <  2:554174/101‑1 (MQ=255)
                                                                          gCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTc                            <  2:3404148/101‑1 (MQ=255)
                                                                           cGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTTGATCGAAACATCTCGGTGCCGTCGCTCGTcc                           <  1:3754861/101‑1 (MQ=255)
                                                                            gTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCa                          <  2:3829519/101‑1 (MQ=255)
                                                                             tGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAg                         <  2:2219896/101‑1 (MQ=255)
                                                                              ggATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGc                        >  2:1836906/1‑101 (MQ=255)
                                                                                  cGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCCGc                    <  1:3128442/101‑1 (MQ=255)
                                                                                   gAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGcc                   >  2:346787/1‑101 (MQ=255)
                                                                                   atctcATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGcc                   <  2:912609/97‑1 (MQ=255)
                                                                                    aaaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCt                  <  2:2438293/101‑1 (MQ=255)
                                                                                     aaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTg                 <  2:2167831/101‑1 (MQ=255)
                                                                                     aaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTg                 <  1:2257267/101‑1 (MQ=255)
                                                                                        aTCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtc              <  2:3556099/101‑1 (MQ=255)
                                                                                          cgccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgt            >  1:2983715/1‑101 (MQ=255)
                                                                                           gccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgtc           <  2:2066604/101‑1 (MQ=255)
                                                                                               cTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTc       <  1:3725765/101‑1 (MQ=255)
                                                                                                tttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTcc      >  1:3211688/1‑101 (MQ=255)
                                                                                                tttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTcc      >  1:1740061/1‑101 (MQ=255)
                                                                                                   ggCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTccccc   >  2:3521967/1‑101 (MQ=255)
                                                                                                   ggCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTccccc   >  1:4171262/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCC  >  NC_014034/1395378‑1395578

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG