Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,839 C→A 13.2% intergenic (+651/‑280) 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,8390C→A13.2% 70.2 ‑1.3 ‑1.6249intergenic (+651/‑280)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (21/12):  ref base (112/104):  total (133/116)
Fisher's exact test strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.69e-02

Only 100 of 287 total aligned reads displayed.

ATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTC  >  NC_014034/1395739‑1395939
                                                                                                       |                                                                                                            
aTCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcc                                                                                                              >  2:1757837/1‑101 (MQ=255)
  cTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccc                                                                                                            <  1:1390594/101‑1 (MQ=255)
    ggCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                                                          <  1:3895310/101‑1 (MQ=255)
    ggCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                                                          <  1:2941816/101‑1 (MQ=255)
       gcAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCcgcg                                                                                                       >  1:3584298/1‑101 (MQ=255)
       gcAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCcgcg                                                                                                       <  1:3274760/101‑1 (MQ=255)
           aCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCggggg                                                                                                   >  1:522976/1‑101 (MQ=255)
            cGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGc                                                                                                  >  2:2550219/1‑101 (MQ=255)
            cGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGc                                                                                                  >  2:4438276/1‑101 (MQ=255)
            cGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGc                                                                                                  >  2:4496992/1‑101 (MQ=255)
            cGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGc                                                                                                  <  1:1017993/101‑1 (MQ=255)
              cggacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCg                                                                                                 >  2:1868981/11‑101 (MQ=0)
               aaaGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTg                                                                                               <  2:3420245/101‑1 (MQ=255)
                gacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCGTg                                                                                               >  2:1254106/9‑101 (MQ=0)
                aaGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTgg                                                                                              <  1:3718644/101‑1 (MQ=255)
                 aGTCTTTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGa                                                                                             <  1:1613199/101‑1 (MQ=255)
                  gTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                                            <  1:2034546/101‑1 (MQ=255)
                  gTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                                            <  1:2500458/101‑1 (MQ=255)
                   tCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                                           <  1:3619447/101‑1 (MQ=255)
                   tCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                                           <  1:2725235/101‑1 (MQ=255)
                    cTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                                          >  2:1169060/1‑101 (MQ=255)
                    cTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                                          <  2:2271939/101‑1 (MQ=255)
                    cTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                                          >  1:2873492/1‑101 (MQ=255)
                     tGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                                         <  2:190521/101‑1 (MQ=255)
                        cgtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                                      >  1:1323786/1‑101 (MQ=255)
                         gtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGCCGTGGATCGAAACt                                                                                     >  2:3108146/1‑101 (MQ=255)
                             cTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgcg.c                                                                                >  2:133945/1‑101 (MQ=255)
                              tCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGcg.cc                                                                               >  2:2086861/1‑101 (MQ=255)
                               ccccccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGcg.ccg                                                                              >  2:4115737/1‑101 (MQ=255)
                                 ccccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCg.ccgc.c                                                                           <  1:1117945/101‑1 (MQ=255)
                                 ccccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCg.ccgc.c                                                                           >  2:4345463/1‑101 (MQ=255)
                                  cccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.Ct                                                                          <  2:2668528/101‑1 (MQ=255)
                                     gcgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATc                                                                       <  1:2667739/101‑1 (MQ=255)
                                      cgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCg                                                                      <  2:3339383/101‑1 (MQ=255)
                                       gggggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGGGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CGGC.CCAtggt                                                                     >  1:3359348/1‑98 (MQ=37)
                                         gggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCAt                                                                   <  1:973256/101‑1 (MQ=255)
                                          ggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATc                                                                  <  1:1104417/101‑1 (MQ=255)
                                             gTGGATCGAAACCTGCG.GCGGGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCgg                                                               >  2:153187/1‑101 (MQ=255)
                                              tGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...t                                                           <  1:4207366/101‑1 (MQ=255)
                                              tGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...t                                                           >  2:3975464/1‑101 (MQ=255)
                                              tGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGAGAGA.CGCCGTGGA.CGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...t                                                           <  1:3978236/101‑1 (MQ=255)
                                               ggATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...Tc                                                          >  1:3088020/1‑101 (MQ=255)
                                                gATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCg                                                         >  1:1856567/1‑101 (MQ=255)
                                                gATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCg                                                         >  2:123555/1‑101 (MQ=255)
                                                     aaaCCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAgggc                                                    >  2:2440965/1‑101 (MQ=255)
                                                        ccTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGcgcc.gc                                                <  1:192333/101‑1 (MQ=255)
                                                         cTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGCCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCc.gcg                                               <  1:820591/101‑1 (MQ=255)
                                                                 gTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGc                                      >  1:2041478/1‑101 (MQ=255)
                                                                  tGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCt                                     <  2:2040900/101‑1 (MQ=255)
                                                                  tGAGTTGCCAGTGATAGA.CCCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCt                                     <  2:2458898/101‑1 (MQ=255)
                                                                     gTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTccc                                  <  1:3638513/101‑1 (MQ=255)
                                                                       tGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTccccc                                <  1:3471444/101‑1 (MQ=255)
                                                                        gCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTcccccc                               >  2:3205004/1‑101 (MQ=255)
                                                                             tGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCCCTCCCCCCGCggg                          >  2:3837015/1‑101 (MQ=255)
                                                                                 aGA.CGCCGGGGA.CGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCg                      <  2:4364221/101‑1 (MQ=255)
                                                                                  gA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCgg                     >  1:409313/1‑100 (MQ=255)
                                                                                     cGCCGTGGA.CGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTgg                   <  2:2218044/101‑1 (MQ=255)
                                                                                      gccaacctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCCCCCGCGGGGGCGTgg                   >  2:2750397/11‑101 (MQ=9)
                                                                                      gCCGTGGA.CGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGa                  <  1:153201/101‑1 (MQ=255)
                                                                                      ccggcccggcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCG.CGGC.GTCCCGAAGACGGTCTTCTGGTGCC.GCT.GC.GTCGCTCCCCCCGCGGGGGCg                      >  1:3759662/10‑101 (MQ=2)
                                                                                        caacctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCCCCCGCGGGGGCGTGGAt                 <  1:3345442/93‑1 (MQ=9)
                                                                                        caacctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.GTCCCGGAGACGC...T.GGGGGCC.TCG.CATGTCGCTCCCCCCGCGGGGGCGTGGAt                 >  1:1431952/9‑101 (MQ=9)
                                                                                        cGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATc                >  1:2295426/1‑101 (MQ=255)
                                                                                         gcccggcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCG.CGGC.GTCCCGAAGACGGTCTTCTGGTGCC.GCT.GC.GTCGCTCCCCCCGCGGGGGCGTgg                   <  2:3650078/95‑1 (MQ=2)
                                                                                         gTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCg               >  2:2897279/1‑101 (MQ=255)
                                                                                         aacctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCCCCCGCGGGGGCGTGGATc                <  2:4074281/94‑1 (MQ=9)
                                                                                          tGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGa              <  1:4438571/101‑1 (MQ=255)
                                                                                          acctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCCCCCGCGGGGGCGTGGATCg               <  1:4672341/95‑1 (MQ=9)
                                                                                           cctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCCCCCGCGGGGGCGTGGATCGa              <  1:3535242/96‑1 (MQ=9)
                                                                                            gA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa            <  1:217556/101‑1 (MQ=255)
                                                                                            ctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCCCCCGCGGGGGCGTGGATCGaa             >  2:3705106/5‑101 (MQ=9)
                                                                                            ctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCCCCCGCGGGGGCGTGGATCGaa             >  1:1397848/5‑101 (MQ=9)
                                                                                             tatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCTCCCGCGGGGGCGTGGGTCGgaa            >  2:4264835/4‑101 (MQ=9)
                                                                                             tatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCCCCCCCGGGGGGGTGGATCGaaa            >  1:2770144/4‑101 (MQ=9)
                                                                                             ggcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCG.CGGC.GTCCCGAAGACGGTCTTCTGGTGCC.GCT.GC.GTCGCTCCCCCCGCGGGGGCGTGGATCg               >  1:4523734/3‑101 (MQ=2)
                                                                                             ggcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCG.CGGC.GTCCCGAAGACGGTCTTCTGGTGCC.GCT.GC.GTCGCTCCCCCCGCGGGGGCGTGGATCg               <  1:1772665/99‑1 (MQ=2)
                                                                                             cccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGTGCGCGCCAT.GAATTC.A...TCGAGCGCCATCGAGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa              <  1:393502/99‑1 (MQ=2)
                                                                                             a.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc           <  2:1378148/101‑1 (MQ=255)
                                                                                              gcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCG.CGGC.GTCCCGAAGACGGTCTTCTGGTGCC.GCT.GC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa              <  2:2152904/100‑1 (MQ=2)
                                                                                              atgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc           >  2:2852548/3‑101 (MQ=9)
                                                                                               tgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAca          >  2:3056295/2‑100 (MQ=9)
                                                                                               tgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAca          <  2:1484605/100‑2 (MQ=9)
                                                                                               cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACg          <  2:2703636/101‑1 (MQ=255)
                                                                                               cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACg          >  1:4627881/1‑101 (MQ=255)
                                                                                                gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GACGCTCCCCCCACGGGGGCGTGGAAGGAAACgc         >  1:1223822/1‑100 (MQ=38)
                                                                                                gTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...AGGGGC.CC.TCG.CACGTCGCTCCCCCCGCCGGGGCGTGGATCGGAAcaa         >  2:742595/1‑99 (MQ=9)
                                                                                                gTCCCTCCCCCCGCGGGGGCGTGGATCGAAACCCCG.CGGC.GTCCCGAAGACGGTCTTCTGGTGCC.GCT.GC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa            <  1:632022/101‑1 (MQ=2)
                                                                                                 ccGCTCCCCCCGCGGGGGCGTTGATCCAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAt        <  2:4637791/100‑1 (MQ=255)
                                                                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATg       >  2:3158240/1‑101 (MQ=255)
                                                                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGc      <  2:1304353/101‑1 (MQ=255)
                                                                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGc      >  1:591223/1‑101 (MQ=255)
                                                                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGc      >  2:3587832/1‑101 (MQ=255)
                                                                                                   gCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGc      >  2:1945729/1‑101 (MQ=9)
                                                                                                   gCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGc      <  1:2608907/101‑1 (MQ=9)
                                                                                                   gCTCACCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGc      >  1:2337891/1‑101 (MQ=9)
                                                                                                    ccccccccGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCg     <  1:3836971/99‑1 (MQ=255)
                                                                                                    cTCGCCCCGCGGGGGCGTGGATCGAAACTGTG.CCGG.TTCCCGGAGACGC...TGGGGC.CC.TCG.CACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCg     <  2:695300/101‑1 (MQ=9)
                                                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGt    >  2:3742678/1‑101 (MQ=255)
                                                                                                      ccccccGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGtt   >  1:1699364/1‑101 (MQ=255)
                                                                                                       cccccGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTc  <  2:475371/101‑1 (MQ=255)
                                                                                                       |                                                                                                            
ATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCG.CCGC.CTATCGCATCCGG...TCGAGCGCC.GCG.CA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTC  >  NC_014034/1395739‑1395939

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG