Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,946 T→G 3.0% 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 3.0% 6.0 ‑0.8 ‑0.0195intergenic (+758/‑173)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/4):  ref base (99/90):  total (101/94)
Fisher's exact test strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.88e-01

Only 100 of 204 total aligned reads displayed.

GGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCG  >  NC_014034/1395846‑1396046
                                                                                                    |                                                                                                    
gggggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCt                                                                                                      >  2:951080/1‑101 (MQ=255)
 ggggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTa                                                                                                     <  2:1570375/101‑1 (MQ=255)
  gggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTAc                                                                                                    >  2:477644/1‑101 (MQ=255)
  gggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTAc                                                                                                    >  1:2032323/1‑101 (MQ=255)
   ggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACg                                                                                                   <  2:3215140/101‑1 (MQ=255)
   ggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACg                                                                                                   <  2:3899934/101‑1 (MQ=255)
     cGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTc                                                                                                 <  2:3450195/101‑1 (MQ=255)
      gTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCg                                                                                                <  2:3000803/101‑1 (MQ=255)
       tGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGCAACGATGCGTTCATGGCCTACGTccg                                                                                               >  2:1826148/1‑99 (MQ=255)
          aTCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCa                                                                                            >  2:3757812/1‑101 (MQ=255)
          aTCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCa                                                                                            >  1:2235464/1‑101 (MQ=255)
            cGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGt                                                                                          <  1:1239611/101‑1 (MQ=255)
            cGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGt                                                                                          >  1:3821589/1‑101 (MQ=255)
            cGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGACTACGTCGGGCAGt                                                                                          >  1:4206819/1‑101 (MQ=255)
              aaaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTgg                                                                                        <  2:56106/101‑1 (MQ=255)
               aaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGa                                                                                       <  1:2453567/101‑1 (MQ=255)
                aCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGAc                                                                                      >  2:4173931/1‑101 (MQ=255)
                 cTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACt                                                                                     >  2:3273359/1‑101 (MQ=255)
                 cTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACt                                                                                     <  1:2396807/101‑1 (MQ=255)
                  tGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTg                                                                                    <  2:1239899/101‑1 (MQ=255)
                  tGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTg                                                                                    >  1:3112920/1‑101 (MQ=255)
                  tGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTg                                                                                    <  1:2615145/101‑1 (MQ=255)
                     gccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTc                                                                                 <  1:4177769/101‑1 (MQ=255)
                      ccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCg                                                                                >  1:1325941/1‑101 (MQ=255)
                          cTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcc                                                                            >  1:4121704/1‑101 (MQ=255)
                           tATCGCATCCGGTCGAGCGCCGGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTccc                                                                           <  1:1826268/101‑1 (MQ=255)
                            aTCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccc                                                                          >  1:2119717/1‑101 (MQ=255)
                            aTCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccc                                                                          <  2:4081659/101‑1 (MQ=255)
                              cGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccccc                                                                        >  1:2761873/1‑101 (MQ=255)
                               gCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCcg                                                                       <  1:493445/101‑1 (MQ=255)
                                 aTCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCcgcg                                                                     >  1:648428/1‑101 (MQ=255)
                                 aTCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCcgcg                                                                     >  2:236290/1‑101 (MQ=255)
                                     ggTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggggg                                                                 <  1:3517226/101‑1 (MQ=255)
                                      gTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGc                                                                >  2:3925998/1‑101 (MQ=255)
                                        cGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCgg                                                              >  1:3584875/1‑100 (MQ=255)
                                        cGAGCGCCACGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCgg                                                              >  2:2715490/1‑100 (MQ=255)
                                         gAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTg                                                             >  1:2812189/1‑101 (MQ=255)
                                          aGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTgg                                                            <  2:2032185/101‑1 (MQ=255)
                                            cgccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                          <  1:3767312/101‑1 (MQ=255)
                                            cgccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                          <  2:2725526/101‑1 (MQ=255)
                                             gccgcGCAGTCCCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                         <  1:1941142/101‑1 (MQ=255)
                                             aacaggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGCTCGACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                         <  2:1308793/93‑1 (MQ=2)
                                              ccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                        <  1:951165/101‑1 (MQ=255)
                                              ccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                        <  2:3486831/101‑1 (MQ=255)
                                                  gcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGGCGCTCCCCCCGCGGGGGCGGGGTTCGAAAc                                                    >  2:1654527/1‑101 (MQ=255)
                                                  gatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGCTCGACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                    <  2:1988191/98‑1 (MQ=2)
                                                   cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                   <  2:2570177/101‑1 (MQ=255)
                                                   cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                   <  1:1219046/101‑1 (MQ=255)
                                                   cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                   <  1:429437/101‑1 (MQ=255)
                                                    aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTc                                                  >  2:3671996/1‑101 (MQ=255)
                                                    aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTc                                                  >  1:453445/1‑101 (MQ=255)
                                                       cGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGCTCGACCGAGGTCACGGCGGGCTATCGTCGCCCCCCCCCCGGGGGCGGGGATCGAAAcagggc                                               >  2:4078251/1‑96 (MQ=2)
                                                       cGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGcc                                               >  1:1988632/1‑101 (MQ=255)
                                                       cGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGcc                                               >  1:587820/1‑101 (MQ=255)
                                                        gCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGccc                                              >  1:2770201/1‑101 (MQ=255)
                                                        gCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGccc                                              <  1:2726090/101‑1 (MQ=255)
                                                        gCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGccc                                              <  2:1717365/101‑1 (MQ=255)
                                                        gCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGccc                                              >  2:427593/1‑101 (MQ=255)
                                                         cTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCt                                             <  1:3273600/101‑1 (MQ=255)
                                                         cTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCt                                             <  1:152801/101‑1 (MQ=255)
                                                          tCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTc                                            <  2:1787259/101‑1 (MQ=255)
                                                          cccccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCTCTc                                            <  2:3584622/100‑1 (MQ=255)
                                                              cccGGGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACTCGCCCTCGATc                                        >  2:621988/1‑101 (MQ=255)
                                                               ccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATcc                                       <  1:2193746/101‑1 (MQ=255)
                                                               ccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATcc                                       >  1:2730971/1‑101 (MQ=255)
                                                               ccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATcc                                       <  2:728786/101‑1 (MQ=255)
                                                               ccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATcc                                       >  2:2858325/1‑101 (MQ=255)
                                                               ccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATcc                                       <  2:893048/101‑1 (MQ=255)
                                                                 gcgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATcaga                                     >  2:3282204/1‑98 (MQ=255)
                                                                   gggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCgcgc                                   >  1:2972893/1‑101 (MQ=255)
                                                                     gggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCtt                                 <  1:837335/101‑1 (MQ=255)
                                                                      ggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTg                                >  1:3331826/1‑101 (MQ=255)
                                                                        cGTGGATCGAAACGACGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGcc                              <  2:1466391/101‑1 (MQ=255)
                                                                         gTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCt                             <  2:2442046/101‑1 (MQ=255)
                                                                         gTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCt                             >  2:4226018/1‑101 (MQ=255)
                                                                           ggATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCa                           <  1:4062206/101‑1 (MQ=255)
                                                                            gATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCAc                          <  1:69424/101‑1 (MQ=255)
                                                                            gATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCAc                          <  2:2761654/101‑1 (MQ=255)
                                                                             aTCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACg                         >  1:2835380/1‑101 (MQ=255)
                                                                             aTCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACg                         >  2:2952092/1‑101 (MQ=255)
                                                                             aTCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCCCGGGGCGGTGGTTCGAAACTGCCCCCCAATCCCCCCTTCCCTCACg                         >  2:1020354/1‑101 (MQ=255)
                                                                                gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCTGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACCCCt                      >  1:3901147/1‑101 (MQ=255)
                                                                                gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                      >  1:801909/1‑101 (MQ=255)
                                                                                gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATAGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                      >  1:3505229/1‑101 (MQ=255)
                                                                                 aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCGCCAATCCGCGCTTGTCTCAAGcccc                     >  2:1503886/1‑99 (MQ=37)
                                                                                 aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTg                     >  2:2935783/1‑101 (MQ=255)
                                                                                  aaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGc                    >  1:3849187/1‑101 (MQ=255)
                                                                                   aCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGcc                   >  2:380639/1‑101 (MQ=255)
                                                                                   aCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGcc                   >  1:234124/1‑101 (MQ=255)
                                                                                     gATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGAGCTTGCCTCACGCCTGCcgc                 >  2:207810/1‑101 (MQ=255)
                                                                                      aTGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCcgcg                >  2:2241609/1‑101 (MQ=255)
                                                                                          gTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGc            >  2:49791/1‑101 (MQ=255)
                                                                                            tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTc          >  1:1669369/1‑101 (MQ=255)
                                                                                            tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTc          <  1:907308/101‑1 (MQ=255)
                                                                                            tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTc          >  1:3974919/1‑101 (MQ=255)
                                                                                             cATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcc         >  2:193331/1‑101 (MQ=255)
                                                                                             cATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcc         >  2:381471/1‑101 (MQ=255)
                                                                                                  ccTACGTC.GGCAGTGGACTGGTCCCTCCCCCCGCGGGGGCGGGGATCGAAACCCGCCCCCGATCCCCGCTTGCCCTCCCCCCGCCGCGTCCCTcccccccc   >  2:779704/1‑99 (MQ=19)
                                                                                                  ccTACGCCGGGGAGTGAACGGGCCGCGCCCCCCCCGGGGGCGTGGATCGAACCTCGCCCTCCACCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCcg    <  1:3604957/101‑1 (MQ=255)
                                                                                                   cTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCcgc   >  2:1898425/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
GGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCG  >  NC_014034/1395846‑1396046

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