Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,737 Δ1 bp 9.9% intergenic (+549/‑382) 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,7370G→. 9.9% 85.5 ‑1.1 ‑0.0292intergenic (+549/‑382)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (11/18):  ref base (128/135):  total (139/153)
Fisher's exact test strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01

Only 100 of 312 total aligned reads displayed.

CTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCT  >  NC_014034/1395637‑1395837
                                                                                                     |                                                                                                           
 cccccccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGGGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                                                            <  1:1760071/100‑1 (MQ=255)
 cccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                                                            <  2:216899/100‑1 (MQ=255)
 cccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                                                            <  2:1599850/100‑1 (MQ=255)
        gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..gc                                                                                                   >  1:1235156/1‑101 (MQ=255)
        gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..gc                                                                                                   <  1:3567444/101‑1 (MQ=255)
             ggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCaaa                                                                                              >  2:521171/1‑101 (MQ=255)
                gTGGATCGAACCATCTGCATGTAGTGGGCGGGCGGTTGGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACgg                                                                                           <  2:4046148/101‑1 (MQ=37)
                gTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAGCGAAACT.GGTATGGATCTG..GCGCAAACgg                                                                                           >  2:479927/1‑101 (MQ=255)
                 tGGCTCGAAACATCTGCTTGTAGGCGGCGGGGAGTGGGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGAGCTG..GCGCAAACGG..a                                                                                        <  1:334782/101‑1 (MQ=25)
                  ggATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..aa                                                                                       >  2:2285565/1‑101 (MQ=255)
                  ggATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..aa                                                                                       <  2:3976818/101‑1 (MQ=255)
                   ggtCGAAGCATCTGGTGGTGGGCGGCGGGCGGTTGGAGCCGGAGTGGCGCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..aaa                                                                                      <  1:3611064/99‑1 (MQ=21)
                   gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..aaa                                                                                      <  1:2981204/101‑1 (MQ=255)
                   gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..aaa                                                                                      >  1:138248/1‑101 (MQ=255)
                     tCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGt                                                                                    >  1:2593528/1‑101 (MQ=255)
                      cGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTc                                                                                   >  2:2767760/1‑101 (MQ=255)
                      cGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTc                                                                                   >  1:516782/1‑101 (MQ=255)
                       gAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCt                                                                                  >  2:1896613/1‑101 (MQ=255)
                        aaaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g                                                                                <  1:89249/101‑1 (MQ=255)
                           cATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.tcg                                                                            >  1:2817565/1‑101 (MQ=255)
                             tCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGCCT.g.ccgtc                                                                          >  2:3679848/1‑101 (MQ=255)
                              cTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.tcgtcg                                                                         >  1:2398007/1‑101 (MQ=255)
                                  aTGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.cc                                                                    >  1:3816673/1‑101 (MQ=255)
                                   ggTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AACGTCT.G.TCGTCGCT.ccc                                                                   <  1:4413511/100‑1 (MQ=255)
                                    gTAGTCGGGGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.cccc                                                                  <  1:1161287/101‑1 (MQ=255)
                                      aGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.cccccc                                                                >  1:73292/1‑101 (MQ=255)
                                       gTCGGGGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCcg                                                               <  1:3210529/101‑1 (MQ=255)
                                         cggcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCcgcg                                                             >  1:2091024/1‑101 (MQ=255)
                                             gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCggggg                                                         >  1:3318497/1‑101 (MQ=255)
                                             gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCggggg                                                         >  1:4037085/1‑101 (MQ=255)
                                               gCAGTGCGAGCCGCAGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCCGCGGGGGc                                                        <  2:2617497/101‑1 (MQ=38)
                                                 aGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTg                                                     <  1:323472/101‑1 (MQ=255)
                                                 aGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTg                                                     >  2:3345035/1‑101 (MQ=255)
                                                  gtgcGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTgg                                                    <  1:2186852/98‑1 (MQ=255)
                                                  gTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTgg                                                    >  2:3049997/1‑101 (MQ=255)
                                                  gTTCGAGCCGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTgg                                                    <  2:4180378/101‑1 (MQ=255)
                                                   ttCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGa                                                   <  1:4930463/101‑1 (MQ=255)
                                                    tcaagattcAGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTg                                                     <  2:3251915/93‑1 (MQ=2)
                                                     caagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTgg                                                    <  1:4882990/94‑1 (MQ=2)
                                                     cGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATc                                                 <  2:721739/101‑1 (MQ=255)
                                                      gcgcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGAt                                                  <  1:98905/93‑1 (MQ=0)
                                                      gAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCg                                                <  1:2328698/101‑1 (MQ=255)
                                                        gcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCg                                                <  1:3977772/95‑1 (MQ=0)
                                                        gattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.GGGTCGCT.CCCCCCGCGGGGGCGTGGATc                                                 >  2:506593/5‑101 (MQ=2)
                                                        gCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGaa                                              <  2:4364384/101‑1 (MQ=255)
                                                         attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCg                                                >  1:212518/4‑101 (MQ=2)
                                                         attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCg                                                <  1:4809219/98‑1 (MQ=2)
                                                          ttcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTACGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGCGGATCGa                                               >  1:4263048/3‑101 (MQ=0)
                                                          ttcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGa                                               <  2:1436274/99‑1 (MQ=2)
                                                          cGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAc                                            <  1:69988/101‑1 (MQ=255)
                                                           tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGaa                                              >  2:3575117/2‑101 (MQ=2)
                                                           ctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGaaa                                             <  2:4923984/98‑1 (MQ=2)
                                                             aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAc                                            <  1:2866293/101‑1 (MQ=2)
                                                              gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAca                                           >  1:1683291/1‑100 (MQ=2)
                                                              gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTgc                                        >  1:1314846/1‑101 (MQ=255)
                                                              gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTgc                                        <  2:3744672/101‑1 (MQ=255)
                                                              gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAGAca                                           <  1:1896147/101‑2 (MQ=2)
                                                               tCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTgcg                                       >  2:123139/1‑101 (MQ=255)
                                                               tCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTgcg                                       <  1:4286753/101‑1 (MQ=255)
                                                                cGCTCCCCCCGCGGGGGTGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCtat                                        <  2:4539446/101‑3 (MQ=2)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCtat                                        <  2:4289950/101‑3 (MQ=2)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTgcgg                                      >  1:452418/1‑101 (MQ=255)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTgcgg                                      <  1:4564855/101‑1 (MQ=255)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTgcgg                                      >  2:4903140/1‑101 (MQ=255)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCCgcgg                                      >  1:2222896/1‑101 (MQ=255)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcaccc                                        >  2:4984721/1‑97 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcaccc                                        <  1:3638760/101‑5 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCtatg                                       <  1:1057341/101‑4 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCtatg                                       <  1:155601/101‑4 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCtatg                                       >  1:3715007/1‑98 (MQ=2)
                                                                  ccccccccGGGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCC.CCCCCCGCGGGGGCGTGGATCGAAAcacccg                                       <  1:2651617/99‑6 (MQ=2)
                                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcacccg                                       >  2:2897307/1‑96 (MQ=2)
                                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCtatga                                      <  1:4987852/101‑5 (MQ=2)
                                                                   tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcacccgg                                      >  2:4618857/1‑95 (MQ=2)
                                                                   tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcgtatgag                                     <  1:2066794/101‑8 (MQ=2)
                                                                   tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCtatgag                                     <  1:3037247/101‑6 (MQ=2)
                                                                   cccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcacccgg                                      <  2:3894399/100‑7 (MQ=0)
                                                                    ccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcacccggc                                     >  2:1542296/1‑94 (MQ=0)
                                                                    ccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcacccggc                                     <  2:4338879/101‑8 (MQ=0)
                                                                       cccGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGGCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGt                               <  1:4623008/101‑1 (MQ=255)
                                                                        ccGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGtt                              <  2:3792730/101‑1 (MQ=255)
                                                                         cgcgGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAgatc                             >  2:5017295/1‑98 (MQ=255)
                                                                             ggggCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAg                         >  2:1964466/1‑101 (MQ=255)
                                                                              gggCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAgg                        >  2:2848348/1‑100 (MQ=255)
                                                                              gggCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGt                        <  2:2090876/101‑1 (MQ=255)
                                                                                 cGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGAt                     <  2:2696841/101‑1 (MQ=255)
                                                                                   tGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAg                   <  1:274270/101‑1 (MQ=255)
                                                                                    ggATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGa                  <  1:2685758/101‑1 (MQ=255)
                                                                                     gATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGACAGAc                 <  2:4273559/101‑1 (MQ=255)
                                                                                         gAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCg             >  2:2298444/1‑101 (MQ=255)
                                                                                         gAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGCGATAGACGACg             >  2:4598871/1‑101 (MQ=255)
                                                                                          aaaCT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGt            >  2:1840929/1‑101 (MQ=255)
                                                                                             cT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGa         >  2:2219080/1‑101 (MQ=255)
                                                                                                 gTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGt      <  2:3621307/101‑1 (MQ=255)
                                                                                                 gTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGt      >  2:1107223/1‑101 (MQ=255)
                                                                                                  tATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTc     >  1:3443522/1‑101 (MQ=255)
                                                                                                   aTGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACgctg    >  1:2485892/1‑98 (MQ=255)
                                                                                                   aTGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCg    <  2:2016822/101‑1 (MQ=255)
                                                                                                     ggATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGTCGTCGCt  >  2:4501529/1‑101 (MQ=255)
                                                                                                     ggATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTCATAGACGCCGTGGACGTCGCt  >  2:3399631/1‑101 (MQ=255)
                                                                                                     |                                                                                                           
CTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCT  >  NC_014034/1395637‑1395837

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