Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,025 C→T 18.5% intergenic (+837/‑94) 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,396,0250C→T18.5% 126.6 0.0 ‑0.4238intergenic (+837/‑94)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (21/23):  ref base (94/100):  total (115/123)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.41e-01

Only 100 of 265 total aligned reads displayed.

CGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATC  >  NC_014034/1395925‑1396125
                                                                                                    |                                                                                                       
cGAAACGATGCGTTCATGGCCTACGTCGGGGAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGcc                                                                                                         >  1:1487176/1‑101 (MQ=255)
cGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGcc                                                                                                         <  2:3010712/101‑1 (MQ=255)
cGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGcc                                                                                                         >  1:2798980/1‑101 (MQ=255)
cGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGcc                                                                                                         >  1:4255778/1‑101 (MQ=255)
 gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                                                                                                        <  1:2393083/101‑1 (MQ=255)
  aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTg                                                                                                       <  1:4131187/101‑1 (MQ=255)
    aCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGc.c                                                                                                    >  2:4380917/1‑101 (MQ=255)
     cGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.cg                                                                                                   <  1:1699028/101‑1 (MQ=255)
      gATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.cgc                                                                                                  <  1:4365902/101‑1 (MQ=255)
       aTGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.cgc..g                                                                                               >  1:4355998/1‑101 (MQ=255)
       aTGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.cgc..g                                                                                               >  1:2734490/1‑101 (MQ=255)
         gCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTc                                                                                             <  2:738869/101‑1 (MQ=255)
             tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTc                                                                                         >  1:2254699/1‑101 (MQ=255)
              cATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGCCCTCGAGCCGCGCTTGCCTCACGCCTGC.CGC..GGCGCTcc                                                                                        >  1:601125/1‑101 (MQ=255)
               aTGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTccc                                                                                       >  1:3974156/1‑101 (MQ=255)
                  gCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcccccc                                                                                    >  2:4175934/1‑101 (MQ=255)
                  gCCTACGTCGGGCAGTGGACTGGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcccccc                                                                                    <  2:3183991/101‑1 (MQ=255)
                    cTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCcgc                                                                                  >  1:2635946/1‑101 (MQ=255)
                     tACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCcgcg                                                                                 >  2:1204448/1‑101 (MQ=255)
                         tCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCggggg                                                                             <  2:3138494/101‑1 (MQ=255)
                            ggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGt                                                                          >  2:3725621/1‑101 (MQ=255)
                              cAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTgg                                                                        <  1:777501/101‑1 (MQ=255)
                                gTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                      >  1:3453939/1‑101 (MQ=255)
                                 ctgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                      >  2:3510725/7‑101 (MQ=2)
                                 ctgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                      >  2:2262128/7‑101 (MQ=9)
                                 ctgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                      >  1:1206894/7‑101 (MQ=2)
                                  ggACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                    >  1:1712442/1‑101 (MQ=255)
                                  ggACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                    <  1:3364206/101‑1 (MQ=255)
                                   gcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                    <  1:447186/97‑1 (MQ=2)
                                   gcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                    <  2:1045072/97‑1 (MQ=9)
                                       ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                >  2:2797503/1‑101 (MQ=2)
                                       ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                               >  1:370351/1‑101 (MQ=255)
                                       ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                               <  2:895378/101‑1 (MQ=255)
                                        gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                               >  1:1662818/1‑100 (MQ=2)
                                        gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                               <  2:2601028/101‑2 (MQ=2)
                                        gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                               <  1:4195105/101‑2 (MQ=2)
                                         tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgg                                                             <  1:3515971/101‑1 (MQ=255)
                                         tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgg                                                             <  1:1388125/101‑1 (MQ=255)
                                          cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtg                                                             <  1:1418259/101‑4 (MQ=2)
                                          cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGgc                                                            >  1:1106700/1‑101 (MQ=255)
                                           gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGgcg                                                           <  2:686919/101‑1 (MQ=255)
                                           ccTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtga                                                            <  1:48787/100‑5 (MQ=2)
                                            ccccccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGgcgc                                                          <  2:3609774/99‑1 (MQ=255)
                                            cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgat                                                           <  1:2511230/101‑6 (MQ=2)
                                             tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCCGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatg                                                          >  1:2160933/1‑95 (MQ=2)
                                             tCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCa                                                         >  1:2622897/1‑101 (MQ=255)
                                             cccccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCa                                                         <  1:3791005/100‑1 (MQ=255)
                                                ccccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatgacc                                                       <  1:3584246/101‑10 (MQ=2)
                                                ccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCt                                                      <  2:2075775/101‑1 (MQ=255)
                                                ccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCt                                                      <  2:530923/101‑1 (MQ=255)
                                                  ccGGGGGGGCGTGGAACGAACCTCGCCCTCGATCCGCAGTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGa                                                    <  2:1076025/101‑1 (MQ=37)
                                                     cgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGGAACAGGCGCATCTGATGt                                                 >  2:2109647/1‑101 (MQ=255)
                                                         ggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTatat                                             >  2:640299/1‑101 (MQ=255)
                                                             tGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGccc                                         >  2:1755252/1‑101 (MQ=255)
                                                             tGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGccc                                         <  1:4176195/101‑1 (MQ=255)
                                                             tGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGccc                                         <  2:1407604/101‑1 (MQ=255)
                                                              ggATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCg                                        <  2:1857078/101‑1 (MQ=255)
                                                              ggATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCACCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCg                                        <  2:3853007/101‑1 (MQ=255)
                                                               gATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCgg                                       >  2:1633484/1‑101 (MQ=255)
                                                                aTCGAAACTCGCCCTCGATCCGCGCTTGCCTCAGGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGt                                      <  1:3858602/101‑1 (MQ=255)
                                                                aTCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGt                                      >  2:2323381/1‑101 (MQ=255)
                                                                   gAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTgag                                   <  2:3282418/101‑1 (MQ=255)
                                                                    aaaCTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTgaga                                  >  1:1050816/1‑101 (MQ=255)
                                                                    aaaCTCGCCCTCGATCCGCCCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTgaga                                  <  2:1417291/101‑1 (MQ=255)
                                                                     aaCTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGAc                                 <  2:1692564/101‑1 (MQ=255)
                                                                      aCTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACg                                <  1:3092080/101‑1 (MQ=255)
                                                                       cTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACgg                               >  1:3237612/1‑101 (MQ=255)
                                                                         cGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCt                             <  1:1590979/101‑1 (MQ=255)
                                                                         cGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCt                             >  2:1140814/1‑101 (MQ=255)
                                                                           cccTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGt                           <  1:2333836/101‑1 (MQ=255)
                                                                           cccTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGt                           <  2:2419983/101‑1 (MQ=255)
                                                                            ccTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTc                          <  2:472324/101‑1 (MQ=255)
                                                                            ccTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTc                          <  2:1360472/101‑1 (MQ=255)
                                                                            ccTCGATCCGCGCTTGCCCCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTc                          <  1:2057623/101‑1 (MQ=255)
                                                                              tCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGc                        >  1:91717/1‑101 (MQ=255)
                                                                              tCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGc                        >  1:4097081/1‑101 (MQ=255)
                                                                                gATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTc                      <  1:2446333/101‑1 (MQ=255)
                                                                                 aTCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcc                     >  1:2944171/1‑101 (MQ=255)
                                                                                  tCCGCGGTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTccc                    <  2:1988359/101‑1 (MQ=255)
                                                                                   ccGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccc                   >  2:1665769/1‑101 (MQ=255)
                                                                                    cgcgCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTccccc                  <  2:4230227/101‑1 (MQ=255)
                                                                                     gcgcTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc                 <  1:1642980/101‑1 (MQ=255)
                                                                                     gcgcTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc                 <  1:2659199/101‑1 (MQ=255)
                                                                                     gcgcTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc                 <  1:402990/101‑1 (MQ=255)
                                                                                      cgcTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTccccccc                >  2:1949959/1‑100 (MQ=255)
                                                                                       gcTTGCCTCCCGCCTGC.CGC..GTCGCTCCCCCCGCGGGCGCGTGGATCGAAACAGGCGCATCTGATGTGTATGCCCGGTGAGACGGCTGTCGCTCCCCCcgc               <  2:2179257/101‑1 (MQ=255)
                                                                                       gcTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGCGAGACGGCTGGCGCTcccccccc               >  1:2873991/1‑99 (MQ=255)
                                                                                         ttGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCgg             <  2:2254554/101‑1 (MQ=255)
                                                                                          tGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCggg            >  1:191382/1‑101 (MQ=255)
                                                                                           gCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCgggg           >  1:1252080/1‑101 (MQ=255)
                                                                                            ccTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCggggg          <  2:2414137/101‑1 (MQ=255)
                                                                                             cgcACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGc         <  2:3966894/99‑1 (MQ=255)
                                                                                             cTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGc         >  1:673173/1‑101 (MQ=255)
                                                                                               cACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGCGCATCTGATGTATATGCCCGGGGAGAGGTCTGTCACTCCCCCCGCGGGGGCGt       >  2:1501019/1‑101 (MQ=37)
                                                                                               cACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGt       <  2:1252913/101‑1 (MQ=255)
                                                                                               cACGCCGGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGt       <  2:2756891/101‑1 (MQ=255)
                                                                                                   ccTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGAt   >  1:821271/1‑101 (MQ=255)
                                                                                                   ccTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGAt   <  1:1320775/101‑1 (MQ=255)
                                                                                                    cTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATc  <  1:2530203/101‑1 (MQ=255)
                                                                                                    cTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATc  <  2:2096846/101‑1 (MQ=255)
                                                                                                    |                                                                                                       
CGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGATC  >  NC_014034/1395925‑1396125

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