Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,964 G→C 3.0% intergenic (+776/‑155) 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,9640G→C 3.0% 4.6 ‑3.5 ‑0.0195intergenic (+776/‑155)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (96/93):  total (96/99)
Fisher's exact test strand distribution p-value = 2.90e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.68e-01

Only 100 of 236 total aligned reads displayed.

TGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA  >  NC_014034/1395864‑1396064
                                                                                                    |                                                                                                        
  cgccgcCTATCGCATCCGGTCGAGGGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGt                                                                                                        <  1:3892858/101‑1 (MQ=255)
  cgccgcCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGt                                                                                                        <  2:526242/101‑1 (MQ=255)
   gccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTc                                                                                                       <  2:1291395/101‑1 (MQ=255)
   gccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTc                                                                                                       <  1:4556401/101‑1 (MQ=255)
    ccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCg                                                                                                      <  1:4283061/101‑1 (MQ=255)
        cTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcc                                                                                                  >  1:1173263/1‑101 (MQ=255)
            cGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccccc                                                                                              >  2:2567574/1‑101 (MQ=255)
            cGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccccc                                                                                              >  2:387858/1‑101 (MQ=255)
             gCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCcg                                                                                             <  1:2867270/101‑1 (MQ=255)
               atacGGACGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCATGCGCTCACAGCCTACGTCGGGCAGCGGACTGGCCGGTCCCCCcgtg                                                                                           >  2:3918164/4‑99 (MQ=21)
                tCCGGTCGAGCTCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCgg                                                                                          <  2:887438/101‑1 (MQ=255)
                 ccGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggg                                                                                         <  1:3748608/101‑1 (MQ=255)
                  cGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCgggg                                                                                        >  2:1884231/1‑101 (MQ=255)
                   ggTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggggg                                                                                       <  1:3187939/101‑1 (MQ=255)
                   ggTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggggg                                                                                       >  1:3066500/1‑101 (MQ=255)
                    gTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGc                                                                                      <  2:627093/101‑1 (MQ=255)
                      cGAGCGCCGCGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.Gt                                                                                   <  1:3272763/101‑1 (MQ=255)
                       gAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTg                                                                                  >  2:1401359/1‑101 (MQ=255)
                         gcgcCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGa                                                                                <  1:2998401/101‑1 (MQ=255)
                             cgcgCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGa                                                                            <  2:3028794/101‑1 (MQ=255)
                              gcgcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGaa                                                                           <  1:2640973/101‑1 (MQ=255)
                              gcgcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGaa                                                                           <  2:410915/101‑1 (MQ=255)
                               cgcAGTCCCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGaaa                                                                          <  2:3354343/101‑1 (MQ=255)
                                gcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAAc                                                                         <  2:243343/101‑1 (MQ=255)
                                  aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGGGGGATCGaaccc                                                                        >  1:686447/1‑98 (MQ=37)
                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGc                                                                     >  2:4139755/1‑101 (MQ=255)
                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GGGGATCGAAACTCGc                                                                     >  1:2131017/1‑101 (MQ=255)
                                        tCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTc                                                                 <  2:86953/101‑1 (MQ=255)
                                         ccccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCg                                                                <  1:1752110/101‑1 (MQ=255)
                                         ccccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGATCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCg                                                                >  1:3548071/1‑101 (MQ=255)
                                          cccccGGGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGa                                                               <  2:3238741/101‑1 (MQ=255)
                                          cccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGa                                                               <  2:3111619/101‑1 (MQ=255)
                                             ccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATcc                                                            <  1:4581888/101‑1 (MQ=255)
                                             ccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATcc                                                            <  1:41308/101‑1 (MQ=255)
                                              cgcgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCcg                                                           <  2:3550244/101‑1 (MQ=255)
                                               gcgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCcgc                                                          >  2:750351/1‑101 (MQ=255)
                                                cgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCTCCCTCGGTCcgcc                                                         >  1:815301/1‑100 (MQ=255)
                                                cgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCcgcg                                                         >  1:3500657/1‑101 (MQ=255)
                                                cgGGGGCGTGGATCGAAACGATGCGTTCAAGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCcgcg                                                         >  2:111859/1‑101 (MQ=255)
                                                  ggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCt                                                       <  2:2900682/101‑1 (MQ=255)
                                                   gggGGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCtt                                                      <  2:4662317/101‑1 (MQ=255)
                                                    ggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTg                                                     >  1:4470961/1‑101 (MQ=255)
                                                       gTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCt                                                  >  2:1813916/1‑101 (MQ=255)
                                                       gTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCt                                                  >  1:2018074/1‑101 (MQ=255)
                                                        tGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                                                 >  2:3736250/1‑101 (MQ=255)
                                                         ggATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCa                                                <  2:2987414/101‑1 (MQ=255)
                                                          gATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTTGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCAc                                               <  2:119661/101‑1 (MQ=255)
                                                          gATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCAc                                               <  2:4380455/101‑1 (MQ=255)
                                                            tCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTccctc                                             >  1:1084115/1‑97 (MQ=255)
                                                             cGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGcc                                            >  2:1004476/1‑101 (MQ=255)
                                                             cGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCCCGATCCGCGCTTGCCTCCCGcc                                            >  1:1202839/1‑101 (MQ=255)
                                                              gAAACGATGCGTTCATGGGCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                                           <  1:4731998/101‑1 (MQ=255)
                                                              gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCCCGCCt                                           >  1:2060740/1‑101 (MQ=255)
                                                              gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                                           >  2:800298/1‑101 (MQ=255)
                                                               aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTg                                          >  2:830817/1‑101 (MQ=255)
                                                                 aCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGc.c                                       <  2:3143583/101‑1 (MQ=255)
                                                                  cGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.cg                                      <  1:1897519/101‑1 (MQ=255)
                                                                  cGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.cg                                      >  1:3112429/1‑101 (MQ=255)
                                                                      gCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTc                                <  2:2312661/101‑1 (MQ=255)
                                                                       cGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCg                               <  2:1088789/101‑1 (MQ=255)
                                                                         ttCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGCTCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCt                             <  2:286958/101‑1 (MQ=255)
                                                                          tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTc                            >  1:614023/1‑101 (MQ=255)
                                                                          tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTc                            >  1:3346228/1‑101 (MQ=255)
                                                                          tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTc                            >  1:3875343/1‑101 (MQ=255)
                                                                           cATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcc                           >  1:1342497/1‑101 (MQ=255)
                                                                           cATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcc                           >  1:859291/1‑101 (MQ=255)
                                                                           cATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcc                           >  1:4407172/1‑101 (MQ=255)
                                                                           cATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcc                           <  2:275834/101‑1 (MQ=255)
                                                                              ggCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTccccc                        >  1:832245/1‑101 (MQ=255)
                                                                                  tACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCcgcg                    <  1:4636794/101‑1 (MQ=255)
                                                                                   aCGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCGGCGCTTGCCTCACGCCTGC.CGC..GTCCCCCCCCCCGCgg                   >  1:2289010/1‑101 (MQ=255)
                                                                                   aCGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCgg                   >  2:1053877/1‑101 (MQ=255)
                                                                                    cGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCggg                  >  1:948638/1‑101 (MQ=255)
                                                                                      tCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCggggg                <  2:401538/101‑1 (MQ=255)
                                                                                       cGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGc               <  2:2525020/101‑1 (MQ=255)
                                                                                         ggCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCgg             >  1:2619908/1‑100 (MQ=255)
                                                                                          gCAGTGGCCTGGTCGCTCCCCCCGCGGGGCC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTg            <  2:2509259/101‑1 (MQ=255)
                                                                                          gCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTg            >  1:848836/1‑101 (MQ=255)
                                                                                            aGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGa          >  2:484475/1‑101 (MQ=255)
                                                                                            aGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGa          <  1:177725/101‑1 (MQ=255)
                                                                                             tctgcgcggTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGa          <  2:241647/94‑1 (MQ=2)
                                                                                             gTGGGCTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGAt         <  1:3981755/101‑1 (MQ=255)
                                                                                             gTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGTAt         <  2:3089516/101‑1 (MQ=255)
                                                                                             gTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGAt         >  2:3939381/1‑101 (MQ=255)
                                                                                             gTGGACTGGTCGCCCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGAt         <  1:1004555/101‑1 (MQ=255)
                                                                                              ctgcgcggTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGAt         >  2:1637025/7‑101 (MQ=2)
                                                                                              ctgcgcggTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGAt         <  2:2240109/95‑1 (MQ=9)
                                                                                              ctgcgcggTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGAt         >  1:1834317/7‑101 (MQ=9)
                                                                                               tgcgcggTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATc        <  2:3794143/96‑1 (MQ=9)
                                                                                               ggACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCg       <  1:2328784/101‑1 (MQ=255)
                                                                                                gcgcggTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCg       <  2:43652/97‑1 (MQ=9)
                                                                                                 gcTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCGCCCGCGGGGGCGTGGATCGaa     <  1:2689300/100‑1 (MQ=255)
                                                                                                 cgcggTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGa      <  1:2937090/98‑1 (MQ=2)
                                                                                                 aCTGGTCGCCCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGaa     <  1:2488023/101‑1 (MQ=255)
                                                                                                   tGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc   >  1:2414293/1‑101 (MQ=255)
                                                                                                    ggTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc   <  1:1385337/101‑1 (MQ=2)
                                                                                                    ggTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc   >  1:717343/1‑101 (MQ=9)
                                                                                                    ggTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa  <  2:2222146/101‑1 (MQ=255)
                                                                                                    ggTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa  <  2:1540336/101‑1 (MQ=255)
                                                                                                    ggTCGCCCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa  <  2:1462523/101‑1 (MQ=255)
                                                                                                    |                                                                                                        
TGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA  >  NC_014034/1395864‑1396064

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG