Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,475 T→G 18.3% intergenic (+287/‑644) 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,4750T→G18.3% 244.1 ‑0.4 0.0499intergenic (+287/‑644)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (42/50):  ref base (197/210):  total (239/260)
Fisher's exact test strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 520 total aligned reads displayed.

CCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCC  >  NC_014034/1395375‑1395575
                                                                                                    |                                                                                                    
ccGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCtt                                                                                                      <  2:33943/101‑1 (MQ=255)
 cgcgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCttt                                                                                                     <  1:7862997/101‑1 (MQ=255)
  gcgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTg                                                                                                    <  2:2303248/101‑1 (MQ=255)
    gggggGGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGc                                                                                                  <  1:2327255/101‑1 (MQ=255)
    gggggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGc                                                                                                  >  1:5773559/1‑101 (MQ=255)
    gggggCGTCGATCGAATCTGCGCCCGCTGGGCCTGCCGGTGATCGGTGATCAGCGGTCGGTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGc                                                                                                  <  1:9151947/101‑1 (MQ=38)
      gggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCgg                                                                                                >  1:819161/1‑101 (MQ=255)
         cGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCg                                                                                             >  2:5076557/1‑101 (MQ=255)
         cGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCg                                                                                             <  1:639011/101‑1 (MQ=255)
         cGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCg                                                                                             <  2:7268588/101‑1 (MQ=255)
           gggATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGAc                                                                                           <  1:802163/100‑1 (MQ=255)
           gggATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGAc                                                                                           <  2:6355547/100‑1 (MQ=255)
            ggATCGAAACTGCCCCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACa                                                                                          <  2:5523758/101‑1 (MQ=255)
            ggATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGAcg                                                                                          >  2:7879826/1‑100 (MQ=255)
            ggATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACa                                                                                          <  1:4050882/101‑1 (MQ=255)
             gATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAg                                                                                         <  1:8323856/101‑1 (MQ=255)
               tCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGc                                                                                       >  1:6917677/1‑101 (MQ=255)
                   aaCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGt                                                                                   <  1:3528374/101‑1 (MQ=255)
                   aaCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGt                                                                                   <  1:6777684/101‑1 (MQ=255)
                       gCACCCGCTGGCCCTGCCGGTGATCGGTTATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGTAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACg                                                                               >  1:6525000/1‑101 (MQ=255)
                       gCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACg                                                                               >  1:5161286/1‑101 (MQ=255)
                       gCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACg                                                                               <  1:1221939/101‑1 (MQ=255)
                        cccccGCTGGCCCCTCCGGTGATCGGTTATCCCGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACgt                                                                              <  2:4255005/99‑1 (MQ=25)
                         ccccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACgtg                                                                             <  1:6321692/100‑1 (MQ=255)
                         ccccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACgtg                                                                             <  1:4388911/100‑1 (MQ=255)
                         ccccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACgtg                                                                             <  2:6744693/100‑1 (MQ=255)
                          cccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACgtgt                                                                            <  1:18147/101‑1 (MQ=255)
                          cccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACgtgt                                                                            <  1:2484454/101‑1 (MQ=255)
                           ccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTc                                                                           <  1:1006881/101‑1 (MQ=255)
                           ccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTc                                                                           <  2:5851731/101‑1 (MQ=255)
                           ccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTc                                                                           <  2:129209/101‑1 (MQ=255)
                            cGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCg                                                                          <  2:2844114/101‑1 (MQ=255)
                            cGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCg                                                                          >  2:3652529/1‑101 (MQ=255)
                             gCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTGACGCGTCGc                                                                         >  1:8576059/1‑101 (MQ=255)
                                ggCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcc                                                                      <  1:2769083/101‑1 (MQ=255)
                                 gCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTccc                                                                     >  2:731126/1‑101 (MQ=255)
                                  cccTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACCTGTCCCTcccc                                                                    <  2:400051/101‑1 (MQ=255)
                                    cTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccccc                                                                  <  2:8602628/101‑1 (MQ=255)
                                     tGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCcg                                                                 >  1:1565879/1‑101 (MQ=255)
                                      gCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCcgc                                                                >  2:7631801/1‑101 (MQ=255)
                                        cggtgaacggtgatCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCgg                                                              >  2:7467915/1‑101 (MQ=255)
                                          gtgatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCgggg                                                            <  2:9157280/101‑1 (MQ=255)
                                          gtgatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCgggg                                                            >  1:2862823/1‑101 (MQ=255)
                                            gatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGc                                                          <  2:6917392/101‑1 (MQ=255)
                                            gatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGc                                                          <  1:5099003/101‑1 (MQ=255)
                                              tcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGt                                                        >  2:2182117/1‑101 (MQ=255)
                                                 gtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACTTGTCGCTTCCCCCGCGGGGGCGTggg                                                     >  2:3777962/1‑100 (MQ=255)
                                                 gctggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                    <  1:902569/92‑1 (MQ=2)
                                                  tgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGAt                                                    <  1:7220204/101‑1 (MQ=255)
                                                  ggatcAGCGGTCGCGCCCCGCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCGCCCCCCGCGGGGGCGTGGAt                                                    <  1:7387187/100‑1 (MQ=255)
                                                  ctggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGGGGATc                                                   >  1:2105526/9‑101 (MQ=2)
                                                   gatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATc                                                   >  1:6036381/1‑101 (MQ=255)
                                                      gcgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                               <  2:1596254/97‑1 (MQ=2)
                                                       aGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                               >  1:7978424/1‑101 (MQ=255)
                                                        gtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                             <  1:4940053/99‑1 (MQ=2)
                                                        gtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                             <  1:887557/99‑1 (MQ=2)
                                                        gtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                             >  2:8739997/3‑101 (MQ=2)
                                                         cGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGGAACGTGTCGCTCACCCCGCGGGGGCGTGGATCGAACCa                                             >  2:9076379/1‑101 (MQ=255)
                                                          ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACagc                                           <  1:3293409/101‑3 (MQ=2)
                                                          ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACagc                                           >  2:6557023/1‑99 (MQ=2)
                                                           gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgct                                          >  1:3626288/1‑101 (MQ=2)
                                                           gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgct                                          <  1:4750253/101‑1 (MQ=2)
                                                            gcGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtct                                          <  1:2671250/100‑1 (MQ=255)
                                                             ctcTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgcttt                                        <  1:6195478/99‑3 (MQ=2)
                                                             cgatCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgcttt                                        <  2:3572540/98‑3 (MQ=2)
                                                             cGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgcttt                                        >  2:1667243/1‑99 (MQ=2)
                                                              gCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttc                                       <  1:4496069/101‑4 (MQ=2)
                                                              gCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttc                                       <  2:4645279/101‑4 (MQ=2)
                                                              gCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttc                                       <  2:8813743/101‑4 (MQ=2)
                                                               cTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCgg                                       <  2:3220184/101‑1 (MQ=255)
                                                                tCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGt                                      >  1:4566010/1‑101 (MQ=255)
                                                                 ccccccGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaat                                    >  1:544274/1‑95 (MQ=0)
                                                                 ccccccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTg                                     <  2:1378752/101‑1 (MQ=255)
                                                                 ccccccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTg                                     >  1:863547/1‑101 (MQ=255)
                                                                        cgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGc                              <  2:6124818/101‑1 (MQ=255)
                                                                        cgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGc                              >  1:3019490/1‑101 (MQ=255)
                                                                          ggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTc                            <  1:4192287/101‑1 (MQ=255)
                                                                          ggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTc                            <  2:2667296/101‑1 (MQ=255)
                                                                           gggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCg                           >  1:7569269/1‑101 (MQ=255)
                                                                           gggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCg                           >  1:5329492/1‑101 (MQ=255)
                                                                            ggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGt                          <  2:5319858/101‑1 (MQ=255)
                                                                            ggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGt                          >  2:7804770/1‑101 (MQ=255)
                                                                               gTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCa                       <  2:8986326/101‑1 (MQ=255)
                                                                               gTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCa                       <  2:4360842/101‑1 (MQ=255)
                                                                                 ggATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGc                     <  1:1329052/101‑1 (MQ=255)
                                                                                 ggATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGc                     >  1:770838/1‑101 (MQ=255)
                                                                                   aTCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCCGCCg                   >  2:4840182/1‑101 (MQ=255)
                                                                                     cGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGc                 >  1:4000803/1‑101 (MQ=255)
                                                                                          cATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGcc            >  2:2022417/1‑100 (MQ=255)
                                                                                          cATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCt            <  1:221034/101‑1 (MQ=255)
                                                                                           aTCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtc           <  2:5797300/101‑1 (MQ=255)
                                                                                             cgccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgt         <  1:5878758/101‑1 (MQ=255)
                                                                                              gccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgtc        >  2:1006205/1‑101 (MQ=255)
                                                                                              gccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgtc        <  2:8532970/101‑1 (MQ=255)
                                                                                              gccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgtc        >  1:138237/1‑101 (MQ=255)
                                                                                               cctatttGGCGGGCGGCAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgtcg       <  2:8141375/97‑1 (MQ=255)
                                                                                               ccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgtcg       <  2:1810746/101‑1 (MQ=255)
                                                                                                 gcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCt     <  2:8040966/101‑1 (MQ=255)
                                                                                                  cTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTc    >  1:3797134/1‑101 (MQ=255)
                                                                                                   tttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTcc   >  1:6572113/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CCGCGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCC  >  NC_014034/1395375‑1395575

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG