Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,670 C→G 9.1% intergenic (+482/‑449) 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,6700C→G 9.1% 79.0 ‑2.2 ‑1.2392intergenic (+482/‑449)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (14/22):  ref base (189/166):  total (203/189)
Fisher's exact test strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.17e-02

Only 100 of 448 total aligned reads displayed.

GCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTC  >  NC_014034/1395570‑1395770
                                                                                                    |                                                                                                        
   ccccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATg                                                                                                       <  2:7257164/101‑1 (MQ=255)
   ccccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATg                                                                                                       <  2:2712063/101‑1 (MQ=255)
   ccccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATg                                                                                                       >  2:3243249/1‑101 (MQ=255)
     ccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTa                                                                                                     >  1:6547022/1‑101 (MQ=255)
     ccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTa                                                                                                     <  2:1988307/101‑1 (MQ=255)
     ccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTa                                                                                                     <  2:3950663/101‑1 (MQ=255)
      cccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGACACAGCTGCATGTAg                                                                                                    >  1:7558207/1‑101 (MQ=255)
       ccGCGGGGGCGTGGATAGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGt                                                                                                   <  1:4619513/101‑1 (MQ=255)
                 gTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGTCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAg                                                                                      >  1:7270667/1‑101 (MQ=255)
                 gTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAg                                                                                      >  1:1982673/1‑101 (MQ=255)
                 gTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAg                                                                                      >  2:8709075/1‑101 (MQ=255)
                  tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGt                                                                                     >  2:746256/1‑101 (MQ=255)
                   ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGtt                                                                                    <  2:2848984/101‑1 (MQ=255)
                   ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGtt                                                                                    >  2:6799099/1‑101 (MQ=255)
                   ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGtt                                                                                    >  1:8490476/1‑101 (MQ=255)
                   ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGtt                                                                                    >  1:8316571/1‑101 (MQ=255)
                   ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGtt                                                                                    >  1:7814055/1‑101 (MQ=255)
                   ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGtt                                                                                    >  2:7731172/1‑101 (MQ=255)
                    gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTc                                                                                   <  2:1858858/101‑1 (MQ=255)
                    gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTc                                                                                   <  1:5260641/101‑1 (MQ=255)
                     aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCg                                                                                  <  2:1362937/101‑1 (MQ=255)
                     aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCg                                                                                  <  1:1006223/101‑1 (MQ=255)
                       cGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.g                                                                               <  2:6036144/101‑1 (MQ=255)
                       cGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.g                                                                               >  1:6600127/1‑101 (MQ=255)
                        gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.Gc                                                                              >  1:2850418/1‑101 (MQ=255)
                         aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.Gcc                                                                             >  1:4824770/1‑101 (MQ=255)
                         aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.Gcc                                                                             >  2:6296632/1‑101 (MQ=255)
                               ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTc                                                                       <  2:7637537/101‑1 (MQ=255)
                                   gCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTc                                                                   <  2:5373687/101‑1 (MQ=255)
                                   gCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTc                                                                   >  2:9046997/1‑101 (MQ=255)
                                     gtaTTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTccc                                                                 <  2:1248803/99‑1 (MQ=255)
                                     ggATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..C.GCAGTTCGA.GCCGCAGTCGCTccc                                                                 <  2:4373872/101‑1 (MQ=255)
                                       gttgcgccTTGCGGGGGGCGAAGCGGCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCTGCATGTAGGCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTccccc                                                               <  1:4263883/94‑1 (MQ=37)
                                       aTTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTccccc                                                               >  1:1659615/1‑101 (MQ=255)
                                        ttGTCCCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTcccccc                                                              <  1:3561818/101‑1 (MQ=255)
                                        tggTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTcccccc                                                              <  1:155663/99‑1 (MQ=255)
                                         tgccccTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                             <  2:819148/96‑1 (MQ=255)
                                          gTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                            <  2:2113328/101‑1 (MQ=255)
                                          gTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTAGCTCCCCCcgc                                                            >  1:555498/1‑101 (MQ=255)
                                          gTCACTTTCGCGGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                            <  2:2437878/101‑1 (MQ=255)
                                            cACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgg                                                          >  1:980219/1‑101 (MQ=255)
                                            cACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgg                                                          >  1:4450982/1‑101 (MQ=255)
                                               tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                       >  2:8475732/1‑101 (MQ=255)
                                                ttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.CCCGCTGGA.GCCGCAGGCGCTTCCCCCGCGGCGGc                                                      >  2:7347798/1‑101 (MQ=25)
                                                   gcgGGGCGAAGCTTCGCTCCCCCCGCGGGGTCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTg                                                   <  1:9047419/101‑1 (MQ=255)
                                                   gcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTg                                                   >  1:3973849/1‑101 (MQ=255)
                                                    cgGGGCGAGGCCTCGCCCCCCCCGCGGGGGCGGGGATCGAAACATCTGCATGTGGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                  <  2:5073187/101‑1 (MQ=38)
                                                    cgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                  <  2:566506/101‑1 (MQ=255)
                                                     gtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                  <  2:4947379/93‑1 (MQ=2)
                                                      gggCGAGGCGTCGCCCCCCCCGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                <  1:1140279/101‑1 (MQ=255)
                                                      gggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                <  2:8503168/101‑1 (MQ=255)
                                                      gggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                >  1:3775268/1‑101 (MQ=255)
                                                       ggCCAACCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                               <  2:5281076/101‑1 (MQ=255)
                                                        cggttgCGTCGCTCCCCCCGCGGGGGCGTGGAGCGAAACATGGCTATGGCGGCAGTGTCA.GGTTTTCCA.GCGGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                               <  2:7054637/96‑1 (MQ=0)
                                                        ccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                               <  1:798771/96‑1 (MQ=2)
                                                        ccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                               <  2:1036298/96‑1 (MQ=2)
                                                         cgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCAGCGGGGGCGTGGATCg                                              <  1:5112703/97‑1 (MQ=2)
                                                          gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                            >  1:8604294/1‑101 (MQ=255)
                                                          gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                            >  2:1153490/1‑101 (MQ=255)
                                                           ctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                            <  2:4601799/99‑1 (MQ=2)
                                                             gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGGTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAACCt                                         >  1:6323400/1‑101 (MQ=255)
                                                             gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                          <  1:2371143/101‑1 (MQ=2)
                                                             gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                          <  2:2355175/101‑1 (MQ=2)
                                                              cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTCCGGAGGTAGTGGCCG..G.GCAGTTCGA.GCGGGAGTGGCGCCCCCCGCGGGGGCGTGGGGCGAAACTg                                        >  2:4252254/1‑101 (MQ=18)
                                                              cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                        <  1:7435450/101‑1 (MQ=255)
                                                                    cccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                  <  1:295057/100‑1 (MQ=255)
                                                                     ccccccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGa                                 <  1:6704232/101‑1 (MQ=255)
                                                                      cccccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGAt                                <  1:2657988/101‑1 (MQ=255)
                                                                      cccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGAt                                <  2:6333356/101‑1 (MQ=255)
                                                                       ccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATc                               <  1:6187466/101‑1 (MQ=255)
                                                                         ccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTg                             <  2:6151622/101‑1 (MQ=255)
                                                                            ggggggCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcg                          <  1:3123282/100‑1 (MQ=255)
                                                                                ggCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAGTGGTATGGATCTGGCGCaaa                      >  1:8807070/1‑101 (MQ=255)
                                                                                 gCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                     <  1:2084392/101‑1 (MQ=255)
                                                                                 gCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                     <  2:5175828/101‑1 (MQ=255)
                                                                                  cGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACg                    >  1:8518422/1‑101 (MQ=255)
                                                                                  cGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACg                    <  2:3497480/101‑1 (MQ=255)
                                                                                    tGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                  >  2:6671083/1‑101 (MQ=255)
                                                                                      gATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa                >  1:227193/1‑101 (MQ=255)
                                                                                      gATCGAAACATCTGCATGTAGCGGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa                <  2:6546762/101‑1 (MQ=255)
                                                                                      gATCGAAACACGTGCAGAGGGTGGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa                <  2:8854089/101‑1 (MQ=25)
                                                                                       aTCGAAACATCTGCATGTGGGCGGCG..G.GCAGTGCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAg               <  2:8600911/101‑1 (MQ=255)
                                                                                          gAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt            <  1:6799370/101‑1 (MQ=255)
                                                                                          gAAACATCTGCATGTAGTCGGCG..G.GCAGTCCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt            <  2:5582109/101‑1 (MQ=255)
                                                                                           aaaCGTCTGCATGTGGCCGGCG..G.GCAGTTGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg           <  2:3019391/101‑1 (MQ=37)
                                                                                           aaaCATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg           >  1:2797769/1‑101 (MQ=255)
                                                                                           aaaCATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg           >  2:7539837/1‑101 (MQ=255)
                                                                                           aaaCATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg           >  1:1380775/1‑101 (MQ=255)
                                                                                           aaaCATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg           >  1:1861645/1‑101 (MQ=255)
                                                                                            aaCATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt          >  2:7996551/1‑101 (MQ=255)
                                                                                            aaCATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt          <  1:5654221/101‑1 (MQ=255)
                                                                                             aCATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtc         >  2:2634241/1‑101 (MQ=255)
                                                                                              cATCTGCATGTAGTGGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcg        <  1:3718546/101‑1 (MQ=255)
                                                                                              cATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcg        <  1:7125536/101‑1 (MQ=255)
                                                                                              cATCTGCATGCAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGCATGGCTCTGGCGCAAACGGAAAGTCGgtcg        >  2:1631007/1‑101 (MQ=37)
                                                                                               ctCTGCATGTGGCCGGCG..G.GCGGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgt       <  2:2240799/100‑1 (MQ=255)
                                                                                               aTCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACTGGTATGGGTCTGGCGCAAACGGAAAGTCTggcgg       >  2:3725528/1‑100 (MQ=255)
                                                                                                tCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc      >  2:1499304/1‑101 (MQ=255)
                                                                                                   gCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCt   >  1:5940120/1‑101 (MQ=255)
                                                                                                    cATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTc  >  2:5752516/1‑101 (MQ=255)
                                                                                                    |                                                                                                        
GCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTC  >  NC_014034/1395570‑1395770

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