Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,839 C→A 12.0% 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→A12.0% 100.1 ‑0.4 ‑1.5411intergenic (+651/‑280)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (24/26):  ref base (186/175):  total (210/201)
Fisher's exact test strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.40e-02

Only 100 of 461 total aligned reads displayed.

ATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTC  >  NC_014034/1395739‑1395939
                                                                                                       |                                                                                                        
aTCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcc                                                                                                          <  1:278120/101‑1 (MQ=255)
 tCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGGGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTccc                                                                                                         <  1:1408833/101‑1 (MQ=255)
  cTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccc                                                                                                        >  2:2678826/1‑101 (MQ=255)
  cTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccc                                                                                                        >  1:5093472/1‑101 (MQ=255)
   tGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTccccc                                                                                                       <  2:5970586/101‑1 (MQ=255)
    ggCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                                                      <  1:470451/101‑1 (MQ=255)
    ggCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                                                      <  1:6838086/101‑1 (MQ=255)
     gcgcAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCcg                                                                                                     <  1:513858/101‑1 (MQ=255)
      cgcAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCcgc                                                                                                    <  1:1004213/101‑1 (MQ=255)
       gcAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCcgcg                                                                                                   >  1:534476/1‑101 (MQ=255)
        cAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                                                                  <  2:1983434/101‑1 (MQ=255)
         aaaCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCggg                                                                                                 >  1:6516110/1‑101 (MQ=255)
          aaCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgggg                                                                                                <  2:5317014/101‑1 (MQ=255)
           aCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCggggg                                                                                               >  1:3062324/1‑101 (MQ=255)
           aCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCggggg                                                                                               <  2:6542426/101‑1 (MQ=255)
            cGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGc                                                                                              <  2:2791213/101‑1 (MQ=255)
            cGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGc                                                                                              >  2:4658427/1‑101 (MQ=255)
             ggAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCg                                                                                             >  1:3676897/1‑101 (MQ=255)
             ggAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCg                                                                                             >  1:4584742/1‑101 (MQ=255)
              gAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCgg                                                                                            >  2:5520865/1‑100 (MQ=255)
               ggacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCGt                                                                                            >  1:3855375/10‑101 (MQ=0)
               aaaGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTg                                                                                           >  1:2616209/1‑101 (MQ=255)
                acgTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTgg                                                                                          <  1:7010043/99‑1 (MQ=255)
                  gTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                                        >  2:3885299/1‑101 (MQ=255)
                     tGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                                     >  1:4385331/1‑101 (MQ=255)
                     tGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                                     <  1:6261863/101‑1 (MQ=255)
                     tGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                                     <  2:5596033/101‑1 (MQ=255)
                      gtcgtcGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                                    >  2:1287157/1‑101 (MQ=255)
                      gtcgtcGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                                    <  2:3595203/101‑1 (MQ=255)
                        cgtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                                  <  1:6019221/101‑1 (MQ=255)
                        cgtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                                  >  1:2067203/1‑101 (MQ=255)
                         gtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGC.CGGCGCTCCCCCCGCGGGGCCGGGGATCCACACt                                                                                 >  1:2286003/1‑101 (MQ=37)
                         gtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                                                 <  2:4476674/101‑1 (MQ=255)
                           cgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgc                                                                               >  1:1272395/1‑101 (MQ=255)
                             cTCCCCCCGGGGGTGCGCGGACCGAAACCGGCG.GCGGGAGTTGCCCGTGATAGC.CGCCGTGGC.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgcgc                                                                             <  2:2307426/101‑1 (MQ=21)
                              tCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgcgac                                                                            >  2:4542251/1‑99 (MQ=255)
                                    cgcgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTAt                                                                      <  1:4760016/101‑1 (MQ=255)
                                    cgcgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTAt                                                                      <  1:2421494/101‑1 (MQ=255)
                                    cgcgGGGGCGTGGAGCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTAt                                                                      <  2:4252732/101‑1 (MQ=255)
                                    cgcgGGGGCGGGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTAt                                                                      <  2:6369910/101‑1 (MQ=255)
                                     gcgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATc                                                                     <  1:5860913/101‑1 (MQ=255)
                                       gggggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGc                                                                   >  2:5376205/1‑101 (MQ=255)
                                         gggCGTTGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCTCTCCCCCCGCGGGGGCGTGGATCGAAGCTGCGCCGCCTATCGCAt                                                                 <  1:6723444/101‑1 (MQ=255)
                                           gCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATcc                                                               <  2:108297/101‑1 (MQ=255)
                                           gCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATcc                                                               <  1:4884305/101‑1 (MQ=255)
                                           gCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATcc                                                               <  1:4055940/101‑1 (MQ=255)
                                            cGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCg                                                              >  2:5857493/1‑101 (MQ=255)
                                              tGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...t                                                         <  1:1791002/101‑1 (MQ=255)
                                               ggATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...Tc                                                        <  2:1537139/101‑1 (MQ=255)
                                               ggATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...Tc                                                        <  2:1126056/101‑1 (MQ=255)
                                                gATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCg                                                       <  2:3370767/101‑1 (MQ=255)
                                                gATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCg                                                       >  1:1954574/1‑101 (MQ=255)
                                                  tCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAg                                                     >  1:351732/1‑101 (MQ=255)
                                                  tCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAg                                                     <  1:1461541/101‑1 (MQ=255)
                                                    gAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAgcg                                                   >  2:1450736/1‑101 (MQ=255)
                                                     aaaCCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAgcgc                                                  >  2:6019968/1‑101 (MQ=255)
                                                     aaaCCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAgagc                                                  >  2:939758/1‑101 (MQ=255)
                                                         cTGCG.GCGTGAGTTGCCAGTCATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCcgcg                                              <  1:5624248/101‑1 (MQ=255)
                                                                 gTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGc                                      >  2:7073200/1‑101 (MQ=255)
                                                                 gTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGc                                      >  1:1683815/1‑101 (MQ=255)
                                                                  tGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCt                                     <  2:1894351/101‑1 (MQ=255)
                                                                  tGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCt                                     >  2:4783940/1‑101 (MQ=255)
                                                                   gAGTTGCCAGTTATAGA.CGCCGTGGA.CGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTc                                    <  2:5688779/101‑1 (MQ=255)
                                                                    aGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTcc                                   <  2:6626703/101‑1 (MQ=255)
                                                                        gCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTcccccc                               <  2:2391465/101‑1 (MQ=255)
                                                                        gCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTcccccc                               >  1:3455718/1‑101 (MQ=255)
                                                                        gCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTcccccc                               <  1:5046987/101‑1 (MQ=255)
                                                                        gCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTcccccc                               >  2:4153204/1‑101 (MQ=255)
                                                                           aGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCcgcg                            <  1:1101713/101‑1 (MQ=255)
                                                                             tGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCggg                          >  2:742705/1‑101 (MQ=255)
                                                                               aTAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCggggg                        >  2:5188918/1‑101 (MQ=255)
                                                                                tAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGc                       >  2:6913077/1‑101 (MQ=255)
                                                                                tAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGc                       >  2:2481263/1‑101 (MQ=255)
                                                                                  gA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCgg                     >  1:5548280/1‑100 (MQ=255)
                                                                                  gA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGt                     <  1:2985900/101‑1 (MQ=255)
                                                                                   a.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTg                    <  1:3772078/101‑1 (MQ=255)
                                                                                     cGCGGTGGA.CGTCCCTCCCCCCGTGGGGGCGGGGATCGAACCTGCGCCACCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTgg                   <  2:4070035/101‑1 (MQ=38)
                                                                                       ccatcctaggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGa                  <  1:3586215/92‑1 (MQ=9)
                                                                                       ccaacctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGa                  <  1:966489/92‑1 (MQ=9)
                                                                                       ccaacctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGa                  <  2:6134033/92‑1 (MQ=9)
                                                                                        caacctatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGAt                 <  2:4999103/93‑1 (MQ=9)
                                                                                        caacctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGAt                 >  2:2986494/9‑101 (MQ=9)
                                                                                         gTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCg               >  1:1897128/1‑101 (MQ=255)
                                                                                         aacctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGATc                >  1:4456625/8‑101 (MQ=9)
                                                                                          tGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGa              >  1:5753739/1‑101 (MQ=255)
                                                                                          acctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGATCg               <  2:2075198/95‑1 (MQ=9)
                                                                                           ccggcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCGCGGCGTCCCGAAGACGGTCTTCTGGTGCCGCTGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                 >  2:5117400/5‑101 (MQ=2)
                                                                                            ctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGATCgga             >  2:1025207/5‑99 (MQ=9)
                                                                                            ctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGaa             >  2:1366227/5‑101 (MQ=9)
                                                                                             tatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.CCTCGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa            >  2:3828038/4‑101 (MQ=9)
                                                                                             tatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.CCTCGCACGTCGCTCCCACCGCGGGGGCGTGGATCGaaa            >  2:3183540/4‑101 (MQ=9)
                                                                                             ggcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCGCGGCGTCCCGAAGACGGTCTTCTGGTGCCGCTGC.GTCGCTCCCCCCGCGGGGGCGTGGATCg               <  2:6689810/99‑1 (MQ=2)
                                                                                              gcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCGCGGCGTCCCGAAGACGGTCTTCTGGTGCCGCTGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa              <  1:1155281/100‑1 (MQ=2)
                                                                                               cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCGCGGCGTCCCGAAGACGGTCTTCTGGTGCCGCTGC.GTCGCTCCCCCCGCGGGGGCGTGGACCGaa             >  2:348098/1‑101 (MQ=2)
                                                                                                 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAt        >  1:3979374/1‑101 (MQ=255)
                                                                                                 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAt        >  2:2632018/1‑101 (MQ=255)
                                                                                                  cGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATg       <  2:1659801/101‑1 (MQ=255)
                                                                                                   gCTCCCCCCGGGGGGGCGGGGATCGAAACGGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGc      <  2:1792732/101‑1 (MQ=255)
                                                                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACTGNGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCg     <  1:3192165/101‑1 (MQ=255)
                                                                                                       cccccGCGGGGGCGTGGATCGAAACTGCGCCCCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTc  <  1:3309835/101‑1 (MQ=255)
                                                                                                       |                                                                                                        
ATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGCCGCGCA.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTC  >  NC_014034/1395739‑1395939

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG