Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,533 T→A 20.9% intergenic (+345/‑586) 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,5330T→A20.9% 230.7 ‑3.3 ‑1.6410intergenic (+345/‑586)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (46/29):  ref base (136/149):  total (207/203)
Fisher's exact test strand distribution p-value = 3.85e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.51e-02

Only 100 of 474 total aligned reads displayed.

GGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCG  >  NC_014034/1395433‑1395633
                                                                                                    |                                                                                                        
  tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctt                                                                                                       >  2:4804208/1‑100 (MQ=2)
  tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtct                                                                                                        >  2:1560552/1‑101 (MQ=255)
  tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtct                                                                                                        <  2:5773985/101‑1 (MQ=255)
   cGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgcttt                                                                                                      >  1:3670604/1‑99 (MQ=2)
   cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtctc                                                                                                       >  1:5623428/1‑101 (MQ=255)
    gCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttc                                                                                                     >  2:4703055/1‑98 (MQ=2)
    gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg                                                                                                      <  2:5632324/101‑1 (MQ=255)
     cTCCCCCCGGGGGGGCGTGGATCGAAACATCGCCGCTTTGGGGGGCGACAGGCTGGTAACGTGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACAtctccg                                                                                                     <  2:4125728/101‑3 (MQ=37)
     cTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttca                                                                                                    >  1:985534/1‑97 (MQ=2)
      tCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg.gg                                                                                                   >  2:3501322/1‑100 (MQ=255)
       ccccccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTg                                                                                                  >  2:2870188/1‑101 (MQ=255)
        cccccGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaatg                                                                                                 <  2:2081577/101‑8 (MQ=2)
         ccccGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaatgg                                                                                                <  2:5556127/101‑9 (MQ=2)
         ccccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGcc                                                                                                >  1:3170801/1‑101 (MQ=255)
          cccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCg                                                                                               <  2:4739575/101‑1 (MQ=255)
          cccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCg                                                                                               >  2:1363716/1‑101 (MQ=255)
            cgcgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTc                                                                                             <  2:3397558/101‑1 (MQ=255)
             gcgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCACCGCGGGGGCGTGGATCGAAACATCTCG.GTGCGGTCg                                                                                            >  1:5560049/1‑101 (MQ=255)
              cgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGc                                                                                           >  1:2720527/1‑101 (MQ=255)
               gggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGcc                                                                                          >  1:4016905/1‑100 (MQ=255)
               gggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCt                                                                                          >  1:1151582/1‑101 (MQ=255)
                ggggCGTTGATCGAAACATCGCCGCTTTGGGGGTCGACAGGCTGGTAACGTGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTc                                                                                         <  1:5436768/101‑1 (MQ=255)
                gcggCGTGGATCGAACCATCGCCACTATGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTc                                                                                         <  2:2649431/99‑1 (MQ=255)
                 gggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCg                                                                                        <  1:1300119/101‑1 (MQ=255)
                  ggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.t                                                                                      <  2:4504263/101‑1 (MQ=255)
                  ggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.t                                                                                      >  1:3377459/1‑101 (MQ=255)
                  ggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.t                                                                                      >  2:2451526/1‑101 (MQ=255)
                    cGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGTATGCCGTcttctgct                                                                                      >  1:5293820/1‑94 (MQ=37)
                       ggATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGc                                                                                 <  2:4554082/101‑1 (MQ=255)
                            gAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.cc                                                                           >  1:1483847/1‑101 (MQ=255)
                             aaaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCt                                                                          >  1:2834114/1‑101 (MQ=255)
                             aaaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCt                                                                          >  2:2234151/1‑101 (MQ=255)
                              aactgcgccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCGCCCCCCGCGGGGGCGTGGATCGCAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.g                                                                        <  2:1277447/96‑1 (MQ=255)
                              aaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.g                                                                        <  1:1326460/101‑1 (MQ=255)
                                   cgccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCtcgt                                                                   <  1:3501539/101‑1 (MQ=255)
                                    gccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCtcgtc                                                                  >  1:4592800/1‑101 (MQ=255)
                                     ccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCtcgtcg                                                                 >  1:3349507/1‑101 (MQ=255)
                                        cTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTc                                                              <  2:2944359/101‑1 (MQ=255)
                                         tttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTcc                                                             <  1:2138662/101‑1 (MQ=255)
                                            ggCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTccccc                                                          <  1:3578388/101‑1 (MQ=255)
                                             gCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTcccccc                                                         <  1:5053436/101‑1 (MQ=255)
                                                  cGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCggg                                                    <  1:2234273/101‑1 (MQ=255)
                                                   gACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCgggg                                                   >  1:1699998/1‑101 (MQ=255)
                                                   gACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCgggg                                                   <  1:5153307/101‑1 (MQ=255)
                                                    aCAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCggggg                                                  <  2:3641430/101‑1 (MQ=255)
                                                     cAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGc                                                 <  1:2870349/101‑1 (MQ=255)
                                                         ccggTACCGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCGAT.GATCAAATCCTCA.ATCGGCCGG.GCTCGATCGTCGCTCCCCCCGCGGGGGCGTg                                              <  1:4350015/99‑1 (MQ=2)
                                                         cTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTAGCTCG.TCCAGCCGG.CCT.GATCGTCGCTCCCCCCCCGGGGGCGTgg                                             >  1:4866181/1‑101 (MQ=255)
                                                         agcggctgcggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAG.CCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTgg                                             >  2:3034590/11‑101 (MQ=0)
                                                          aggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCG.GCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGa                                            <  1:2499469/92‑1 (MQ=0)
                                                           ggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCG.GCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGAt                                           <  1:2458198/93‑1 (MQ=2)
                                                           ggTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGAt                                           >  2:1060318/1‑101 (MQ=255)
                                                            gTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATc                                          <  2:464786/101‑1 (MQ=255)
                                                             tgccgcgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGAT.GATCAAATCCGCA.ATCGGCCGG.GCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATc                                          <  2:827328/95‑1 (MQ=2)
                                                             tAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                         >  1:611848/1‑101 (MQ=255)
                                                               ggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCG.GCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                       >  1:3473522/5‑101 (MQ=2)
                                                                cacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                      >  2:5455354/4‑101 (MQ=2)
                                                                cacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                      <  2:654011/98‑1 (MQ=2)
                                                                cGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                      <  2:5019911/101‑1 (MQ=255)
                                                                cGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                      <  1:1877176/101‑1 (MQ=255)
                                                                  tgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacc                                    >  2:4372877/1‑101 (MQ=255)
                                                                  ggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAG.CCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacc                                    <  2:149606/100‑1 (MQ=2)
                                                                  cgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAact                                    <  2:1758998/100‑2 (MQ=2)
                                                                   gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacca                                   >  1:847506/1‑101 (MQ=255)
                                                                   gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAG.CCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccc                                   >  1:2968089/1‑100 (MQ=2)
                                                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAG.CCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGTGCGTGGATCGAAacccc                                  >  1:682591/1‑99 (MQ=2)
                                                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAG.CCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacccc                                  <  2:3970681/101‑3 (MQ=2)
                                                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAG.CCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacccc                                  <  1:767547/101‑3 (MQ=2)
                                                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAG.CCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCACGCGGGGGCGTGGCTCGCAacacc                                  >  2:2036539/1‑98 (MQ=0)
                                                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCCTGGATCGAAACCACCt                                >  2:5231980/1‑101 (MQ=255)
                                                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgt                                <  1:471534/101‑6 (MQ=2)
                                                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGGCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgt                                >  2:2333436/1‑96 (MQ=2)
                                                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAG.CCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgc                                >  1:4472075/1‑97 (MQ=2)
                                                                        ccccccgGCGGGGGCGGGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgc                              <  1:2411658/100‑8 (MQ=0)
                                                                         ccccccGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCttcc                             >  2:4918852/1‑99 (MQ=255)
                                                                         ccccccGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgcg                             <  1:262140/101‑9 (MQ=2)
                                                                         ccccccGCGGGGGCGTGGATCGAAACA.CCCG.GCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctgg                             >  2:4085895/1‑93 (MQ=2)
                                                                          cccccGCGGGGGCGTGGATCGAAACAAAAAG.CCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcgccg                            >  2:766770/1‑93 (MQ=0)
                                                                          cccccGCGCGGGCTTGGGTCGAAGCGTCTCG.ATGCCGTCGCCCT.TCCAGCCGG.CCT.GCTGGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCg                            <  2:5293530/101‑1 (MQ=25)
                                                                           cccgGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgcgcg                           <  2:1691699/101‑11 (MQ=0)
                                                                           ccccGCGGGGGCGTGGATCGAAACAAAAAG.CCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcggcgt                           <  1:3755437/101‑10 (MQ=2)
                                                                           ccccGCGGGGGCGTGGATCGAAACA.CCCG.GCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggct                           >  2:1168536/1‑91 (MQ=2)
                                                                            cccGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGa                          <  2:1983109/101‑1 (MQ=255)
                                                                                 gggggCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                     >  2:2557926/1‑101 (MQ=255)
                                                                                 gggggCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                     >  2:1854467/1‑101 (MQ=255)
                                                                                    ggCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACt                  >  2:3776007/1‑101 (MQ=255)
                                                                                    ggCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACt                  >  1:2204582/1‑101 (MQ=255)
                                                                                    ggCGTGGATCGAAACATCTCG.GTGCCGTCGCTCC.TCCCGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACt                  <  1:4954920/101‑1 (MQ=255)
                                                                                     gCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACtt                 <  1:4165263/101‑1 (MQ=255)
                                                                                     gCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACtt                 <  1:2300608/101‑1 (MQ=255)
                                                                                     gCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACtt                 <  1:1006253/101‑1 (MQ=255)
                                                                                         ggATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcgc             >  2:5073312/1‑101 (MQ=255)
                                                                                         ggATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCAATTTcgc             >  1:1564393/1‑101 (MQ=255)
                                                                                             cGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCgggg         >  1:4664472/1‑101 (MQ=255)
                                                                                             cGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCgggg         <  1:391367/101‑1 (MQ=255)
                                                                                               aaaCATTTCG.CTGCCCTGGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCCTTGCGGATTGTCACTTTCGCGGGGCg       <  2:2792827/101‑1 (MQ=37)
                                                                                               aaaCATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCg       <  2:431107/101‑1 (MQ=255)
                                                                                                 aCATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGaa     <  1:2055132/101‑1 (MQ=255)
                                                                                                 aCATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGaa     <  2:92620/101‑1 (MQ=255)
                                                                                                    tctcG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCg  <  1:1087997/101‑1 (MQ=255)
                                                                                                    |                                                                                                        
GGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCG  >  NC_014034/1395433‑1395633

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