Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,830 +G 3.0% intergenic (+642/‑289) 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,8301.→G 3.0% 13.1 ‑1.3 ‑0.0233intergenic (+642/‑289)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/5):  ref base (121/105):  total (123/110)
Fisher's exact test strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.46e-01

Only 100 of 268 total aligned reads displayed.

CTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC  >  NC_014034/1395730‑1395930
                                                                                                       |                                                                                                    
cTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                                                                                       >  1:4036796/1‑101 (MQ=255)
cTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                                                                                       >  2:2219080/1‑101 (MQ=255)
cTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                                                                                       >  2:4311946/1‑101 (MQ=255)
 tGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.c                                                                                                     <  1:3282583/101‑1 (MQ=255)
 tGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.c                                                                                                     >  1:1873195/1‑101 (MQ=255)
 tGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCAGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.c                                                                                                     >  1:4332930/1‑101 (MQ=255)
   gTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGt                                                                                                   <  2:3621307/101‑1 (MQ=255)
   gTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGt                                                                                                   <  1:1608947/101‑1 (MQ=255)
    tATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTc                                                                                                  >  1:3443522/1‑101 (MQ=255)
      tGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGc                                                                                                >  2:4436599/1‑101 (MQ=255)
       ggATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGT.CGTCGCt                                                                                               >  2:4501529/1‑101 (MQ=255)
       ggATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCt                                                                                               >  1:3807805/1‑101 (MQ=255)
       ggATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTCATAGA.CGCCGTGGA.CGTCGCt                                                                                               >  2:3399631/1‑101 (MQ=255)
            tGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTccccc                                                                                          <  2:3862515/101‑1 (MQ=255)
            tGGCGCAAACGGAAAGTCTGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTccccc                                                                                          <  1:521250/101‑1 (MQ=255)
             ggCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                                         <  2:2233526/101‑1 (MQ=255)
             ggCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                                         <  1:4949833/101‑1 (MQ=255)
              gcgcAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCcg                                                                                        <  2:369860/101‑1 (MQ=255)
               cgcAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCcgc                                                                                       <  1:4801551/101‑1 (MQ=255)
                 cAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                                                     >  1:1500376/1‑101 (MQ=255)
                 cAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                                                     >  1:4738905/1‑101 (MQ=255)
                  aaaCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCggg                                                                                    >  1:267426/1‑101 (MQ=255)
                   aaCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgggg                                                                                   <  2:3585135/101‑1 (MQ=255)
                    aCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCggggg                                                                                  >  2:4205555/1‑101 (MQ=255)
                     cGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGc                                                                                 <  1:1008604/101‑1 (MQ=255)
                     cGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGc                                                                                 >  1:1230956/1‑101 (MQ=255)
                        ggccccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCGt                                                                               <  1:2042812/92‑1 (MQ=0)
                        ggacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCgg                                                                               >  2:4864894/10‑100 (MQ=0)
                        aaaGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTg                                                                              <  1:2752851/101‑1 (MQ=255)
                         aaGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTgg                                                                             <  1:3665525/101‑1 (MQ=255)
                         aaGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTgg                                                                             <  2:3329786/101‑1 (MQ=255)
                            tCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                          >  2:647316/1‑101 (MQ=255)
                            tCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGGGATAGA.CGCCGTGGA.CGTCGCCCCCCCCCCGGGGGCGTGGATc                                                                          >  1:1781594/1‑101 (MQ=255)
                              tGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                        <  1:2285729/101‑1 (MQ=255)
                              tGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                        >  1:2872901/1‑101 (MQ=255)
                               gtcgtcGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                       >  2:21099/1‑101 (MQ=255)
                                tcgtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                      <  2:4783071/101‑1 (MQ=255)
                                  gtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGaaccg                                                                    >  2:4059311/1‑98 (MQ=255)
                                  gtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                                    <  1:386729/101‑1 (MQ=255)
                                  gtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                                    >  2:4961735/1‑101 (MQ=255)
                                   tcgCCCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                                                   <  1:1085362/101‑1 (MQ=255)
                                    cgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTGGCTCCCCCCGCGGGGGCGTGGATCGAAACAgc                                                                  >  2:1535704/1‑101 (MQ=255)
                                    cgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgc                                                                  <  1:2986346/101‑1 (MQ=255)
                                    cgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgc                                                                  >  2:662734/1‑101 (MQ=255)
                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgcg                                                                 >  1:4391848/1‑101 (MQ=255)
                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgcg                                                                 >  1:3373364/1‑101 (MQ=255)
                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgcg                                                                 >  2:1627830/1‑101 (MQ=255)
                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgcg                                                                 <  1:984384/101‑1 (MQ=255)
                                      cTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCCCGGGGGCGTGGAtgggaccacca                                                                 >  1:1381833/1‑91 (MQ=0)
                                       tCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGcgcc                                                               <  2:2817352/101‑1 (MQ=255)
                                          ccccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCgccgcc                                                            <  1:4993005/101‑1 (MQ=255)
                                          ccccGCGGGGCCGTGGATGGAAGCCTGCG.GGGGGAGGTGCCAGTGATAGA.CGCCGTGGA.CGCCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCgccgcc                                                            <  2:812530/101‑1 (MQ=37)
                                              gcgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCTGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATc                                                        <  2:2593347/101‑1 (MQ=255)
                                              gcgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATc                                                        <  1:4155781/101‑1 (MQ=255)
                                                gggggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGc                                                      <  1:2298608/101‑1 (MQ=255)
                                                 ggggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCa                                                     >  1:2016962/1‑101 (MQ=255)
                                                  gggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCAt                                                    >  2:4552547/1‑101 (MQ=255)
                                                   ggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATc                                                   >  1:1640122/1‑101 (MQ=255)
                                                    gCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATcc                                                  <  2:2347770/101‑1 (MQ=255)
                                                     cGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCg                                                 <  1:2495295/101‑1 (MQ=255)
                                                       gggATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGt                                               <  2:2485718/100‑1 (MQ=255)
                                                          aTCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGCA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGa                                            >  2:1595270/1‑101 (MQ=255)
                                                           tCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAg                                           <  1:2772901/101‑1 (MQ=255)
                                                            cGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAgc                                          >  1:3131049/1‑101 (MQ=255)
                                                                   tGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCgcgc                                   <  2:1773835/101‑1 (MQ=255)
                                                                    gcg.gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCa                                  <  2:2891851/101‑1 (MQ=255)
                                                                     cg.gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCgcgccg                                 >  2:4265924/1‑99 (MQ=255)
                                                                        gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTc                               >  1:2107206/1‑101 (MQ=255)
                                                                        gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTc                               >  1:3620172/1‑101 (MQ=255)
                                                                          gTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGc                             >  2:1119825/1‑101 (MQ=255)
                                                                               ttGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTcccc                        >  1:398629/1‑101 (MQ=255)
                                                                               tggCCAGTGATGGG.CGCGGTGGC.CGCCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTcccc                        <  2:4542717/99‑1 (MQ=37)
                                                                                tGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTccccc                       <  2:3087176/101‑1 (MQ=255)
                                                                                tGCCAGTGATAGA.CGCCGGGGA.CGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTccccc                       <  2:2319511/101‑1 (MQ=255)
                                                                                 gCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTcccccc                      <  2:3266435/101‑1 (MQ=255)
                                                                                  gcAGTGATAGG.CGCCGTGGA.CGTCGCTCCCCACGCGGGGGCGTGAATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCcg                     <  2:4600809/100‑1 (MQ=255)
                                                                                   cAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCcgc                    <  2:2048089/101‑1 (MQ=255)
                                                                                   cAGTGATAGA.CGCCGTGGA.CGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCcgc                    <  2:1677516/101‑1 (MQ=255)
                                                                                       gATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCgggg                >  1:4147111/1‑101 (MQ=255)
                                                                                        aTAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCggggg               >  2:2009993/1‑101 (MQ=255)
                                                                                        aTAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCggggg               >  2:96073/1‑101 (MQ=255)
                                                                                         tAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGc              >  2:1838052/1‑101 (MQ=255)
                                                                                          aGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCg             >  2:3161475/1‑101 (MQ=255)
                                                                                          aGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCg             >  1:1829408/1‑101 (MQ=255)
                                                                                          aGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCg             >  1:4235641/1‑101 (MQ=255)
                                                                                           gA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCgg            >  2:1402514/1‑100 (MQ=255)
                                                                                           gA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGt            >  1:2720519/1‑101 (MQ=255)
                                                                                            a.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTg           <  1:3345298/101‑1 (MQ=255)
                                                                                               cccGTGGA.CGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGa         <  1:72297/100‑1 (MQ=255)
                                                                                                 cGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc       <  1:2355471/101‑1 (MQ=255)
                                                                                                  gTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg      >  2:2979389/1‑101 (MQ=255)
                                                                                                  gTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg      >  2:2675619/1‑101 (MQ=255)
                                                                                                    ggA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa    >  1:3080000/1‑101 (MQ=255)
                                                                                                    ggA.CGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa    <  1:2595230/101‑1 (MQ=255)
                                                                                                     gA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa   >  2:3973886/1‑101 (MQ=255)
                                                                                                     gA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa   <  2:1230839/101‑1 (MQ=255)
                                                                                                     gA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa   <  1:2121695/101‑1 (MQ=255)
                                                                                                     gA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGAGCGaaa   >  2:4381606/1‑101 (MQ=255)
                                                                                                     gA.CGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa   <  1:4673698/101‑1 (MQ=255)
                                                                                                     gA.CGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa   <  2:1563313/101‑1 (MQ=255)
                                                                                                       |                                                                                                    
CTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC  >  NC_014034/1395730‑1395930

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