Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,763 C→G 17.7% intergenic (+575/‑356) 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,7630C→G17.7% 134.6 ‑2.2 0.0326intergenic (+575/‑356)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (20/38):  ref base (124/144):  total (144/182)
Fisher's exact test strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 372 total aligned reads displayed.

ACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC  >  NC_014034/1395663‑1395863
                                                                                                            |                                                                                                       
aCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.g.tc                                                                                                         <  1:1854577/101‑1 (MQ=255)
  aTCTGCAGGTAGTCGGCGGGCAGTCCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.g.tcgt                                                                                                       <  2:4097254/101‑1 (MQ=255)
   tCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.g.tcgtc                                                                                                      <  2:3016019/101‑1 (MQ=255)
   tCTGCATGTAGGCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.g.tcgtc                                                                                                      <  1:3792757/101‑1 (MQ=255)
    cTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.gcgtcg                                                                                                     >  2:2000575/1‑101 (MQ=255)
      gctcgtggTCGGCGGGGGGCTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCt                                                                                                   <  2:5325610/94‑1 (MQ=37)
      gCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGcc                                                                                                   >  1:3102790/1‑100 (MQ=255)
      gCATGGAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCt                                                                                                   <  1:4846620/101‑1 (MQ=255)
       cATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTc                                                                                                  >  1:536504/1‑101 (MQ=255)
         tGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTccc                                                                                                <  2:4073200/101‑1 (MQ=255)
          gTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGGCT.G.TCGTCGCTcccc                                                                                               >  2:2154292/1‑101 (MQ=255)
           tAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTccccc                                                                                              >  1:3696244/1‑101 (MQ=255)
            aGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccccc                                                                                             <  1:928122/101‑1 (MQ=255)
               cggcggGCAGTGCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCcgcg                                                                                          <  1:5530866/101‑1 (MQ=255)
                   gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCggggg                                                                                      >  2:560231/1‑101 (MQ=255)
                       aGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTg                                                                                  >  2:3793058/1‑101 (MQ=255)
                       aGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTg                                                                                  >  1:4009694/1‑101 (MQ=255)
                        tgtgcaGCCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTgg                                                                                 <  2:3618546/96‑1 (MQ=255)
                        gTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTgg                                                                                 >  1:196773/1‑101 (MQ=255)
                         ttCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGa                                                                                <  2:5761060/101‑1 (MQ=255)
                           caagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTgg                                                                                 >  2:1821220/8‑101 (MQ=2)
                           cGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                              <  2:5226279/101‑1 (MQ=255)
                           cGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                              >  1:472028/1‑101 (MQ=255)
                            gcgcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.G.ACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                               <  1:4890815/93‑1 (MQ=0)
                            aagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGa                                                                                <  1:3569634/95‑1 (MQ=2)
                              gcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.G.ACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                             >  1:2555406/7‑101 (MQ=0)
                               ccGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                          <  2:1663106/101‑1 (MQ=255)
                               attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                             <  1:1194636/98‑1 (MQ=2)
                               attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                             <  2:3171379/98‑1 (MQ=2)
                               attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                             <  1:2679353/98‑1 (MQ=2)
                                 tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           <  2:3757600/100‑1 (MQ=2)
                                 gCAGTCGCTCCCACCGCGGGGGCGTGGATCGAAACT.GGTATGG.CTCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                                                        >  1:3133615/1‑101 (MQ=255)
                                  tggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.G.ACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                         <  2:4029717/99‑1 (MQ=2)
                                   ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.G.ACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                                                        >  2:5300578/2‑101 (MQ=2)
                                   aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                         >  1:2374453/1‑101 (MQ=2)
                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAca                                                                        <  1:4525489/101‑2 (MQ=2)
                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAca                                                                        <  2:1047468/101‑2 (MQ=2)
                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.ATGGCCAGAT..CCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCt                                                                       <  2:696738/101‑1 (MQ=2)
                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.G.ACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCt                                                                       <  2:5368573/101‑1 (MQ=2)
                                    ctCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAca                                                                        <  2:499125/100‑2 (MQ=2)
                                     tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcac                                                                       <  2:2724002/101‑3 (MQ=2)
                                     tCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.G.ACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCta                                                                      >  2:789097/1‑100 (MQ=2)
                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.G.ACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtat                                                                     >  1:4591808/1‑99 (MQ=2)
                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.G.ACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtat                                                                     >  1:1687567/1‑99 (MQ=2)
                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcg.g                                                                  >  2:3314413/1‑101 (MQ=255)
                                      cGCTCCCCCCGCGGGGGCGGGGATCGAAACC.GGTACGG.ATCTG..GGGTAAACGG..AAGGGCT.C.GAGTGGCTCCCCACGCGGGGGGGGGGGGCGAAAcaagcag                                                                   >  2:1203382/1‑95 (MQ=12)
                                      cGCGCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.G.ACCAGACCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtat                                                                     <  2:4955516/101‑3 (MQ=2)
                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.G.ACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                                                    >  1:1576866/1‑98 (MQ=2)
                                         tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGGGGATCGAAAcacccgg                                                                   >  2:4887347/1‑95 (MQ=2)
                                         tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.G.ACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatgag                                                                  >  1:606988/1‑96 (MQ=2)
                                         tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.G.ACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatgag                                                                  <  2:5061920/101‑6 (MQ=2)
                                         tCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGt                                                               <  2:2913854/101‑1 (MQ=255)
                                          ccccccGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggc                                                                  <  1:4785185/101‑8 (MQ=0)
                                          ccccccGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggc                                                                  >  1:3869007/1‑94 (MQ=0)
                                           cccccGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGa                                                             <  1:5097431/101‑1 (MQ=255)
                                             cccGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggcgcc                                                               <  1:3957830/101‑11 (MQ=2)
                                             cccGCGGGGGCGTGGATCGAAACTCGGCAAGG.ACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggcgcc                                                               <  2:1499589/101‑11 (MQ=2)
                                              ccGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGtt                                                          >  2:882449/1‑101 (MQ=255)
                                                  gggggCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTTGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCa                                                      >  2:238112/1‑101 (MQ=255)
                                                    gggCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGt                                                    <  2:4623966/101‑1 (MQ=255)
                                                    gggCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGt                                                    <  1:1518457/101‑1 (MQ=255)
                                                     ggCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTg                                                   <  1:2922227/101‑1 (MQ=255)
                                                     ggCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTg                                                   >  2:1086876/1‑101 (MQ=255)
                                                      gCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGa                                                  >  2:3422218/1‑101 (MQ=255)
                                                      gCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGa                                                  <  1:1361793/101‑1 (MQ=255)
                                                        gTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATa                                                >  2:5705846/1‑101 (MQ=255)
                                                        gTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATa                                                <  1:1086951/101‑1 (MQ=255)
                                                         tGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAg                                               >  1:2220312/1‑101 (MQ=255)
                                                             tCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGc                                          <  1:3194997/101‑1 (MQ=255)
                                                               gAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCg                                        >  2:3267762/1‑101 (MQ=255)
                                                                 aaCT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTg                                      >  1:3271256/1‑101 (MQ=255)
                                                                  aCT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTgg                                     <  2:2460646/101‑1 (MQ=255)
                                                                   cT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                    >  1:2095772/1‑101 (MQ=255)
                                                                   cT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                    <  1:560281/101‑1 (MQ=255)
                                                                   cT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                    <  1:4146816/101‑1 (MQ=255)
                                                                    t.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.c                                  <  2:257769/101‑1 (MQ=255)
                                                                        tATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTc                               >  2:2291426/1‑101 (MQ=255)
                                                                         aTGG.ATCTG..GCTCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCg                              <  2:4717644/101‑1 (MQ=255)
                                                                           gg.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGaagatcgg                             >  2:450290/1‑94 (MQ=255)
                                                                                cTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGATGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccc                        >  1:4939063/1‑101 (MQ=255)
                                                                                 tG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTccccc                       >  2:3693871/1‑101 (MQ=255)
                                                                                  g..gCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                      <  1:798477/101‑1 (MQ=255)
                                                                                       gcAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCcccg                   >  1:3373228/1‑101 (MQ=255)
                                                                                        cAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                  >  2:421550/1‑101 (MQ=255)
                                                                                        cAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                  >  2:3450623/1‑101 (MQ=255)
                                                                                         aaaCGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCggg                 >  2:3063034/1‑101 (MQ=255)
                                                                                         aaaCGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCggg                 >  2:1386939/1‑101 (MQ=255)
                                                                                          aaCGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTccccccagatcg                >  2:657197/1‑95 (MQ=255)
                                                                                          aaCGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgggg                <  2:3139127/101‑1 (MQ=255)
                                                                                            cGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGc              <  1:665394/101‑1 (MQ=255)
                                                                                             gg..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCg             <  1:1716095/101‑1 (MQ=255)
                                                                                                 aaaGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTg           >  1:1023535/1‑101 (MQ=255)
                                                                                                   ggacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCGt            <  1:2484643/92‑1 (MQ=0)
                                                                                                   aGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGa         <  1:4392092/101‑1 (MQ=255)
                                                                                                    gacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCGTg           <  1:2996774/93‑1 (MQ=0)
                                                                                                    gacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCGTg           >  1:5039368/9‑101 (MQ=0)
                                                                                                    gTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGGGGAt        >  2:5530636/1‑101 (MQ=255)
                                                                                                      cT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCg      <  1:1243645/101‑1 (MQ=255)
                                                                                                         g.tcgtcGCTCCCCCCGCGGGGGCGTGGATCGAAAGCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGaa    >  1:5761005/1‑101 (MQ=255)
                                                                                                            cgtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc  >  1:5303001/1‑101 (MQ=255)
                                                                                                            |                                                                                                       
ACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGG.ATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC  >  NC_014034/1395663‑1395863

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