Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,943 G→A 3.5% intergenic (+755/‑176) 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,9430G→A 3.5% 8.4 ‑2.1 ‑0.0200intergenic (+755/‑176)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (1/6):  ref base (98/94):  total (99/101)
Fisher's exact test strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.91e-01

Only 100 of 212 total aligned reads displayed.

CGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCC  >  NC_014034/1395843‑1396043
                                                                                                    |                                                                                                     
cgggggggCGTGGGTCGGAGCTGGGCCGCCTATCGCATCCGGTCGAGCGCCGGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATgg                                                                                                       <  1:2965116/98‑1 (MQ=37)
cgcgGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATgg                                                                                                       >  2:2328651/1‑101 (MQ=255)
  cgGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGcc                                                                                                     >  1:1814253/1‑101 (MQ=255)
   gggggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCt                                                                                                    >  1:2701549/1‑101 (MQ=255)
   gggggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCt                                                                                                    >  1:404846/1‑101 (MQ=255)
   gggggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCt                                                                                                    >  1:1461662/1‑101 (MQ=255)
   gggggCGTGGATCGAAACTGCGCCGCCTAGCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCt                                                                                                    <  1:2409565/101‑1 (MQ=255)
     gggAGTGGATCGAACCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTAc                                                                                                  <  1:324238/101‑1 (MQ=255)
      ggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACg                                                                                                 >  1:50225/1‑101 (MQ=255)
      ggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACg                                                                                                 <  2:681734/101‑1 (MQ=255)
       gCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGt                                                                                                >  1:86969/1‑101 (MQ=255)
        cGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTc                                                                                               <  2:175458/101‑1 (MQ=255)
         gTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGGGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGGTCATGGCCTACGTCg                                                                                              >  1:851861/1‑101 (MQ=255)
           ggATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCggg                                                                                            >  2:4646752/1‑101 (MQ=255)
            gATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGc                                                                                           <  2:3952884/101‑1 (MQ=255)
            gATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGc                                                                                           >  2:2655055/1‑101 (MQ=255)
             aTCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCa                                                                                          <  2:3636679/101‑1 (MQ=255)
               cGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAgg                                                                                        >  2:4015885/1‑100 (MQ=255)
               cGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAgg                                                                                        >  2:2156573/1‑100 (MQ=255)
               cGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGt                                                                                        >  2:2763740/1‑101 (MQ=255)
                gAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTg                                                                                       >  2:4170018/1‑101 (MQ=255)
                gAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTg                                                                                       >  1:3935720/1‑101 (MQ=255)
                  aaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGa                                                                                     <  2:273828/101‑1 (MQ=255)
                  aaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGa                                                                                     >  1:793757/1‑101 (MQ=255)
                    cTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACt                                                                                   <  1:1562980/101‑1 (MQ=255)
                    cTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACt                                                                                   <  2:3296748/101‑1 (MQ=255)
                     tGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTg                                                                                  <  1:2499404/101‑1 (MQ=255)
                      gcgcCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTgg                                                                                 <  2:1590081/101‑1 (MQ=255)
                       cgccgcCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGt                                                                                <  2:526242/101‑1 (MQ=255)
                        gccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTc                                                                               <  1:4556401/101‑1 (MQ=255)
                        gccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTc                                                                               <  2:1291395/101‑1 (MQ=255)
                         ccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCg                                                                              <  2:3068396/101‑1 (MQ=255)
                            ccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTc                                                                           <  1:420453/101‑1 (MQ=255)
                            ccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTc                                                                           >  2:829270/1‑101 (MQ=255)
                             cTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcc                                                                          >  1:1173263/1‑101 (MQ=255)
                              tATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTccc                                                                         >  1:3409275/1‑101 (MQ=255)
                              tATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTccc                                                                         >  2:1363095/1‑101 (MQ=255)
                                 cGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccccc                                                                      >  2:2567574/1‑101 (MQ=255)
                                  gCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCcg                                                                     <  1:2867270/101‑1 (MQ=255)
                                    atacGGACGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCATGCGCTCACAGCCTACGTCGGGCAGCGGACTGGCCGGTCCCCCcgtg                                                                   >  2:3918164/4‑99 (MQ=21)
                                      ccGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggg                                                                 <  1:3748608/101‑1 (MQ=255)
                                      ccGGTCGAGCCCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggg                                                                 <  1:2048461/101‑1 (MQ=255)
                                       cGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCgggg                                                                <  1:2567726/101‑1 (MQ=255)
                                        ggTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggggg                                                               >  1:4570409/1‑101 (MQ=255)
                                        ggTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggggg                                                               >  1:3066500/1‑101 (MQ=255)
                                         gTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGTGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGc                                                              <  1:2324274/101‑1 (MQ=255)
                                         gTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGc                                                              <  2:627093/101‑1 (MQ=255)
                                           cGAGCGCCGCGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.Gt                                                           <  1:3272763/101‑1 (MQ=255)
                                              gcgcCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGa                                                        >  1:1191625/1‑101 (MQ=255)
                                              gcgcCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGa                                                        <  1:2998401/101‑1 (MQ=255)
                                               gaacaggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGCTCGACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGC.GTGGAt                                                       <  1:3163728/92‑1 (MQ=2)
                                                gccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATc                                                      <  1:2092797/101‑1 (MQ=255)
                                                gccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATc                                                      <  1:1259858/101‑1 (MQ=255)
                                                  caggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGCTCGACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGC.GTGGATCGa                                                    <  2:1680610/95‑1 (MQ=2)
                                                   gcgcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGaa                                                   <  1:2681727/101‑1 (MQ=255)
                                                   aggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGCTCGACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGC.GTGGATCGaa                                                   <  2:4311944/96‑1 (MQ=2)
                                                   aggatgTCGCGCCCCCCGCGGGGGCGTAGATCGAAACTATGCTGGCTCGACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGC.GTGGATCGaa                                                   <  1:1994921/96‑1 (MQ=0)
                                                    cgcAGTCCCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGaaa                                                  <  2:3354343/101‑1 (MQ=255)
                                                       aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGGGGGATCGaaccc                                                >  1:686447/1‑98 (MQ=37)
                                                       aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GGGGGTTGAAACtg                                               >  2:2284172/1‑100 (MQ=255)
                                                        gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCg                                              >  1:4381871/1‑101 (MQ=255)
                                                          cGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGcc                                            >  2:107422/1‑101 (MQ=255)
                                                            gtCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCt                                          <  2:1113036/100‑1 (MQ=255)
                                                             tCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTc                                         <  2:86953/101‑1 (MQ=255)
                                                              ccccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCg                                        <  2:1461568/101‑1 (MQ=255)
                                                              ccccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCg                                        <  1:1752110/101‑1 (MQ=255)
                                                              ccccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGATCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCg                                        >  1:3548071/1‑101 (MQ=255)
                                                                ccccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGAt                                      >  2:758937/1‑101 (MQ=255)
                                                                 cccGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATc                                     <  2:1884033/101‑1 (MQ=255)
                                                                   cgcgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCcg                                   >  2:388903/1‑101 (MQ=255)
                                                                    gcgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCcgc                                  <  1:953119/101‑1 (MQ=255)
                                                                     cgGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCcgcg                                 <  1:729362/101‑1 (MQ=255)
                                                                       ggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCt                               <  1:2823714/101‑1 (MQ=255)
                                                                       ggggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCt                               <  2:2900682/101‑1 (MQ=255)
                                                                         ggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTg                             >  1:4470961/1‑101 (MQ=255)
                                                                          gCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTAGATCCGCGCTTGc                            >  2:4122472/1‑101 (MQ=255)
                                                                           cGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGcc                           >  1:2173617/1‑101 (MQ=255)
                                                                            gTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCt                          >  2:1813916/1‑101 (MQ=255)
                                                                             tGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                         >  2:3736250/1‑101 (MQ=255)
                                                                              ggATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCa                        <  2:2987414/101‑1 (MQ=255)
                                                                              ggATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCa                        >  2:2510229/1‑101 (MQ=255)
                                                                               gATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTTGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCAc                       <  2:119661/101‑1 (MQ=255)
                                                                                 tCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTccctc                     >  1:1084115/1‑97 (MQ=255)
                                                                                 tCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGc                     >  1:3537587/1‑101 (MQ=255)
                                                                                  cGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGcc                    >  2:1004476/1‑101 (MQ=255)
                                                                                   gAAACGATGCGTTCATGGGCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                   <  1:4731998/101‑1 (MQ=255)
                                                                                   gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                   >  1:2890823/1‑101 (MQ=255)
                                                                                    aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTg                  >  2:3049513/1‑101 (MQ=255)
                                                                                      aCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGcc                <  2:3143583/101‑1 (MQ=255)
                                                                                             gTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGc         >  2:3216156/1‑101 (MQ=255)
                                                                                              ttCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGCTCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCt        <  2:286958/101‑1 (MQ=255)
                                                                                               tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTc       >  1:100765/1‑101 (MQ=255)
                                                                                               tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTc       >  1:3346228/1‑101 (MQ=255)
                                                                                               tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTc       >  1:3875343/1‑101 (MQ=255)
                                                                                               tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTc       >  1:614023/1‑101 (MQ=255)
                                                                                                 aTGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTccc     >  2:40314/1‑101 (MQ=255)
                                                                                                  tGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccc    >  1:4595307/1‑101 (MQ=255)
                                                                                                  tGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccc    >  2:2472438/1‑101 (MQ=255)
                                                                                                   ggCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTccccc   >  1:832245/1‑101 (MQ=255)
                                                                                                    gCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccccc  <  2:470829/101‑1 (MQ=255)
                                                                                                    |                                                                                                     
CGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGC.GTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCC  >  NC_014034/1395843‑1396043

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