Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,946 T→G 2.9% intergenic (+758/‑173) 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,9460T→G 2.9% 7.8 ‑1.3 ‑0.0241intergenic (+758/‑173)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/5):  ref base (121/112):  total (123/118)
Fisher's exact test strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.37e-01

Only 100 of 263 total aligned reads displayed.

GGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCG  >  NC_014034/1395846‑1396046
                                                                                                    |                                                                                                    
gggggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGCGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGACGCGATCATGGCCt                                                                                                      >  2:4275150/1‑101 (MQ=255)
gggggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCt                                                                                                      <  1:2979617/101‑1 (MQ=255)
gggggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCt                                                                                                      >  2:1140408/1‑101 (MQ=255)
 gtggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTa                                                                                                     <  1:1974485/99‑1 (MQ=255)
   ggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACg                                                                                                   <  1:352989/101‑1 (MQ=255)
   ggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACg                                                                                                   <  1:427807/101‑1 (MQ=255)
   ggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCCACg                                                                                                   >  2:888358/1‑101 (MQ=255)
      gTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCg                                                                                                >  1:860014/1‑101 (MQ=255)
       tGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCgg                                                                                               >  1:5018700/1‑101 (MQ=255)
        ggATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCggg                                                                                              >  1:396042/1‑101 (MQ=255)
         gATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGCTCATGGCCTACGTCGTGc                                                                                             >  1:1841863/1‑101 (MQ=255)
          aTCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGcc                                                                                            >  2:1793991/1‑100 (MQ=255)
           tCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAg                                                                                           <  2:4391529/101‑1 (MQ=255)
             gAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTg                                                                                         <  1:2051392/101‑1 (MQ=255)
              aaaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTgg                                                                                        <  1:4048777/101‑1 (MQ=255)
              aaaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTgg                                                                                        >  1:3130263/1‑101 (MQ=255)
              aaaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCGCCGCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTgg                                                                                        >  1:793342/1‑101 (MQ=255)
               aaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGa                                                                                       <  2:3327301/101‑1 (MQ=255)
               aaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGa                                                                                       <  2:1020171/101‑1 (MQ=255)
               aaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGa                                                                                       >  1:1153853/1‑101 (MQ=255)
               aaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTAGa                                                                                       <  1:2884298/101‑1 (MQ=255)
                  tGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTg                                                                                    <  2:3488435/101‑1 (MQ=255)
                  tGCGCCGCCTAACGCATCCGGTCGAGCGCCGCGCAGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTg                                                                                    <  2:793240/101‑1 (MQ=255)
                   gcgcCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTgg                                                                                   <  2:2585300/101‑1 (MQ=255)
                    cgccgcCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGt                                                                                  <  2:3101028/101‑1 (MQ=255)
                       cgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGc                                                                               <  1:21109/101‑1 (MQ=255)
                        cccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCt                                                                              <  2:376160/100‑1 (MQ=255)
                           tATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTccc                                                                           >  2:2679476/1‑101 (MQ=255)
                           tATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTccc                                                                           <  2:3079762/101‑1 (MQ=255)
                            aTCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccc                                                                          >  1:2425144/1‑101 (MQ=255)
                             tCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTccccc                                                                         <  2:2637019/101‑1 (MQ=255)
                             tCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTccccc                                                                         >  1:3356076/1‑101 (MQ=255)
                             tCGCATCCGGTCGAGCGCCGCGCAGTCGCTACCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGGCGCTccccc                                                                         >  1:4492324/1‑101 (MQ=255)
                              cGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccccc                                                                        >  1:2317295/1‑101 (MQ=255)
                              cGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccccc                                                                        <  1:2775648/101‑1 (MQ=255)
                               gCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGGGGACTGGGCGCCccccccc                                                                       >  1:3460235/1‑100 (MQ=38)
                                    cGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGTCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCgggg                                                                  >  1:3065895/1‑101 (MQ=255)
                                      gTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCCCGGGGGc                                                                >  1:4415115/1‑101 (MQ=255)
                                        cGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCgg                                                              >  2:2854719/1‑100 (MQ=255)
                                        cGAGCGCCGCGCAGTCGCTCCCCCAGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGt                                                              <  2:2971172/101‑1 (MQ=255)
                                         gAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGGGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCCGTGGCCTGGTCGCCCCCCCCGCGGGGGCggg                                                             >  1:1194560/1‑99 (MQ=38)
                                         gAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTg                                                             >  1:1687840/1‑101 (MQ=255)
                                          aGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTgg                                                            >  2:3683644/1‑101 (MQ=255)
                                           gcgcCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGa                                                           <  1:3678662/101‑1 (MQ=255)
                                           gcgcCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGa                                                           <  1:4903467/101‑1 (MQ=255)
                                            gaacaggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGCTCGACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                          <  1:4626638/92‑1 (MQ=2)
                                             gccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                         <  2:1840048/101‑1 (MQ=255)
                                             gccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                         <  2:2784360/101‑1 (MQ=255)
                                             gccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGACCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                         >  2:4366721/1‑101 (MQ=255)
                                             gccgcGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                         <  2:4004429/101‑1 (MQ=255)
                                                gcgcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                      <  2:728739/101‑1 (MQ=255)
                                                 cgcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                     <  1:1599981/101‑1 (MQ=255)
                                                  gcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                    >  1:4806706/1‑101 (MQ=255)
                                                   cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                   <  2:168745/101‑1 (MQ=255)
                                                   aggTCGCGCCACCCGGGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                   <  2:2258586/99‑1 (MQ=255)
                                                    tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGCTCGACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgg                                                  <  1:631125/100‑3 (MQ=2)
                                                    aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTc                                                  <  2:3195784/101‑1 (MQ=255)
                                                    aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTc                                                  <  1:1916348/101‑1 (MQ=255)
                                                     gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCg                                                 >  1:712807/1‑101 (MQ=255)
                                                      tCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGc                                                <  2:2185279/101‑1 (MQ=255)
                                                       cGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGcc                                               >  1:1620817/1‑101 (MQ=255)
                                                         nnnncnnnNNGGNGGCGNGNANCGAAACGATCCGTNCATGGCCTACGTCGGTCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCt                                             <  2:2397843/97‑1 (MQ=37)
                                                         cTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCt                                             >  2:369610/1‑101 (MQ=255)
                                                         cTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCt                                             >  2:464191/1‑101 (MQ=255)
                                                         cTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCCCGGGGCCG.GGGTTGGAACCCGCCCTc                                            >  2:1350029/1‑101 (MQ=37)
                                                          tCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTGGGGCAGCGGACGGGTCGCCACCCCCGCGGGGGCGTCGATCGAAAGCCGcccac                                            >  2:5027274/1‑99 (MQ=25)
                                                          tCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTc                                            <  1:2854932/101‑1 (MQ=255)
                                                           ccccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCg                                           <  1:726125/101‑1 (MQ=255)
                                                           ccccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCg                                           >  1:4857594/1‑101 (MQ=255)
                                                                cgcgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcg                                      <  1:2183926/101‑1 (MQ=255)
                                                                 gcgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgc                                     >  1:2723576/1‑101 (MQ=255)
                                                                  cgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgcg                                    >  2:1852836/1‑101 (MQ=255)
                                                                  cgGGGGAGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgcg                                    >  2:2904946/1‑101 (MQ=255)
                                                                    ggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCt                                  >  1:4152917/1‑101 (MQ=255)
                                                                    ggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGAAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCt                                  >  2:4063597/1‑101 (MQ=255)
                                                                       gCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGc                               <  1:1013914/101‑1 (MQ=255)
                                                                          tGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                            >  1:3727840/1‑101 (MQ=255)
                                                                          tGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCCCGCTTCCCTc                            >  2:556709/1‑101 (MQ=255)
                                                                           ggATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCa                           >  1:2932487/1‑101 (MQ=255)
                                                                            gATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCAc                          <  2:4532085/101‑1 (MQ=255)
                                                                             aTCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACg                         >  2:4415350/1‑101 (MQ=255)
                                                                              tCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGc                        <  1:1041999/101‑1 (MQ=255)
                                                                                gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                      >  2:3455376/1‑101 (MQ=255)
                                                                                gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                      >  2:1564135/1‑101 (MQ=255)
                                                                                 aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTg                     >  1:1252000/1‑101 (MQ=255)
                                                                                    cGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCcg                  <  1:4381923/101‑1 (MQ=255)
                                                                                     gATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCcgc                 <  1:4312263/101‑1 (MQ=255)
                                                                                      aTGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCcgct                >  2:4951606/1‑100 (MQ=255)
                                                                                       tGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGt               <  1:228689/101‑1 (MQ=255)
                                                                                         cGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCg             <  2:4836104/101‑1 (MQ=255)
                                                                                           ttCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCt           <  2:2017453/101‑1 (MQ=255)
                                                                                            tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTc          >  2:3013636/1‑101 (MQ=255)
                                                                                             cATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcc         >  1:1375475/1‑101 (MQ=255)
                                                                                               tGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccc       <  1:1766122/101‑1 (MQ=255)
                                                                                               tGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccc       >  1:1492346/1‑101 (MQ=255)
                                                                                                ggCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTccccc      <  2:4211254/101‑1 (MQ=255)
                                                                                                 gCCTACGGCGGGCAGTGGACTGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccccc     <  1:1412277/101‑1 (MQ=255)
                                                                                                  ccTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCcg    >  2:3259117/1‑101 (MQ=255)
                                                                                                   cTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCcgc   >  2:4433299/1‑101 (MQ=255)
                                                                                                   cTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCcgc   <  1:974746/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
GGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCG  >  NC_014034/1395846‑1396046

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