Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,688 T→A 9.6% intergenic (+500/‑431) 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,6880T→A 9.6% 45.6 ‑1.8 ‑0.0221intergenic (+500/‑431)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/13):  ref base (108/90):  total (118/103)
Fisher's exact test strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.46e-01

Only 100 of 261 total aligned reads displayed.

TGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGA  >  NC_014034/1395588‑1395788
                                                                                                       |                                                                                                     
tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGt                                                                                                       >  2:2087524/1‑101 (MQ=255)
 ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGtt                                                                                                      >  1:4216302/1‑101 (MQ=255)
 ggATCGAAAACACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGtt                                                                                                      >  2:3704904/1‑101 (MQ=255)
  gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTc                                                                                                     <  2:3362515/101‑1 (MQ=255)
  gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTc                                                                                                     <  1:4564826/101‑1 (MQ=255)
   aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCg                                                                                                    <  1:1499738/101‑1 (MQ=255)
     cGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.g                                                                                                 >  2:1916211/1‑101 (MQ=255)
      gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.Gc                                                                                                >  1:1563443/1‑101 (MQ=255)
      gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.Gc                                                                                                >  2:3124007/1‑101 (MQ=255)
       aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.Gca                                                                                               >  1:2145383/1‑100 (MQ=255)
        aaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCg                                                                                              <  2:3815001/101‑1 (MQ=255)
        aaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCg                                                                                              <  1:2624974/101‑1 (MQ=255)
         accaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGGG..G.GCAGTTCGA.GCCGc                                                                                             >  1:1953396/1‑101 (MQ=255)
         accaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGc                                                                                             >  1:4538627/1‑101 (MQ=255)
         accaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGc                                                                                             >  2:4399477/1‑101 (MQ=255)
           ccccTTGCGGATTTTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAg                                                                                           <  2:1677199/99‑1 (MQ=255)
            cccTTGCGGATTGTCACTTTCGCGGGGCGAAGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGt                                                                                          <  1:2087942/100‑1 (MQ=255)
             ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTc                                                                                         <  2:5004547/101‑1 (MQ=255)
                tGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCt                                                                                      <  2:2908564/101‑1 (MQ=255)
                 gCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTc                                                                                     >  1:2694637/1‑101 (MQ=255)
                  cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTcc                                                                                    <  1:2733478/101‑1 (MQ=255)
                   ggATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTccc                                                                                   >  1:4403890/1‑101 (MQ=255)
                    gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTcccc                                                                                  >  2:1759941/1‑101 (MQ=255)
                      ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTcccccc                                                                                <  2:4331844/101‑1 (MQ=255)
                       tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                                               >  2:1144047/1‑101 (MQ=255)
                       tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                                               <  2:1076791/101‑1 (MQ=255)
                        gTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                                              >  1:325436/1‑101 (MQ=255)
                           ccTTTCGCGGGGCGAGGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggg                                                                           <  2:4443916/100‑1 (MQ=255)
                            cTTTCGCGGGGCGAGGCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgggg                                                                          <  2:4204907/101‑1 (MQ=255)
                             tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                                         >  1:2062927/1‑101 (MQ=255)
                              ttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGc                                                                        >  1:1020281/1‑101 (MQ=255)
                                 tggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGt                                                                      <  1:31130/91‑1 (MQ=2)
                                 tggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGt                                                                      <  1:293838/91‑1 (MQ=2)
                                  gggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                                    <  2:2146594/100‑1 (MQ=255)
                                  gggggCGAAGCGTCGCCCCCCCCGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                                    <  1:8585/100‑1 (MQ=255)
                                  gggggCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                                    <  1:1100674/100‑1 (MQ=255)
                                   gtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                                    <  1:2882074/93‑1 (MQ=2)
                                   ggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGa                                                                   >  1:2585481/1‑101 (MQ=255)
                                     ggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATTTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                 <  2:3336160/101‑1 (MQ=255)
                                     ggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                 >  1:3128158/1‑101 (MQ=255)
                                     gccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                  >  1:1698323/7‑101 (MQ=2)
                                      gCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                >  1:3793015/1‑101 (MQ=255)
                                      ccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                 <  2:2527541/96‑1 (MQ=2)
                                      ccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                 <  1:2234570/96‑1 (MQ=2)
                                      ccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                 <  2:1007073/96‑1 (MQ=2)
                                       cGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                               >  1:721828/1‑101 (MQ=255)
                                         ctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                              <  1:2523865/99‑1 (MQ=2)
                                         aaGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                             <  1:2306476/101‑1 (MQ=255)
                                          tgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGGTATGGCGTCGGTGTCA.GGTTTTCCA.GCGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                             >  2:804286/2‑101 (MQ=2)
                                          aGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                            >  2:393591/1‑101 (MQ=255)
                                           gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                           <  1:59184/101‑1 (MQ=255)
                                           gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                            >  2:2608704/1‑101 (MQ=2)
                                            ccgcGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                                          <  2:1389154/98‑1 (MQ=255)
                                              tcccTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGt                                                        <  1:3689995/98‑1 (MQ=255)
                                                gCTCCCCCCGCGGGGGCGTGGATCGAAACCCCCGGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcatcga                                                       <  1:4308045/101‑6 (MQ=2)
                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATg                                                     <  2:3657155/101‑1 (MQ=255)
                                                  tCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGGGGGGGTCGAAAcacccatg                                                     >  1:1714131/1‑94 (MQ=2)
                                                  cccccccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                                    <  2:4906255/100‑1 (MQ=255)
                                                  cccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                                    <  2:1599850/100‑1 (MQ=255)
                                                    cccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGAt                                                  <  1:3621762/101‑1 (MQ=255)
                                                     ccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATc                                                 <  1:340998/101‑1 (MQ=255)
                                                      cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCt                                                <  1:2805292/101‑1 (MQ=255)
                                                        cgcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTgg                                              >  2:2836856/1‑101 (MQ=255)
                                                         gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                                             >  1:1235156/1‑101 (MQ=255)
                                                         gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                                             <  1:3567444/101‑1 (MQ=255)
                                                         gcgGGGGCGGGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                                             <  2:325383/101‑1 (MQ=255)
                                                          ggggggCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcg                                            <  1:3230763/100‑1 (MQ=255)
                                                           gggggCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcgc                                           <  1:5029480/101‑1 (MQ=255)
                                                            ggggCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCa                                          >  2:3915629/1‑101 (MQ=255)
                                                             gggCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaa                                         <  2:351360/101‑1 (MQ=255)
                                                              ggggTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaaa                                        <  1:3068081/98‑1 (MQ=255)
                                                               gCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                                       <  2:2212242/101‑1 (MQ=255)
                                                                 gTGGATCGAACCATCTGCATGTAGTGGGCG..G.GCGGTTGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                                     <  2:4046148/101‑1 (MQ=37)
                                                                 gTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAGCGAAACTGGTATGGATCTGGCGCAAACgg                                     >  2:479927/1‑101 (MQ=255)
                                                                  tGGCTCGAAACATCTGCTTGTAGGCGGCG..G.GGAGTGGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGAGCTGGCGCAAACGGa                                    <  1:334782/101‑1 (MQ=25)
                                                                   ggATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa                                   <  1:4171172/101‑1 (MQ=255)
                                                                    gATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa                                  >  1:138248/1‑101 (MQ=255)
                                                                     aTCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAg                                 <  2:1310250/101‑1 (MQ=255)
                                                                     aTCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAg                                 <  2:3590980/101‑1 (MQ=255)
                                                                       cGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc                               >  2:2767760/1‑101 (MQ=255)
                                                                       cGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc                               >  2:4527218/1‑101 (MQ=255)
                                                                        gAAACATCTGCATGTAGTGGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGCCt                              >  2:2138198/1‑101 (MQ=255)
                                                                        gAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt                              >  2:1940270/1‑101 (MQ=255)
                                                                         aaaCATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg                             <  1:220889/101‑1 (MQ=255)
                                                                         aaaCATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg                             <  1:89249/101‑1 (MQ=255)
                                                                          aaCATCTGCATGTGGTGGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt                            <  1:2106180/101‑1 (MQ=255)
                                                                            cATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcg                          <  1:4909450/101‑1 (MQ=255)
                                                                              tCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc                        >  2:2678765/1‑101 (MQ=255)
                                                                                tGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGc                      >  2:2260291/1‑101 (MQ=255)
                                                                                     gTAGTCGGGG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcccc                 <  1:1161287/101‑1 (MQ=255)
                                                                                     gTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcccc                 <  2:4599327/101‑1 (MQ=255)
                                                                                     gTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcccc                 >  1:4734729/1‑101 (MQ=255)
                                                                                     gTAGTCGGCG..G.GAAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcccc                 <  2:4215985/101‑1 (MQ=255)
                                                                                       ggTCGGCG..G.GCAGTCCGA.GCCGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcccccc               <  2:516703/100‑1 (MQ=255)
                                                                                       cggcggcg..g.GCAGTTCGA.CCCGCTGTCGCTGCCCCCGCGGGGGCGTGGATCTAAACTGGTTTGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcccccc               <  1:798152/98‑1 (MQ=38)
                                                                                       aGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcccccc               >  2:55441/1‑101 (MQ=255)
                                                                                              g..g.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCggggg        >  1:4037085/1‑101 (MQ=255)
                                                                                              g..g.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCggggg        >  1:4274480/1‑101 (MQ=255)
                                                                                                 g.gCGGTCCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGc       <  1:3888888/101‑1 (MQ=255)
                                                                                                      gTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTgg   >  2:3049997/1‑101 (MQ=255)
                                                                                                       |                                                                                                     
TGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGA  >  NC_014034/1395588‑1395788

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 40 ≤  ATCG