Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,698 Δ1 bp 5.3% intergenic (+510/‑421) 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,6980A→. 5.3% 29.9 ‑1.4 ‑0.0245intergenic (+510/‑421)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (5/8):  ref base (134/98):  total (139/106)
Fisher's exact test strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.89e-01

Only 100 of 262 total aligned reads displayed.

CCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTG  >  NC_014034/1395598‑1395798
                                                                                                        |                                                                                                       
 cgccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAg                                                                                                        <  1:3229415/99‑1 (MQ=255)
 caccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAg                                                                                                        <  1:4240796/101‑1 (MQ=255)
 caccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAg                                                                                                        >  1:2184386/1‑101 (MQ=255)
   ccTTGCGGATTGTCCCCTTCGCGGGGCGAAGCGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTc                                                                                                      <  1:1424475/101‑1 (MQ=255)
    cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCg                                                                                                     >  2:2829198/1‑101 (MQ=255)
     ttGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGc                                                                                                    >  2:1897122/1‑101 (MQ=255)
      tggggATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCt                                                                                                   <  2:3863349/98‑1 (MQ=255)
      tGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCt                                                                                                   <  1:4413393/101‑1 (MQ=255)
       gCGGATTGTCCCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTc                                                                                                  <  1:93345/101‑1 (MQ=255)
          gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTcccc                                                                                               >  2:465538/1‑101 (MQ=255)
          gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTcccc                                                                                               >  2:4133052/1‑101 (MQ=255)
          gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTcccc                                                                                               >  1:4092808/1‑101 (MQ=255)
           aTTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTccccc                                                                                              <  1:1519990/101‑1 (MQ=255)
           aTTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTccccc                                                                                              >  1:190545/1‑101 (MQ=255)
              gTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                                                           <  2:2048398/101‑1 (MQ=255)
                 ccTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggg                                                                                        <  1:1624763/100‑1 (MQ=255)
                  cTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgggg                                                                                       >  2:2725037/1‑101 (MQ=255)
                   tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                                                      >  2:3945711/1‑101 (MQ=255)
                      cgcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..T.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGt                                                                                   <  2:1151018/101‑1 (MQ=255)
                       gcgGGGCGAAGCGTCCCCCCCCCCGCGGGGGCGTGGATCCAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTg                                                                                  <  2:2365053/101‑1 (MQ=255)
                        cgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGCgg                                                                                 >  1:102714/1‑101 (MQ=255)
                         gtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                                                 <  2:370638/93‑1 (MQ=2)
                          tgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                                >  1:901269/8‑101 (MQ=0)
                           ggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                              >  1:3742416/1‑101 (MQ=255)
                              gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           >  1:1311469/1‑101 (MQ=255)
                              gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           <  1:1719810/101‑1 (MQ=255)
                              gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           >  1:525475/1‑101 (MQ=255)
                               ctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           <  1:4388192/99‑1 (MQ=2)
                               ctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           <  1:4626361/99‑1 (MQ=2)
                               ctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           <  2:1510206/99‑1 (MQ=2)
                               aaGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                          >  1:3595838/1‑101 (MQ=255)
                                 gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCGGCAGTGTCA.GGTTTTCCA.GCGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                         >  2:171963/1‑101 (MQ=2)
                                 gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                         <  2:1202685/101‑1 (MQ=2)
                                  cgccGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                      <  1:3299792/98‑1 (MQ=255)
                                  cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                      >  2:489493/1‑101 (MQ=255)
                                  cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                      >  1:1269588/1‑101 (MQ=255)
                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGt                                                                    >  2:1572560/1‑101 (MQ=255)
                                      gccccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTAt                                                                  <  1:3509010/98‑1 (MQ=255)
                                       ccccccccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATg                                                                 <  2:2462964/99‑1 (MQ=255)
                                        tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.gggctgg                                                                >  2:3718231/1‑96 (MQ=255)
                                        cccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                <  2:92621/100‑1 (MQ=255)
                                         ccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGa                                                               >  2:2570285/1‑101 (MQ=255)
                                           ccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATc                                                             <  1:2749405/101‑1 (MQ=255)
                                             ccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTg                                                           >  1:2046872/1‑101 (MQ=255)
                                              cgggggggGGTGGATCGAAACATCTGCATGTAGTCGGGG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTgg                                                          <  2:3200289/98‑1 (MQ=255)
                                                     gCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAAc                                                   >  1:4164529/1‑101 (MQ=255)
                                                     gCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAAc                                                   >  1:2491245/1‑101 (MQ=255)
                                                      cGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACg                                                  >  1:2124910/1‑101 (MQ=255)
                                                       gTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACgg                                                 >  2:4544386/1‑101 (MQ=255)
                                                       gTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACgg                                                 >  2:1390499/1‑101 (MQ=255)
                                                        tGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..a                                              >  2:3619179/1‑101 (MQ=255)
                                                         ggATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..aa                                             >  2:4207078/1‑101 (MQ=255)
                                                          gATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..aaa                                            >  2:3274515/1‑101 (MQ=255)
                                                          gATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..aaa                                            >  1:196784/1‑101 (MQ=255)
                                                          gATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCCGGCGCAAACGG..caa                                            >  2:3115960/1‑101 (MQ=255)
                                                           aTCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAg                                           <  1:3653671/101‑1 (MQ=255)
                                                            tCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGt                                          >  1:259671/1‑101 (MQ=255)
                                                            tCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGt                                          >  2:2665598/1‑101 (MQ=255)
                                                              gAAACATCTGCATGTAGTCGGGG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCt                                        <  1:2078454/101‑1 (MQ=255)
                                                              gAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCt                                        <  1:1056151/101‑1 (MQ=255)
                                                              gAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCt                                        >  1:1759182/1‑101 (MQ=255)
                                                               aaaCATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTg                                       >  2:3580190/1‑101 (MQ=255)
                                                                aaCATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgt                                      <  1:3855499/101‑1 (MQ=255)
                                                                  cATCTGCATGTGGGCGGGG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcg                                    <  1:397770/101‑1 (MQ=255)
                                                                  cATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTcagatc                                    >  1:2410716/1‑96 (MQ=255)
                                                                  cATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcg                                    >  2:4617126/1‑101 (MQ=255)
                                                                  cATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcg                                    >  1:1927152/1‑101 (MQ=255)
                                                                   gtCTGCATGTAGGCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcgt                                   <  1:240478/100‑1 (MQ=255)
                                                                   aTCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcgt                                   <  2:2213704/101‑1 (MQ=255)
                                                                     cTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGtcgtcg                                 <  2:4412411/101‑1 (MQ=255)
                                                                     cTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGtcgtcg                                 >  1:3119137/1‑101 (MQ=255)
                                                                      tGCATGTAGTCGGGG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGc                                <  2:4596815/101‑1 (MQ=255)
                                                                       gCATGTAGTCGGGG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCCCt                               <  1:2665783/101‑1 (MQ=255)
                                                                        cATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTc                              >  1:3003576/1‑101 (MQ=255)
                                                                           gTAGTCGTCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTcccc                           >  2:3977951/1‑101 (MQ=255)
                                                                            tAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTccccc                          >  1:4517277/1‑101 (MQ=255)
                                                                            tAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTccccc                          >  1:4301493/1‑101 (MQ=255)
                                                                            tAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTccccc                          >  2:1697106/1‑101 (MQ=255)
                                                                            gaGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTccccc                          <  2:541183/100‑1 (MQ=255)
                                                                                  gcg..g.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GATATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCggg                    >  2:2951661/1‑101 (MQ=255)
                                                                                   gg..g.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCgggg                   <  1:1489895/100‑1 (MQ=255)
                                                                                       g.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATGTGGCGCAAACGG..AAAGTCTGTCCTCGCTCCCCCCGCGGGGGc                 >  2:3451924/1‑101 (MQ=255)
                                                                                       g.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGc                 >  2:4349142/1‑101 (MQ=255)
                                                                                       g.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGc                 <  1:3615158/101‑1 (MQ=255)
                                                                                        gggagtgggA.GCCGCGGTCGCTCCCCCCGGGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAAAGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGc                 <  1:3166716/93‑1 (MQ=37)
                                                                                         gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCg                >  2:1110743/1‑101 (MQ=255)
                                                                                          cAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCgg               >  2:2036165/1‑100 (MQ=255)
                                                                                             ttCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGa            >  2:686877/1‑101 (MQ=255)
                                                                                               cGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATc          <  1:479351/101‑1 (MQ=255)
                                                                                                gA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCg         <  2:102700/101‑1 (MQ=255)
                                                                                                caagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTgg             <  1:2860000/94‑1 (MQ=2)
                                                                                                 a.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGa        >  2:86100/1‑101 (MQ=255)
                                                                                                    tcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa      <  2:3542269/100‑1 (MQ=255)
                                                                                                    ccGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa      <  2:3475789/101‑1 (MQ=255)
                                                                                                    ccGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGGCGCTCCCCCCGCGGGGGCGTGGATCGaaa      <  2:1446838/101‑1 (MQ=255)
                                                                                                    attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCg         <  1:3400467/98‑1 (MQ=2)
                                                                                                     ttcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGa        >  1:2222085/3‑101 (MQ=2)
                                                                                                      gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc    <  1:4137029/101‑1 (MQ=255)
                                                                                                       cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa      <  1:2800545/101‑1 (MQ=2)
                                                                                                       cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATAGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCCCGGGGGCGTGGATCGgaa      >  2:760971/1‑101 (MQ=2)
                                                                                                        |                                                                                                       
CCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTG  >  NC_014034/1395598‑1395798

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