Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,011 Δ1 bp 13.0% intergenic (+823/‑108) 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,396,0110C→.13.0% 93.8 ‑0.4 ‑0.4231intergenic (+823/‑108)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (13/17):  ref base (96/105):  total (109/122)
Fisher's exact test strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.58e-01

Only 100 of 261 total aligned reads displayed.

GCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCG  >  NC_014034/1395911‑1396111
                                                                                                    |                                                                                                       
  gggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGTGGGGGCGTGGATCGAAACTCGCCCTCGATCCgcgc                                                                                                       >  1:746439/1‑101 (MQ=255)
  gggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCgcgc                                                                                                       <  2:2066226/101‑1 (MQ=255)
  gggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCgcgc                                                                                                       <  1:2441131/101‑1 (MQ=255)
   ggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCt                                                                                                      >  1:1130848/1‑101 (MQ=255)
   ggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCt                                                                                                      >  1:3668144/1‑101 (MQ=255)
   ggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCt                                                                                                      <  1:2793985/101‑1 (MQ=255)
    gggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCtt                                                                                                     <  2:499114/101‑1 (MQ=255)
    gggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCtt                                                                                                     >  1:3462238/1‑101 (MQ=255)
    gggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCtt                                                                                                     <  2:3106823/101‑1 (MQ=255)
     ggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTg                                                                                                    <  2:1470601/101‑1 (MQ=255)
     ggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTAGCTCCCCCAGAGGGGGCATGGATCGAAAATCGCCCTACATCAGCGCTTg                                                                                                    >  1:1532505/1‑101 (MQ=255)
         tGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                                                                                                >  1:108379/1‑101 (MQ=255)
         tGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                                                                                                <  1:2701909/101‑1 (MQ=255)
             tCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGc                                                                                            <  1:2015402/101‑1 (MQ=255)
              cGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGcc                                                                                           >  2:1042557/1‑101 (MQ=255)
              cGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGcc                                                                                           >  2:254475/1‑101 (MQ=255)
               gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                                                                                          <  2:3342911/101‑1 (MQ=255)
                aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTg                                                                                         >  1:2084398/1‑101 (MQ=255)
                aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACCCCtc                                                                                         >  1:4041580/1‑100 (MQ=255)
                 aaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGc                                                                                        >  1:4653511/1‑101 (MQ=255)
                 aaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGc                                                                                        <  2:3003347/101‑1 (MQ=255)
                   cGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.cg                                                                                     <  2:2567666/101‑1 (MQ=255)
                       gCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTc                                                                               <  2:743191/101‑1 (MQ=255)
                       gCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..CTc                                                                               <  2:270055/101‑1 (MQ=255)
                        cGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCg                                                                              <  2:1699253/101‑1 (MQ=255)
                          ttCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCt                                                                            >  1:2514409/1‑101 (MQ=255)
                             aTGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTccc                                                                         <  2:1540790/101‑1 (MQ=255)
                              tGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcccc                                                                        >  1:3320483/1‑101 (MQ=255)
                                 ccTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCcg                                                                     <  2:3607987/101‑1 (MQ=255)
                                   tACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCcgcg                                                                   <  1:2550393/101‑1 (MQ=255)
                                    aCGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCgg                                                                  >  1:3145741/1‑101 (MQ=255)
                                      gTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCCCACGCCTGC.CGC..GTCGCTCCCCCCGCgggg                                                                <  1:2868459/101‑1 (MQ=255)
                                        cGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGc                                                              >  1:4064866/1‑101 (MQ=255)
                                         gggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCg                                                             <  2:3790043/101‑1 (MQ=255)
                                         gggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCg                                                             <  2:4041295/101‑1 (MQ=255)
                                          ggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGt                                                            <  2:4563317/101‑1 (MQ=255)
                                            gaGTGGCCTGGTCGCCCCCCCCGCGGGGGCNTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTgg                                                          <  1:1371432/100‑1 (MQ=255)
                                            cAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTgg                                                          <  2:215747/101‑1 (MQ=255)
                                            cAGTGGACTGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTgg                                                          <  1:3190873/101‑1 (MQ=255)
                                              tctgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGa                                                         <  1:3425018/94‑1 (MQ=2)
                                              gTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGAt                                                        >  1:3576604/1‑101 (MQ=255)
                                               tGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATc                                                       <  2:1905792/101‑1 (MQ=255)
                                               tGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATc                                                       <  1:1288211/101‑1 (MQ=255)
                                               tGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATc                                                       <  1:4469809/101‑1 (MQ=255)
                                                 gcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                      <  2:1002675/97‑1 (MQ=2)
                                                   gcggTCGCTTCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                    >  2:4461188/3‑101 (MQ=2)
                                                   gcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                    >  2:935276/3‑101 (MQ=2)
                                                   cTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                   <  1:919845/101‑1 (MQ=255)
                                                    tGGTCGCTCCCCCGGCGGGGGCGTGGATCGAAACTCGCCCTCCATCCACGCTTGCCTCACGCCTGC.CCC..GTCGCTCCACCCGCGGGGGCGTGGagctgaac                                                  >  2:405819/1‑94 (MQ=37)
                                                    tGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                  >  1:958748/1‑101 (MQ=255)
                                                      gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                 <  2:4136029/101‑2 (MQ=2)
                                                      gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAg                                                >  1:1779890/1‑101 (MQ=255)
                                                       gcGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgg                                               <  2:2940811/100‑1 (MQ=255)
                                                        cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtg                                               >  1:1243905/1‑98 (MQ=2)
                                                           tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatg                                            >  1:4586988/1‑95 (MQ=2)
                                                           tCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCTCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCa                                           <  1:2909626/101‑1 (MQ=255)
                                                           tCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCa                                           <  2:1262980/101‑1 (MQ=255)
                                                            ccccccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatga                                           <  1:4268064/101‑8 (MQ=2)
                                                            ccccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCAt                                          >  2:288134/1‑101 (MQ=255)
                                                            ccccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCAt                                          <  1:2190562/101‑1 (MQ=255)
                                                            ccccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCAt                                          <  2:3277937/101‑1 (MQ=255)
                                                            ccccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCAt                                          <  2:1958666/101‑1 (MQ=255)
                                                             cccccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatgac                                          >  2:3520337/1‑93 (MQ=2)
                                                              ccccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatgacc                                         <  2:1370798/101‑10 (MQ=2)
                                                              ccccGCGGGGGCGTGGATCGAAACTCGCCCTAGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCt                                        >  2:626216/1‑101 (MQ=255)
                                                               cccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatgacca                                        >  2:279094/1‑91 (MQ=2)
                                                                 cgcgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGAt                                     <  1:37383/101‑1 (MQ=255)
                                                                 cgcgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGAt                                     <  2:2431869/101‑1 (MQ=255)
                                                                  gcgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATg                                    <  2:2084258/101‑1 (MQ=255)
                                                                   cgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGt                                   >  1:4477014/1‑101 (MQ=255)
                                                                   cgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGt                                   >  2:1232704/1‑101 (MQ=255)
                                                                    gggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGta                                  >  2:2617216/1‑101 (MQ=255)
                                                                     ggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtat                                 >  1:1646869/1‑101 (MQ=255)
                                                                      gggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..TTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtata                                <  1:150192/101‑1 (MQ=255)
                                                                      gggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtata                                <  2:2305256/101‑1 (MQ=255)
                                                                      gggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtata                                <  1:261038/101‑1 (MQ=255)
                                                                          gTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGcc                            >  2:700290/1‑101 (MQ=255)
                                                                          gTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGcc                            <  1:1270220/101‑1 (MQ=255)
                                                                           tGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGccc                           <  2:215449/101‑1 (MQ=255)
                                                                            ggATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCg                          <  1:1940198/101‑1 (MQ=255)
                                                                             gATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCgg                         >  2:4355629/1‑101 (MQ=255)
                                                                                 gAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTgag                     <  2:4214528/101‑1 (MQ=255)
                                                                                    aCTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACg                  <  2:3919809/101‑1 (MQ=255)
                                                                                     cTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACgg                 >  2:3330220/1‑101 (MQ=255)
                                                                                       cGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCt               <  2:1646758/101‑1 (MQ=255)
                                                                                        gCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTg              >  1:2236078/1‑101 (MQ=255)
                                                                                           cTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCg           <  1:2277716/101‑1 (MQ=255)
                                                                                             cGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGACACAGGCGCATCTCATGTATATGCCCGGTGGGACGGCtcgcgcc         >  2:3870468/1‑95 (MQ=38)
                                                                                             cGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCt         <  1:1809717/101‑1 (MQ=255)
                                                                                             cGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCt         >  1:2161361/1‑101 (MQ=255)
                                                                                             cGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCt         >  1:2220908/1‑101 (MQ=255)
                                                                                             cGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCt         >  1:802656/1‑101 (MQ=255)
                                                                                              gATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTc        <  1:626277/101‑1 (MQ=255)
                                                                                               aTCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcc       <  1:1252636/101‑1 (MQ=255)
                                                                                                tCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTccc      >  1:3286095/1‑101 (MQ=255)
                                                                                                 ccGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccc     <  2:1130759/101‑1 (MQ=255)
                                                                                                  cgcgCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTccccc    >  1:4359143/1‑101 (MQ=255)
                                                                                                   gcgcTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc   >  1:3865750/1‑101 (MQ=255)
                                                                                                   gcgcTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc   <  2:4476719/101‑1 (MQ=255)
                                                                                                   gcgcTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc   >  1:4090804/1‑101 (MQ=255)
                                                                                                    |                                                                                                       
GCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCG  >  NC_014034/1395911‑1396111

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG