Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,958 T→G 3.0% intergenic (+770/‑161) 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,9580T→G 3.0% 8.0 ‑1.9 ‑0.1239intergenic (+770/‑161)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/6):  ref base (125/105):  total (127/112)
Fisher's exact test strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.14e-01

Only 100 of 254 total aligned reads displayed.

CGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGATC  >  NC_014034/1395858‑1396058
                                                                                                    |                                                                                                    
 gAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTg                                                                                                     <  1:2331528/101‑1 (MQ=255)
 gAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTg                                                                                                     >  1:4445293/1‑101 (MQ=255)
 gAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTg                                                                                                     <  1:2051392/101‑1 (MQ=255)
  aaaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTgg                                                                                                    >  1:3130263/1‑101 (MQ=255)
  aaaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTgg                                                                                                    <  1:123155/101‑1 (MQ=255)
  aaaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTgg                                                                                                    <  2:398564/101‑1 (MQ=255)
  aaaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGGgg                                                                                                    >  2:4938969/1‑101 (MQ=255)
  aaaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCGCCGCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTgg                                                                                                    >  1:793342/1‑101 (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)
       gcgcCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTgg                                                                                               <  2:2585300/101‑1 (MQ=255)
         cccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTc                                                                                             <  2:3688301/100‑1 (MQ=255)
          ccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCg                                                                                            >  1:4976927/1‑101 (MQ=255)
            gccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCt                                                                                          >  2:3103112/1‑101 (MQ=255)
               tATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTccc                                                                                       <  2:3079762/101‑1 (MQ=255)
                aTCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccc                                                                                      <  2:2878121/101‑1 (MQ=255)
                 tCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTccccc                                                                                     <  2:2637019/101‑1 (MQ=255)
                  cGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccccc                                                                                    <  1:2775648/101‑1 (MQ=255)
                   gCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGGGGACTGGGCGCCccccccc                                                                                   >  1:3460235/1‑100 (MQ=38)
                    cATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCcgc                                                                                  <  1:976344/101‑1 (MQ=255)
                          gTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGc                                                                            >  1:1290091/1‑101 (MQ=255)
                           tCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCg                                                                           >  2:4367195/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)
                              aGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTgg                                                                        >  2:3683644/1‑101 (MQ=255)
                               gcgcCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGa                                                                       <  1:3678662/101‑1 (MQ=255)
                                cgccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                      <  2:2997020/101‑1 (MQ=255)
                                cgccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                      <  1:1794124/101‑1 (MQ=255)
                                 gccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                     <  2:2784360/101‑1 (MQ=255)
                                 gccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                     <  1:1402640/101‑1 (MQ=255)
                                 gccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                     <  1:4399796/101‑1 (MQ=255)
                                 gccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                     <  2:1840048/101‑1 (MQ=255)
                                 gccgcGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                     <  1:479992/101‑1 (MQ=255)
                                  ccgcGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                    <  1:4599193/101‑1 (MQ=255)
                                    gcgcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                  <  2:728739/101‑1 (MQ=255)
                                      gcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                <  1:2553952/101‑1 (MQ=255)
                                       atgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGCTCGACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGGGGGTCGAAAcg                                                               >  1:2967903/3‑100 (MQ=2)
                                       aggTCGCGCCACCCGGGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                               <  2:2258586/99‑1 (MQ=255)
                                        tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGCTCGACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgg                                                              <  1:631125/100‑3 (MQ=2)
                                        aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTc                                                              <  2:3195784/101‑1 (MQ=255)
                                         gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCg                                                             >  1:712807/1‑101 (MQ=255)
                                          tCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGc                                                            <  2:2185279/101‑1 (MQ=255)
                                           cGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGcc                                                           >  1:1997426/1‑101 (MQ=255)
                                           cGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGcc                                                           >  1:1620817/1‑101 (MQ=255)
                                             nnnncnnnNNGGNGGCGNGNANCGAAACGATCCGTNCATGGCCTACGTCGGTCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCt                                                         <  2:2397843/97‑1 (MQ=37)
                                             cTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCCCGGGGCCG.GGGTTGGAACCCGCCCTc                                                        >  2:1350029/1‑101 (MQ=37)
                                              tCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTGGGGCAGCGGACGGGTCGCCACCCCCGCGGGGGCGTCGATCGAAAGCCGcccac                                                        >  2:5027274/1‑99 (MQ=25)
                                              tCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTc                                                        <  1:3262248/101‑1 (MQ=255)
                                              tCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTc                                                        <  1:2854932/101‑1 (MQ=255)
                                              gccccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTc                                                        <  2:1687711/100‑1 (MQ=255)
                                              cccccccGTGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTc                                                        <  1:4363598/100‑1 (MQ=255)
                                               ccccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGCTCGAAACTCGCCCTCg                                                       >  1:837314/1‑101 (MQ=255)
                                               ccccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCg                                                       <  1:726125/101‑1 (MQ=255)
                                               ccccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCg                                                       >  1:4211571/1‑101 (MQ=255)
                                                cccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGa                                                      <  1:4855296/101‑1 (MQ=255)
                                                     gcgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgc                                                 >  1:2095565/1‑101 (MQ=255)
                                                      cgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgcg                                                >  2:1852836/1‑101 (MQ=255)
                                                         gggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCtt                                             >  2:3262068/1‑101 (MQ=255)
                                                         gggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCtt                                             <  1:2182988/101‑1 (MQ=255)
                                                           gCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGc                                           <  1:1013914/101‑1 (MQ=255)
                                                             gTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCt                                         <  2:56780/101‑1 (MQ=255)
                                                             gTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCACTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGACt                                         >  2:4328012/1‑101 (MQ=255)
                                                               ggATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCa                                       >  1:150255/1‑101 (MQ=255)
                                                                 aTCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACg                                     <  2:4126614/101‑1 (MQ=255)
                                                                  tCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGc                                    <  1:1041999/101‑1 (MQ=255)
                                                                   cGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGcc                                   >  1:191017/1‑101 (MQ=255)
                                                                    gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                                  >  2:3455376/1‑101 (MQ=255)
                                                                    gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                                  >  1:3080522/1‑101 (MQ=255)
                                                                     aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTg                                 >  2:4522897/1‑101 (MQ=255)
                                                                     aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTg                                 <  1:170377/101‑1 (MQ=255)
                                                                          aTGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCcgcg                            >  2:523347/1‑101 (MQ=255)
                                                                           tGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGt                           <  1:228689/101‑1 (MQ=255)
                                                                            gCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTc                          <  1:4367512/101‑1 (MQ=255)
                                                                            gCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTc                          <  2:1446199/101‑1 (MQ=255)
                                                                             cGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCg                         <  2:4836104/101‑1 (MQ=255)
                                                                             cGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCg                         >  1:3027118/1‑101 (MQ=255)
                                                                              gTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGc                        >  2:816159/1‑101 (MQ=255)
                                                                              gTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGc                        >  1:4878176/1‑101 (MQ=255)
                                                                               ttCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCt                       <  2:2017453/101‑1 (MQ=255)
                                                                                tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTc                      >  1:1413128/1‑101 (MQ=255)
                                                                                    ggCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTccccc                  >  1:2956853/1‑101 (MQ=255)
                                                                                    ggCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTccccc                  >  1:3504857/1‑101 (MQ=255)
                                                                                    ggCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTccccc                  <  2:4211254/101‑1 (MQ=255)
                                                                                      ccTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCcg                >  2:2674060/1‑101 (MQ=255)
                                                                                      ccTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCcg                <  1:4007528/101‑1 (MQ=255)
                                                                                       cTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCcgc               >  2:4433299/1‑101 (MQ=255)
                                                                                        tACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCcgcg              >  2:4411338/1‑101 (MQ=255)
                                                                                          cGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCggg            >  2:3292698/1‑101 (MQ=255)
                                                                                           gTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCgggg           >  1:2618753/1‑101 (MQ=255)
                                                                                           gTCGGGCAGTGGACTGGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGCCTTCGATCCGCACTGGCCGCACGCATGCCGCGGCGGTCCCCCCGCgggg           >  2:2755143/1‑101 (MQ=25)
                                                                                             cGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGGGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGc         <  1:4501849/101‑1 (MQ=255)
                                                                                             cGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGc         >  2:3828170/1‑101 (MQ=255)
                                                                                                gCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTg      <  1:38598/101‑1 (MQ=255)
                                                                                                gCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTg      >  1:1632106/1‑101 (MQ=255)
                                                                                                  aGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGa    >  1:363901/1‑101 (MQ=255)
                                                                                                    tGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGATc  <  1:1566101/101‑1 (MQ=255)
                                                                                                    tGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGATc  <  2:1678815/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
CGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGGATC  >  NC_014034/1395858‑1396058

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