Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,931 G→A 8.9% intergenic (+743/‑188) 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,9310G→A 8.9% 55.3 ‑0.4 ‑0.3279intergenic (+743/‑188)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (12/13):  ref base (134/116):  total (146/133)
Fisher's exact test strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.75e-01

Only 100 of 292 total aligned reads displayed.

CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCG  >  NC_014034/1395831‑1396030
                                                                                                     |                                                                                                   
  ccGCTCCCCCCGCGGGGGCGTTGATCCAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAt                                                                                                   <  2:4637791/100‑1 (MQ=255)
   cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATg                                                                                                  >  2:3158240/1‑101 (MQ=255)
   cGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGCTGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATg                                                                                                  >  2:1288220/1‑101 (MQ=9)
    gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGc                                                                                                 >  1:591223/1‑101 (MQ=255)
    gCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGCTGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGc                                                                                                 <  1:2608907/101‑1 (MQ=9)
    gCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGCTGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGc                                                                                                 >  1:2337891/1‑101 (MQ=9)
    cccccccccGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGc                                                                                                 <  1:3355303/98‑1 (MQ=255)
     ccccccccGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCg                                                                                                <  1:3836971/99‑1 (MQ=255)
     cTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCg                                                                                                >  2:2696962/1‑101 (MQ=255)
     cTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCg                                                                                                >  1:2133240/1‑101 (MQ=255)
       ccccccGCGGGGGCGTGGATCGAAACTGCGCCGGCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGtt                                                                                              <  1:3975749/101‑1 (MQ=255)
       ccccccGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGtt                                                                                              <  2:3624020/101‑1 (MQ=255)
       ccccccGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGtt                                                                                              <  2:3087776/101‑1 (MQ=255)
       caccccGCGGGGCCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGCTGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCggc                                                                                              >  1:2533379/3‑99 (MQ=2)
        cccccGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTc                                                                                             <  2:475371/101‑1 (MQ=255)
         ccccGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCa                                                                                            <  2:1856441/101‑1 (MQ=255)
          cccGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAt                                                                                           <  2:2493510/101‑1 (MQ=255)
           ccGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGCTGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGCGGTGGATCGAAACAATGCggcgcgg                                                                                          >  2:2147493/1‑95 (MQ=2)
           ccGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATg                                                                                          <  1:1742006/101‑1 (MQ=255)
           ccGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATg                                                                                          <  2:2348411/101‑1 (MQ=255)
           ccGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATg                                                                                          <  1:1354207/101‑1 (MQ=255)
            cgcgGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATgg                                                                                         >  2:2050680/1‑101 (MQ=255)
            cgcgGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATgg                                                                                         >  1:387395/1‑101 (MQ=255)
             gcgGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGCTGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCggcgcgggc                                                                                        >  2:4523972/1‑93 (MQ=2)
             gcgGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGc                                                                                        <  2:736191/101‑1 (MQ=255)
              cgGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGcc                                                                                       >  1:3670097/1‑101 (MQ=255)
              cgGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGcc                                                                                       <  2:373830/101‑1 (MQ=255)
               gggggCGTGTATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCt                                                                                      <  2:908516/101‑1 (MQ=255)
               gggggCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGCTGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCggcgcgggcct                                                                                      <  2:2736497/101‑11 (MQ=2)
                 gggCGTGGATCGAATCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTAc                                                                                    <  2:3001182/101‑1 (MQ=255)
                 gggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTAc                                                                                    >  1:2068326/1‑101 (MQ=255)
                  ggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACg                                                                                   >  1:3086284/1‑101 (MQ=255)
                  ggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACg                                                                                   <  2:1436700/101‑1 (MQ=255)
                   gCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGt                                                                                  <  2:2566770/101‑1 (MQ=255)
                    cGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTc                                                                                 >  1:1007820/1‑101 (MQ=255)
                     gTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCg                                                                                >  2:4381946/1‑101 (MQ=255)
                     gTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCg                                                                                >  1:148965/1‑101 (MQ=255)
                      tGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGGTCATGGCCTACGTCgg                                                                               >  2:4322938/1‑101 (MQ=255)
                       ggATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCggg                                                                              >  1:2034581/1‑101 (MQ=255)
                        gATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGc                                                                             <  1:4345755/101‑1 (MQ=255)
                           cGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGt                                                                          >  2:3245083/1‑101 (MQ=255)
                            gAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTg                                                                         >  1:3278182/1‑101 (MQ=255)
                            gAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTg                                                                         <  1:4591890/101‑1 (MQ=255)
                               aCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGAc                                                                      <  1:3578708/101‑1 (MQ=255)
                               aCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGAc                                                                      <  1:2879540/101‑1 (MQ=255)
                                cTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACt                                                                     <  1:3742953/101‑1 (MQ=255)
                                cTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACt                                                                     >  1:1905927/1‑101 (MQ=255)
                                 tGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTg                                                                    >  2:2868243/1‑101 (MQ=255)
                                 tGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTg                                                                    <  1:3108390/101‑1 (MQ=255)
                                  gcgcCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGGGGACTgg                                                                   >  2:1371337/1‑101 (MQ=255)
                                   cgccgcCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGt                                                                  <  1:2753563/101‑1 (MQ=255)
                                      cgccTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGc                                                               <  2:196760/101‑1 (MQ=255)
                                       gccTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCt                                                              <  1:727739/101‑1 (MQ=255)
                                          tATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTccc                                                           >  2:3348191/1‑101 (MQ=255)
                                             cGCATCCGGTCGAGGGCCGCGCA.GTCGGTCCCCCCGCGGGGGCGGGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccccc                                                        <  2:1223730/101‑1 (MQ=255)
                                               cATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCcgc                                                      >  1:1722311/1‑101 (MQ=255)
                                                 tCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCgg                                                    >  1:145263/1‑101 (MQ=255)
                                                  ccGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggg                                                   <  1:4443757/101‑1 (MQ=255)
                                                   cGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCgggg                                                  >  1:121300/1‑101 (MQ=255)
                                                    ggTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggggg                                                 >  1:1164052/1‑101 (MQ=255)
                                                    ggTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggggg                                                 >  2:3587378/1‑101 (MQ=255)
                                                     gTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGc                                                >  2:2190412/1‑101 (MQ=255)
                                                       cGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGt                                              >  2:3451097/1‑101 (MQ=255)
                                                        gAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGGGGACTGGGCGCTCCCCCCGCGGGGGCGTg                                             >  2:3983329/1‑101 (MQ=255)
                                                           cgccgcGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGAt                                          >  2:1573161/1‑101 (MQ=255)
                                                           cgccgcGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGAt                                          <  1:1120712/101‑1 (MQ=255)
                                                            gccgcGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                         <  1:1572664/101‑1 (MQ=255)
                                                            gccgcGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                         >  2:945959/1‑101 (MQ=255)
                                                              cgcgCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                       >  1:28986/1‑101 (MQ=255)
                                                                cgcA.GTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                     <  1:3212863/101‑1 (MQ=255)
                                                                 ggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGCTCGACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                     <  1:2216660/97‑1 (MQ=2)
                                                                    aggcGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTc                                  <  2:2891010/98‑1 (MQ=255)
                                                                        gCTCCCCCCGCGGGGGCGTGGATCTAAACGATGCGTTCATGGCCTACGGCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGccc                              <  1:2705082/101‑1 (MQ=255)
                                                                           ccccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCg                           >  2:261006/1‑101 (MQ=255)
                                                                           ccccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCg                           >  1:2259655/1‑101 (MQ=255)
                                                                            cccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGa                          <  2:3910269/101‑1 (MQ=255)
                                                                             ccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGCCCCCGAt                         >  2:193935/1‑101 (MQ=255)
                                                                              cccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATc                        <  2:2394232/101‑1 (MQ=255)
                                                                                 gcgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgc                     <  1:4116023/101‑1 (MQ=255)
                                                                                  cgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgcg                    >  1:1251371/1‑101 (MQ=255)
                                                                                  cgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgcg                    >  2:653032/1‑101 (MQ=255)
                                                                                   gggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgccg                   >  2:1241600/1‑99 (MQ=255)
                                                                                   gggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCgcgc                   <  2:2066226/101‑1 (MQ=255)
                                                                                   gggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCgcgc                   >  1:1286013/1‑101 (MQ=255)
                                                                                   gggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCgcgc                   <  1:2441131/101‑1 (MQ=255)
                                                                                    ggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCt                  >  2:3317456/1‑101 (MQ=255)
                                                                                     gggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCtt                 >  1:3462238/1‑101 (MQ=255)
                                                                                      ggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTg                >  1:936474/1‑101 (MQ=255)
                                                                                       gCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGc               >  2:1665282/1‑101 (MQ=255)
                                                                                        cGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCCCGATCCGCGCTTGcc              >  2:2909397/1‑101 (MQ=255)
                                                                                           ggATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTcc           >  1:3158860/1‑100 (MQ=255)
                                                                                           ggATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCa           >  2:4229782/1‑101 (MQ=255)
                                                                                           ggATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCa           <  1:782550/101‑1 (MQ=255)
                                                                                             aTCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACg         >  2:3753561/1‑101 (MQ=255)
                                                                                              tCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGc        <  2:2295264/101‑1 (MQ=255)
                                                                                              tCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGc        <  1:2015402/101‑1 (MQ=255)
                                                                                               cGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGcc       >  2:1042557/1‑101 (MQ=255)
                                                                                                 aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTg     >  2:4469514/1‑101 (MQ=255)
                                                                                                  aaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGc    >  1:4653511/1‑101 (MQ=255)
                                                                                                  aaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGc    >  2:2448755/1‑101 (MQ=255)
                                                                                                     |                                                                                                   
CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCG  >  NC_014034/1395831‑1396030

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG