Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,674 T→C 5.5% intergenic (+486/‑445) 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,6740T→C 5.5% 18.7 ‑1.0 ‑0.5217intergenic (+486/‑445)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/4):  ref base (102/101):  total (110/107)
Fisher's exact test strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.45e-01

Only 100 of 246 total aligned reads displayed.

CCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCC  >  NC_014034/1395574‑1395773
                                                                                                    |                                                                                                      
  cccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAg                                                                                                      <  2:2546423/101‑1 (MQ=255)
  cccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAg                                                                                                      <  1:2772661/101‑1 (MQ=255)
  cccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAg                                                                                                      >  1:157611/1‑101 (MQ=255)
    cgcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGGTCGAAACATCTGCATGTAGTc                                                                                                    >  1:954095/1‑101 (MQ=255)
     gcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAATCTGCATGTAGTcg                                                                                                   >  2:2137601/1‑101 (MQ=255)
      cgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcgg                                                                                                  >  1:809245/1‑101 (MQ=255)
       gggggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggc                                                                                                 >  2:32946/1‑101 (MQ=255)
        ggggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggcg                                                                                                >  1:68600/1‑101 (MQ=255)
          ggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggcg.gc                                                                                             <  1:2572949/101‑2 (MQ=255)
           gCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.Gc                                                                                           <  2:3340658/101‑1 (MQ=255)
             gTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAg                                                                                         >  2:2535801/1‑101 (MQ=255)
              tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGt                                                                                        <  2:1553762/101‑1 (MQ=255)
               ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGtt                                                                                       >  2:2352148/1‑101 (MQ=255)
                gATCGAAACCACCTTGCGGATTGGCACTTTCGCGGTGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTc                                                                                      <  2:527116/101‑1 (MQ=255)
                 aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCg                                                                                     <  2:4143959/101‑1 (MQ=255)
                 aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGAGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCg                                                                                     >  2:1146289/1‑101 (MQ=255)
                   cGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.g                                                                                  >  1:2237007/1‑101 (MQ=255)
                   cGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.g                                                                                  >  1:1802406/1‑101 (MQ=255)
                     aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.Gcc                                                                                >  2:2835884/1‑101 (MQ=255)
                     aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.Gcc                                                                                >  1:3094899/1‑101 (MQ=255)
                       accaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGc                                                                              >  2:1864242/1‑101 (MQ=255)
                        cccccTTGCGGATTGTCCCTTTCGCGGGGCGAAGCGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCa                                                                             <  1:1325651/98‑1 (MQ=255)
                        ccaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCa                                                                             <  2:194520/101‑1 (MQ=255)
                         caccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAg                                                                            >  2:621506/1‑101 (MQ=255)
                          cccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGt                                                                           <  1:4004078/100‑1 (MQ=255)
                          accTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGt                                                                           >  1:3439271/1‑101 (MQ=255)
                          accTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGt                                                                           <  1:1040768/101‑1 (MQ=255)
                           ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTc                                                                          <  1:2520057/101‑1 (MQ=255)
                              tGCGGATTGTCCCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCt                                                                       <  1:1428999/101‑1 (MQ=255)
                               gCGGATTGTCGCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTc                                                                      <  1:754805/101‑1 (MQ=255)
                               ccGGATTGTCCCTTGCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTc                                                                      <  2:4190437/100‑1 (MQ=255)
                                cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTcc                                                                     >  1:4014569/1‑101 (MQ=255)
                                  gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTcccc                                                                   <  1:2058529/101‑1 (MQ=255)
                                  gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTcccc                                                                   <  1:427003/101‑1 (MQ=255)
                                       tCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgcg                                                              >  1:3501222/1‑101 (MQ=255)
                                       tCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgcg                                                              <  1:500662/101‑1 (MQ=255)
                                        cACTTTCGCGGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgg                                                             <  1:2981603/101‑1 (MQ=255)
                                        cACTTTCGCGGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgg                                                             <  1:3878685/101‑1 (MQ=255)
                                         aCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggg                                                            >  1:1931550/1‑101 (MQ=255)
                                           tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                          >  2:712163/1‑101 (MQ=255)
                                            ttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGc                                                         <  2:3951089/101‑1 (MQ=255)
                                              cgcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCgg                                                       >  2:4280485/1‑100 (MQ=255)
                                               gcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTg                                                      <  2:2236862/101‑1 (MQ=255)
                                                gggggCGACGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                     <  1:1265649/100‑1 (MQ=255)
                                                cgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                     <  1:3678340/101‑1 (MQ=255)
                                                cgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                     <  2:2745645/101‑1 (MQ=255)
                                                cgGGGCGAAGCGTCGCCCCCCCCGCGCGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                     <  1:3735942/101‑1 (MQ=255)
                                                 gtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                     <  2:489962/93‑1 (MQ=2)
                                                  tgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGa                                                    >  2:1717133/8‑101 (MQ=0)
                                                  gggGGAGGCGCCGCGCCCCCCGCGGGGGCGTGGATCGAAACTTCTGCATCGGGTGGGGG.G.GCGGTCCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                   <  2:3506204/101‑1 (MQ=21)
                                                  gggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                   <  2:896881/101‑1 (MQ=255)
                                                    gCGAGGCGTCCCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                 <  2:1581173/101‑1 (MQ=255)
                                                    gCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                 <  1:1096445/101‑1 (MQ=255)
                                                    gCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                 <  1:3910691/101‑1 (MQ=255)
                                                    gCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                 <  1:2359705/101‑1 (MQ=255)
                                                    ccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                  <  2:3976626/96‑1 (MQ=2)
                                                     cGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                >  1:2300615/1‑101 (MQ=255)
                                                      gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCCGCTCCA.GCCGCAGTCGCTCCCCCCCCGGGGGCGGGGATCGaa                                               >  1:4089719/1‑101 (MQ=37)
                                                      gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGGGGATCGaa                                               >  1:1869936/1‑101 (MQ=255)
                                                      cgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAA.GCTGGTCCC.GTCGCCGTCGGTCCCCCCGCGGGGGCGTGGATCGa                                                <  2:154366/96‑1 (MQ=2)
                                                       aaGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                              >  2:4071815/1‑101 (MQ=255)
                                                         gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                             <  2:4005160/101‑1 (MQ=2)
                                                           gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgg                                          <  2:3788204/101‑1 (MQ=255)
                                                            gcGCCCCCCCCGGCGGGGGGTGGATCAACACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGt                                         <  1:32959/100‑1 (MQ=25)
                                                              gCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTAt                                       >  2:2973063/1‑101 (MQ=255)
                                                               ccccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATg                                      <  1:217975/99‑1 (MQ=255)
                                                                tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGttcgg                                     >  2:3422239/1‑97 (MQ=255)
                                                                 ccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGa                                    <  2:1389129/101‑1 (MQ=255)
                                                                        cgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcg                             >  1:2722511/1‑101 (MQ=255)
                                                                            ggCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaaa                         <  1:3362370/101‑1 (MQ=255)
                                                                            ggCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaaa                         <  1:621574/101‑1 (MQ=255)
                                                                            ggCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaaa                         <  2:1531056/101‑1 (MQ=255)
                                                                              cGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACg                       >  1:602825/1‑101 (MQ=255)
                                                                              cGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACg                       >  1:2661199/1‑101 (MQ=255)
                                                                               gtagATGGAAACATATGCATGTGGCCGGGG.G.GGAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                      <  1:1193797/98‑1 (MQ=25)
                                                                               gTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                      >  1:2037351/1‑101 (MQ=255)
                                                                               gTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                      >  1:3768994/1‑101 (MQ=255)
                                                                                tGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                     >  2:829063/1‑101 (MQ=255)
                                                                                tGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                     >  2:1281276/1‑101 (MQ=255)
                                                                                tGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                     <  2:1743662/101‑1 (MQ=255)
                                                                                  gATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa                   <  1:472489/101‑1 (MQ=255)
                                                                                   aTCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAg                  >  2:3603562/1‑101 (MQ=255)
                                                                                   aTCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAg                  >  1:2134865/1‑101 (MQ=255)
                                                                                    tCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGt                 <  2:321850/101‑1 (MQ=255)
                                                                                    tCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTAGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGt                 >  1:3108166/1‑101 (MQ=255)
                                                                                     cGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTAGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc                <  1:594622/101‑1 (MQ=255)
                                                                                       aaaCATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg              >  2:1691194/1‑101 (MQ=255)
                                                                                       aaaCATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg              <  2:3165424/101‑1 (MQ=255)
                                                                                          cATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcg           <  1:3632364/101‑1 (MQ=255)
                                                                                          cATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcg           >  2:4056125/1‑101 (MQ=255)
                                                                                          cATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCAGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcg           >  1:3958474/1‑101 (MQ=255)
                                                                                           agcTGCATGGAGTCGGCG.G.GGAGTTCGA.CCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgt          <  2:3044277/99‑1 (MQ=255)
                                                                                           aTCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgt          >  2:1172522/1‑101 (MQ=255)
                                                                                            tCTGGAGGTGGTGGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc         <  2:4341927/101‑1 (MQ=37)
                                                                                            tCTGCATGTGGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc         <  1:2672198/101‑1 (MQ=255)
                                                                                            tCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc         <  1:3172250/101‑1 (MQ=255)
                                                                                            tCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc         <  2:2453155/101‑1 (MQ=255)
                                                                                               tgaTGTAGGGGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGAGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCt      <  1:1000128/99‑1 (MQ=255)
                                                                                                  tGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTccc   <  1:3039365/101‑1 (MQ=255)
                                                                                                   gTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcccc  <  2:2937660/101‑1 (MQ=255)
                                                                                                    |                                                                                                      
CCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCC  >  NC_014034/1395574‑1395773

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