Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,489 G→T 19.2% intergenic (+301/‑630) 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,4890G→T19.2% 247.1 ‑0.6 ‑0.0474intergenic (+301/‑630)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (42/49):  ref base (190/192):  total (232/242)
Fisher's exact test strand distribution p-value = 5.61e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.96e-01

Only 100 of 491 total aligned reads displayed.

ATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTG  >  NC_014034/1395389‑1395589
                                                                                                       |                                                                                                    
aTCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAgg                                                                                                      >  2:6675991/1‑101 (MQ=255)
aTCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAgg                                                                                                      >  1:7468508/1‑101 (MQ=255)
  cGAAACTGCCCCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCt                                                                                                    <  1:6118568/101‑1 (MQ=255)
  cGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCt                                                                                                    <  1:6329333/101‑1 (MQ=255)
  cGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCt                                                                                                    >  1:1938844/1‑101 (MQ=255)
  cGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCt                                                                                                    <  1:7616801/101‑1 (MQ=255)
   gAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTg                                                                                                   >  1:8533355/1‑101 (MQ=255)
    aaaCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTgg                                                                                                  >  1:5106944/1‑101 (MQ=255)
    aaaCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTgg                                                                                                  <  1:11225/101‑1 (MQ=255)
     aaCTGCCCCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGt                                                                                                 <  1:1293168/101‑1 (MQ=255)
     aaCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGt                                                                                                 <  1:6777684/101‑1 (MQ=255)
     aaCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGt                                                                                                 >  2:4859363/1‑101 (MQ=255)
      aCTGCCCCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTa                                                                                                <  2:7816931/101‑1 (MQ=255)
       cTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTaa                                                                                               <  2:9103021/101‑1 (MQ=255)
       cTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTaa                                                                                               <  1:863149/101‑1 (MQ=255)
         gCACCCGCTGGCCCTGCCGGTGATCGGTTATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGTAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACg                                                                                             >  1:6525000/1‑101 (MQ=255)
         gCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACg                                                                                             >  1:5161286/1‑101 (MQ=255)
          cccccGCTGGCCCCTCCGGTGATCGGTTATCCCGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACgt                                                                                            <  2:4255005/99‑1 (MQ=25)
             ccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTc                                                                                         <  1:7993846/101‑1 (MQ=255)
             ccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTc                                                                                         <  1:5856007/101‑1 (MQ=255)
             ccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTc                                                                                         <  2:129209/101‑1 (MQ=255)
             ccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTc                                                                                         <  1:1006881/101‑1 (MQ=255)
                cTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCt                                                                                      <  2:4712010/101‑1 (MQ=255)
                cTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCt                                                                                      <  1:4094774/101‑1 (MQ=255)
                 tgccccCGCCGGTGATGGGTGATCAGCCGTGGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTc                                                                                     <  1:420131/98‑1 (MQ=37)
                 tGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTc                                                                                     >  1:3557764/1‑101 (MQ=255)
                 tGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGGAACGTGTCGCTc                                                                                     >  2:3216240/1‑101 (MQ=255)
                  ggCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTcc                                                                                    <  1:2769083/101‑1 (MQ=255)
                    cccTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACCTGTCCCTcccc                                                                                  <  2:400051/101‑1 (MQ=255)
                      cTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccccc                                                                                >  2:1690274/1‑101 (MQ=255)
                      cTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccccc                                                                                <  2:8602628/101‑1 (MQ=255)
                       tGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCcg                                                                               <  2:6391801/101‑1 (MQ=255)
                           cgtgatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCggg                                                                           <  1:4983395/100‑1 (MQ=255)
                            gtgatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCgggg                                                                          <  2:5508702/101‑1 (MQ=255)
                            gtgatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCgggg                                                                          >  1:2862823/1‑101 (MQ=255)
                              gatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGc                                                                        <  2:2027817/101‑1 (MQ=255)
                                  ggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTgg                                                                    <  1:4717616/101‑1 (MQ=255)
                                  cgctggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGa                                                                   >  2:5002965/11‑101 (MQ=2)
                                   gtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGa                                                                   >  2:9163434/1‑101 (MQ=255)
                                   gctggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                  <  1:1365540/92‑1 (MQ=2)
                                   gctggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                  <  1:5712136/92‑1 (MQ=2)
                                    ctggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                 >  2:7672361/9‑101 (MQ=0)
                                     tggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                >  2:275063/8‑101 (MQ=2)
                                      ggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                               >  2:4496453/7‑101 (MQ=2)
                                       tcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                               <  1:3119142/101‑1 (MQ=255)
                                        gcgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                             <  2:1596254/97‑1 (MQ=2)
                                         aGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                             >  1:1551480/1‑101 (MQ=255)
                                         aGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                             >  2:7278805/1‑101 (MQ=255)
                                          gtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                           >  2:8739997/3‑101 (MQ=2)
                                          gtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACa                                                           >  1:3628405/3‑101 (MQ=2)
                                            ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                                          <  2:598392/101‑1 (MQ=255)
                                            ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                                          <  1:6027888/101‑1 (MQ=255)
                                            ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                                          >  2:184960/1‑101 (MQ=255)
                                               cgatCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgcttt                                                      <  2:3572540/98‑3 (MQ=2)
                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgcttt                                                      <  1:7803580/101‑3 (MQ=2)
                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtctc                                                       >  1:8168283/1‑101 (MQ=255)
                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtctc                                                       >  2:3349824/1‑101 (MQ=255)
                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtctc                                                       <  2:8469131/101‑1 (MQ=255)
                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttca                                                    <  2:3878235/101‑5 (MQ=2)
                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttca                                                    <  1:3442166/101‑5 (MQ=2)
                                                   ccccccGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaat                                                  <  2:3662026/101‑7 (MQ=0)
                                                   ccccccGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaat                                                  >  1:544274/1‑95 (MQ=0)
                                                   ccccccGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaat                                                  <  1:2902881/101‑7 (MQ=0)
                                                   ccccccGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTg                                                   <  2:1378752/101‑1 (MQ=255)
                                                   ccccccGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTg                                                   >  2:6793337/1‑101 (MQ=255)
                                                     ccccGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGcc                                                 <  2:2313766/101‑1 (MQ=255)
                                                      cccGCGGGGGCGTGGATCGAAACATCGCCG..CT.GTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCg                                                <  1:4551165/101‑1 (MQ=255)
                                                       ccGCGGGGGCGTGGATCGAAACACCCCGGACCTGTCGGCGG.CGA.ACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtcgatgatcaa                                                <  1:6389223/101‑10 (MQ=0)
                                                         gcgGGGGCGGGGATCTAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCg                                             <  2:7602785/101‑1 (MQ=255)
                                                          cgGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGc                                            >  1:2618791/1‑101 (MQ=255)
                                                          cgGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGc                                            <  2:6124818/101‑1 (MQ=255)
                                                            ggggCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTc                                          >  1:2417184/1‑101 (MQ=255)
                                                               gCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTc                                       >  1:4414197/1‑101 (MQ=255)
                                                               gCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTc                                       <  1:8314668/101‑1 (MQ=255)
                                                                 gTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCa                                     <  1:7901781/101‑1 (MQ=255)
                                                                 gTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCa                                     >  1:7112904/1‑101 (MQ=255)
                                                                 gTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCa                                     >  2:8798040/1‑101 (MQ=255)
                                                                 gTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCa                                     <  2:8986326/101‑1 (MQ=255)
                                                                  tGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATTTCGGGGCCGTCGCTCGTCCAg                                    >  1:9129070/1‑101 (MQ=255)
                                                                   ggATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGc                                   >  2:7435126/1‑101 (MQ=255)
                                                                   ggATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGc                                   >  1:770838/1‑101 (MQ=255)
                                                                   ggATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCTCGGTGCCGTCCCTCGTCCCGc                                   >  1:873700/1‑101 (MQ=255)
                                                                       cGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGc                               >  1:6242553/1‑101 (MQ=255)
                                                                       cGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGc                               <  2:5161117/101‑1 (MQ=255)
                                                                        gAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGcc                              >  1:9046347/1‑101 (MQ=255)
                                                                            cATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCt                          <  1:221034/101‑1 (MQ=255)
                                                                             aTCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtc                         <  2:5797300/101‑1 (MQ=255)
                                                                             aTCGCCG..CT.TTGGCGATCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtc                         >  1:7257468/1‑101 (MQ=255)
                                                                               cgccg..cT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgt                       >  1:6912365/1‑101 (MQ=255)
                                                                                gccg..cT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgtc                      >  2:1006205/1‑101 (MQ=255)
                                                                                   g..cT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCt                   <  1:4859523/101‑1 (MQ=255)
                                                                                   g..cT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCt                   <  2:8040966/101‑1 (MQ=255)
                                                                                       t.ttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTcc                 >  1:6572113/1‑101 (MQ=255)
                                                                                           ggCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTccccc              >  1:5322001/1‑101 (MQ=255)
                                                                                             cGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCcg            <  1:1689490/101‑1 (MQ=255)
                                                                                               gTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCcgcg          >  1:3874324/1‑101 (MQ=255)
                                                                                                    cAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGc     <  2:7159632/101‑1 (MQ=255)
                                                                                                     aGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCg    >  1:450028/1‑101 (MQ=255)
                                                                                                       gCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTg  >  1:4016161/1‑101 (MQ=255)
                                                                                                       gCTGGGAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTg  <  2:2984236/101‑1 (MQ=255)
                                                                                                       |                                                                                                    
ATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTG  >  NC_014034/1395389‑1395589

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG