Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,672 T→A 8.8% intergenic (+484/‑447) 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,6720T→A 8.8% 49.3 ‑1.8 ‑0.1241intergenic (+484/‑447)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/13):  ref base (121/97):  total (130/111)
Fisher's exact test strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.48e-01

Only 100 of 260 total aligned reads displayed.

TCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCC  >  NC_014034/1395572‑1395772
                                                                                                    |                                                                                                       
  cccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGt                                                                                                       <  1:656153/101‑1 (MQ=255)
  cccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGt                                                                                                       <  1:2248046/101‑1 (MQ=255)
   ccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTa                                                                                                      <  1:2645655/101‑1 (MQ=255)
    cccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAg                                                                                                     <  1:4438630/101‑1 (MQ=255)
     ccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGt                                                                                                    >  1:4588781/1‑101 (MQ=255)
     ccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGt                                                                                                    <  1:3548680/101‑1 (MQ=255)
      cgcgGGGGCGTGGATCGAAGCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTc                                                                                                   <  1:2499542/101‑1 (MQ=255)
      cgcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTc                                                                                                   <  2:2300482/101‑1 (MQ=255)
        cgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCGACATGTGGTcgg                                                                                                 >  2:642142/1‑101 (MQ=38)
         gggggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggc                                                                                                >  2:2679407/1‑101 (MQ=255)
           gggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggcgg                                                                                              >  1:5004877/1‑101 (MQ=255)
            ggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCgg..g                                                                                           >  1:2617682/1‑101 (MQ=255)
             gCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..Gc                                                                                          <  2:222938/101‑1 (MQ=255)
              cGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCa                                                                                         >  2:299656/1‑101 (MQ=255)
                tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGt                                                                                       >  1:1488893/1‑101 (MQ=255)
                 ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGtt                                                                                      <  2:1509675/101‑1 (MQ=255)
                 ggATCGAAAACACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGtt                                                                                      >  2:3704904/1‑101 (MQ=255)
                   aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCg                                                                                    >  2:133545/1‑101 (MQ=255)
                   aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCg                                                                                    <  1:1499738/101‑1 (MQ=255)
                   aTCGAAACCACCTTGCGGATTGTCACNNNCNCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCg                                                                                    <  1:2696737/101‑1 (MQ=255)
                      taaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.Gc                                                                                <  2:2145235/100‑1 (MQ=255)
                      gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.Gc                                                                                <  1:522476/101‑1 (MQ=255)
                      gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.Gc                                                                                >  2:3124007/1‑101 (MQ=255)
                      gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.Gc                                                                                >  2:2026657/1‑101 (MQ=255)
                      gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.Gc                                                                                >  1:1563443/1‑101 (MQ=255)
                       aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.Gcc                                                                               <  2:725296/101‑1 (MQ=255)
                       aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.Gca                                                                               >  1:2145383/1‑100 (MQ=255)
                        aaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCGGTTCGA.GCCg                                                                              >  2:2553771/1‑101 (MQ=255)
                         accaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGc                                                                             >  1:4538627/1‑101 (MQ=255)
                         accaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGc                                                                             >  1:4305370/1‑101 (MQ=255)
                             ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTc                                                                         <  2:5004547/101‑1 (MQ=255)
                              cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGCTCGA.GCCGCAGTCg                                                                        >  1:4360436/1‑101 (MQ=255)
                                tGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCt                                                                      <  2:2908564/101‑1 (MQ=255)
                                 gCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTc                                                                     >  2:4212090/1‑101 (MQ=255)
                                  cGGATTGTCACTTTCGCGGGGCGAGGCGTCGCACCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGCCGG..GCAGTTCGA.GCCGCAGTCGCTcc                                                                    <  2:4747033/101‑1 (MQ=255)
                                   ggATTGTCCCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTccc                                                                   <  1:119712/101‑1 (MQ=255)
                                    gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTcccc                                                                  >  2:1759941/1‑101 (MQ=255)
                                      ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTcccccc                                                                >  2:837391/1‑101 (MQ=255)
                                      ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTcccccc                                                                <  2:4331844/101‑1 (MQ=255)
                                       tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                               <  2:1076791/101‑1 (MQ=255)
                                       tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                               >  2:1144047/1‑101 (MQ=255)
                                       tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                               >  2:2108197/1‑101 (MQ=255)
                                        gTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                              <  2:4026298/101‑1 (MQ=255)
                                           ccTTTCGCGGGGCGAGGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggg                                                           <  2:4443916/100‑1 (MQ=255)
                                             tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                         >  1:1310371/1‑101 (MQ=255)
                                             tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                         >  2:4564534/1‑101 (MQ=255)
                                              ttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGc                                                        >  1:1020281/1‑101 (MQ=255)
                                                 tggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGt                                                      <  1:31130/91‑1 (MQ=2)
                                                 gcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTg                                                     <  2:2207154/101‑1 (MQ=255)
                                                  gggggCGAAGCGTCGCCCCCCCCGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                    <  1:8585/100‑1 (MQ=255)
                                                  gggggCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                    <  1:1100674/100‑1 (MQ=255)
                                                   ggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGa                                                   <  1:3762838/101‑1 (MQ=255)
                                                   ggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGa                                                   >  1:2585481/1‑101 (MQ=255)
                                                     ggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATTTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                 <  2:3336160/101‑1 (MQ=255)
                                                     gccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                  >  1:1698323/7‑101 (MQ=2)
                                                      ccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                 <  2:1007073/96‑1 (MQ=2)
                                                       cGAAGCGTCGCTCCCGCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTTGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGGTCGa                                               >  1:4260684/1‑101 (MQ=255)
                                                        gctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                               >  1:4118330/4‑101 (MQ=2)
                                                          aGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                            >  1:2233681/1‑101 (MQ=255)
                                                          aGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                            >  2:393591/1‑101 (MQ=255)
                                                          aGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                            >  2:3435013/1‑101 (MQ=255)
                                                           gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGGAACt                                           >  2:3517962/1‑101 (MQ=255)
                                                           gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                            >  2:2608704/1‑101 (MQ=2)
                                                              ccGCGCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGt                                        <  1:3641943/100‑1 (MQ=255)
                                                                  tCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGGGGGGGTCGAAAcacccatg                                     >  1:1714131/1‑94 (MQ=2)
                                                                  cccccccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                    <  2:4906255/100‑1 (MQ=255)
                                                                   ccccccGCGGGGGCGTGGATCGAAACATGGCTATGGCGGCAGCGGCAGGGTTTCCA.GCGGC.GTCCCTCCCCCCCCGGGGGCGGGGATCGAAAcccgccccg                                    >  1:2049403/1‑93 (MQ=0)
                                                                    cccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGAt                                  <  1:2087672/101‑1 (MQ=255)
                                                                     ccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATc                                 <  1:340998/101‑1 (MQ=255)
                                                                     ccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATc                                 <  2:1935384/101‑1 (MQ=255)
                                                                        cgcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTgg                              >  2:2836856/1‑101 (MQ=255)
                                                                         gcgGGGGCGGGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                             <  2:325383/101‑1 (MQ=255)
                                                                         gcaggggCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                             >  1:2892064/4‑101 (MQ=255)
                                                                          ggggggCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcg                            <  1:3230763/100‑1 (MQ=255)
                                                                           gggggGCTGGACCGAAATATCTGCGTGTAGTCGGGGG..GCTGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcgc                           <  2:3835088/101‑1 (MQ=25)
                                                                              ggCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaaa                        >  1:728828/1‑101 (MQ=255)
                                                                              ggCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaaa                        >  2:521171/1‑101 (MQ=255)
                                                                                  tGGCTCGAAACATCTGCTTGTAGGCGGCGG..GGAGTGGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGAGCTGGCGCAAACGGa                    <  1:334782/101‑1 (MQ=25)
                                                                                  tGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                    >  1:4838796/1‑101 (MQ=255)
                                                                                   ggATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa                   >  2:2285565/1‑101 (MQ=255)
                                                                                   ggATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa                   <  1:4171172/101‑1 (MQ=255)
                                                                                   ggATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa                   >  2:170361/1‑101 (MQ=255)
                                                                                    gATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa                  >  2:3603901/1‑101 (MQ=255)
                                                                                    gATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa                  >  2:3665253/1‑101 (MQ=255)
                                                                                    gATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa                  <  1:2981204/101‑1 (MQ=255)
                                                                                       cGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc               >  2:2767760/1‑101 (MQ=255)
                                                                                        gAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt              >  2:1896613/1‑101 (MQ=255)
                                                                                         aaaCATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg             >  1:2734176/1‑101 (MQ=255)
                                                                                         aaaCATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg             >  1:3585406/1‑101 (MQ=255)
                                                                                         aaaCATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg             <  1:220889/101‑1 (MQ=255)
                                                                                         aaaCATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg             <  1:89249/101‑1 (MQ=255)
                                                                                         aaaCATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg             <  2:4403571/101‑1 (MQ=255)
                                                                                          aaCATCTGCATGTGGTGGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt            <  1:2106180/101‑1 (MQ=255)
                                                                                           aCATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtc           >  1:1928921/1‑101 (MQ=255)
                                                                                            cATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcg          >  1:2817565/1‑101 (MQ=255)
                                                                                            cATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcg          >  2:4155464/1‑101 (MQ=255)
                                                                                             aTCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgt         <  2:1616498/101‑1 (MQ=255)
                                                                                               cTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGtcgtcg       >  1:2398007/1‑101 (MQ=255)
                                                                                                   aTGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcc   >  1:3816673/1‑101 (MQ=255)
                                                                                                   aTGTAGTCGGCGG..GAAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcc   <  1:3124245/101‑1 (MQ=255)
                                                                                                    |                                                                                                       
TCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCC  >  NC_014034/1395572‑1395772

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG