Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,670 C→G 8.7% intergenic (+482/‑449) 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,6700C→G 8.7% 41.8 ‑1.8 ‑0.1243intergenic (+482/‑449)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/13):  ref base (122/99):  total (131/112)
Fisher's exact test strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.76e-01

Only 100 of 266 total aligned reads displayed.

GCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGC  >  NC_014034/1395570‑1395768
                                                                                                    |                                                                                                     
gCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGc                                                                                                       >  1:3744950/1‑101 (MQ=255)
gCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGc                                                                                                       >  2:3610793/1‑101 (MQ=255)
gCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAATCTGc                                                                                                       >  1:812632/1‑101 (MQ=255)
 cTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCCAAGCGTCGCTCCCCCCGCGGGGGAGTGGCTGGAAACCTCTGCa                                                                                                      >  2:4118483/1‑101 (MQ=38)
     ccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTa                                                                                                  <  1:2645655/101‑1 (MQ=255)
      cccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAg                                                                                                 <  1:4438630/101‑1 (MQ=255)
       ccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGt                                                                                                <  1:3548680/101‑1 (MQ=255)
        cgcgGGGGCGTGGATCGAAGCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTc                                                                                               <  1:2499542/101‑1 (MQ=255)
         gcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcg                                                                                              >  2:1744567/1‑101 (MQ=255)
          cgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcgg                                                                                             <  1:2212940/101‑1 (MQ=255)
          cgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCGACATGTGGTcgg                                                                                             >  2:642142/1‑101 (MQ=38)
           gggggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggc                                                                                            >  2:3481978/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)
                 gTGGATCGAAGCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAg                                                                                    >  1:4364701/1‑101 (MQ=255)
                  tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGt                                                                                   >  1:1488893/1‑101 (MQ=255)
                  tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAgg                                                                                   <  2:2370631/101‑2 (MQ=255)
                   ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGtt                                                                                  <  2:1509675/101‑1 (MQ=255)
                   ggATCGAAAACACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGtt                                                                                  >  2:3704904/1‑101 (MQ=255)
                     aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCg                                                                                <  1:1499738/101‑1 (MQ=255)
                     aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCg                                                                                >  2:133545/1‑101 (MQ=255)
                     aTCGAAACCACCTTGCGGATTGTCACNNNCNCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCg                                                                                <  1:2696737/101‑1 (MQ=255)
                        taaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.Gc                                                                            <  2:2145235/100‑1 (MQ=255)
                        gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.Gc                                                                            >  1:1563443/1‑101 (MQ=255)
                         aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.Gcc                                                                           <  2:725296/101‑1 (MQ=255)
                           accaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGc                                                                         >  2:4399477/1‑101 (MQ=255)
                           accaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGc                                                                         >  1:4305370/1‑101 (MQ=255)
                              cccTTGCGGATTGTCACTTTCGCGGGGCGAAGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGt                                                                      <  1:2087942/100‑1 (MQ=255)
                                cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCg                                                                    >  2:3218459/1‑101 (MQ=255)
                                  tGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCt                                                                  <  2:2908564/101‑1 (MQ=255)
                                   gCGGATTGTCTCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTc                                                                 <  1:2165062/101‑1 (MQ=255)
                                   gCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTc                                                                 >  2:4212090/1‑101 (MQ=255)
                                   gCGGATTGCGACTTTGGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAGCATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTc                                                                 >  2:334729/1‑101 (MQ=38)
                                    cGGATTGTCACTTTCGCGGGGCGAGGCGTCGCACCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGCCGG..GCAGTTCGA.GCCGCAGTCGCTcc                                                                <  2:4747033/101‑1 (MQ=255)
                                     ggATTGTCCCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTccc                                                               <  1:119712/101‑1 (MQ=255)
                                     ggATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTccc                                                               >  1:2291507/1‑101 (MQ=255)
                                        ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTcccccc                                                            >  2:837391/1‑101 (MQ=255)
                                         tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                           >  2:2108197/1‑101 (MQ=255)
                                         tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                           <  2:2587669/101‑1 (MQ=255)
                                         tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                           >  2:1144047/1‑101 (MQ=255)
                                         tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                           <  2:1076791/101‑1 (MQ=255)
                                          gTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                          <  2:4026298/101‑1 (MQ=255)
                                           tCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcgcg                                                         >  1:4180695/1‑101 (MQ=255)
                                             ccTTTCGCGGGGCGAGGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggg                                                       <  2:4443916/100‑1 (MQ=255)
                                               tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                     >  1:1310371/1‑101 (MQ=255)
                                                ttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGc                                                    >  1:1020281/1‑101 (MQ=255)
                                                  cgcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGt                                                  >  1:4545643/1‑101 (MQ=255)
                                                    gggggCGAAGCGTCGCCCCCCCCGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                <  1:8585/100‑1 (MQ=255)
                                                    gggggCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                <  1:1100674/100‑1 (MQ=255)
                                                    cgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                >  2:4949503/1‑101 (MQ=255)
                                                     gtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                <  1:2882074/93‑1 (MQ=2)
                                                     ggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGa                                               >  1:2585481/1‑101 (MQ=255)
                                                       ggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                             >  1:3128158/1‑101 (MQ=255)
                                                       gccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                              >  1:1698323/7‑101 (MQ=2)
                                                        gCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                            <  1:720017/101‑1 (MQ=255)
                                                        gCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                            >  1:3793015/1‑101 (MQ=255)
                                                        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)
                                                           ctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                          <  2:2016318/99‑1 (MQ=2)
                                                            aGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                        >  2:393591/1‑101 (MQ=255)
                                                            aGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                        >  1:2233681/1‑101 (MQ=255)
                                                             gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGGAACt                                       >  2:3517962/1‑101 (MQ=255)
                                                             gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                       <  1:59184/101‑1 (MQ=255)
                                                              cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                      >  2:4326350/1‑101 (MQ=255)
                                                               gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgg                                     >  2:2690386/1‑101 (MQ=255)
                                                                ccGCGCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGt                                    <  1:3641943/100‑1 (MQ=255)
                                                                   cTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATg                                 >  1:4599649/1‑101 (MQ=255)
                                                                    cccccccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGGGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                <  1:1760071/100‑1 (MQ=255)
                                                                    cccccccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                <  2:4906255/100‑1 (MQ=255)
                                                                    cccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                <  2:216899/100‑1 (MQ=255)
                                                                     ccccccGCGGGGGCGTGGATCGAAACATGGCTATGGCGGCAGCGGCAGGGTTTCCA.GCGGC.GTCCCTCCCCCCCCGGGGGCGGGGATCGAAAcccgccccg                                >  1:2049403/1‑93 (MQ=0)
                                                                       ccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATc                             <  1:340998/101‑1 (MQ=255)
                                                                       ccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATc                             <  2:1935384/101‑1 (MQ=255)
                                                                        cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCt                            <  1:2805292/101‑1 (MQ=255)
                                                                           gcgGGGGCGGGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                         <  2:325383/101‑1 (MQ=255)
                                                                            ggggggCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcg                        <  1:3230763/100‑1 (MQ=255)
                                                                             gggggGCTGGACCGAAATATCTGCGTGTAGTCGGGGG..GCTGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcgc                       <  2:3835088/101‑1 (MQ=25)
                                                                             gggggCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcgc                       <  1:5029480/101‑1 (MQ=255)
                                                                               gggCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaa                     <  1:3254446/101‑1 (MQ=255)
                                                                                ggggTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaaa                    <  1:3068081/98‑1 (MQ=255)
                                                                                ggCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaaa                    >  2:521171/1‑101 (MQ=255)
                                                                                    tGGATCGAAACATCTGCATGTAGTGGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                <  2:1860631/101‑1 (MQ=255)
                                                                                    tGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                >  1:4838796/1‑101 (MQ=255)
                                                                                     ggATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa               >  2:170361/1‑101 (MQ=255)
                                                                                     ggATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa               >  2:2285565/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:4527218/1‑101 (MQ=255)
                                                                                         cGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc           >  2:2767760/1‑101 (MQ=255)
                                                                                           aaaCATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg         >  1:3585406/1‑101 (MQ=255)
                                                                                           aaaCATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg         <  1:89249/101‑1 (MQ=255)
                                                                                           aaaCATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg         >  1:2734176/1‑101 (MQ=255)
                                                                                            aaCATCTGCATGTGGTGGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt        <  1:2106180/101‑1 (MQ=255)
                                                                                             aCATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtc       >  1:1928921/1‑101 (MQ=255)
                                                                                               aTCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgt     <  2:1616498/101‑1 (MQ=255)
                                                                                                tCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc    >  2:2678765/1‑101 (MQ=255)
                                                                                                 cTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGtcgtcg   >  1:2398007/1‑101 (MQ=255)
                                                                                                  tGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGc  >  2:2857642/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
GCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGC  >  NC_014034/1395570‑1395768

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