Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,955 C→G 2.6% intergenic (+767/‑164) 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,9550C→G 2.6% 5.4 ‑0.9 ‑0.2236intergenic (+767/‑164)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/4):  ref base (121/106):  total (124/112)
Fisher's exact test strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.80e-01

Only 100 of 250 total aligned reads displayed.

GATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGG  >  NC_014034/1395855‑1396055
                                                                                                        |                                                                                                    
gATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGc                                                                                                      <  2:4006879/101‑1 (MQ=255)
gATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GCTCAT..GGCCTACGTCGTGc                                                                                                      >  1:1841863/1‑101 (MQ=255)
  tCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAg                                                                                                    <  2:4391529/101‑1 (MQ=255)
   cGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGt                                                                                                   <  1:3518232/101‑1 (MQ=255)
    gAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTg                                                                                                  >  1:3671290/1‑101 (MQ=255)
     aaaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTgg                                                                                                 <  1:123155/101‑1 (MQ=255)
      aaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGa                                                                                                >  1:1153853/1‑101 (MQ=255)
      aaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGa                                                                                                <  2:1357382/101‑1 (MQ=255)
      aaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGa                                                                                                <  2:3327301/101‑1 (MQ=255)
      aaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGa                                                                                                >  1:1678944/1‑101 (MQ=255)
      aaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTAGa                                                                                                <  1:2884298/101‑1 (MQ=255)
        cTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACt                                                                                              <  2:2652816/101‑1 (MQ=255)
         tGCGCCGCCTAACGCATCCGGTCGAGCGCCGCGCAGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTg                                                                                             <  2:793240/101‑1 (MQ=255)
            gccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTc                                                                                          <  2:3685649/101‑1 (MQ=255)
             ccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCg                                                                                         >  1:4976927/1‑101 (MQ=255)
               gccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCt                                                                                       <  1:1627962/101‑1 (MQ=255)
               gccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCt                                                                                       >  2:3103112/1‑101 (MQ=255)
               cccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCt                                                                                       <  2:376160/100‑1 (MQ=255)
                 cTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTcc                                                                                     <  1:2857855/101‑1 (MQ=255)
                  tATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTccc                                                                                    >  2:2679476/1‑101 (MQ=255)
                    tCGCATCCGGTCGAGCGCCGCGCAGTCGCTACCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGGCGCTccccc                                                                                  >  1:4492324/1‑101 (MQ=255)
                     cGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTcccccc                                                                                 <  1:2775648/101‑1 (MQ=255)
                          ccGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggg                                                                            >  1:3830483/1‑101 (MQ=255)
                          ccGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggg                                                                            >  2:1565971/1‑101 (MQ=255)
                          ccGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggg                                                                            >  2:719639/1‑101 (MQ=255)
                           cGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GTCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCgggg                                                                           >  1:3065895/1‑101 (MQ=255)
                            ggTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggggg                                                                          >  1:56790/1‑101 (MQ=255)
                             gTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGc                                                                         >  1:1290091/1‑101 (MQ=255)
                             gTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGc                                                                         <  2:1841730/101‑1 (MQ=255)
                             gTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCCCGGGGGc                                                                         >  1:4415115/1‑101 (MQ=255)
                              tCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCg                                                                        >  2:4367195/1‑101 (MQ=255)
                              tCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCg                                                                        <  1:3987524/101‑1 (MQ=255)
                               cGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGt                                                                       <  2:2985292/101‑1 (MQ=255)
                               cGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGt                                                                       >  2:2756360/1‑101 (MQ=255)
                                gAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGGGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCCGTGGCCTGGTCGCCCCCCCCGCGGGGGCggg                                                                      >  1:1194560/1‑99 (MQ=38)
                                   gtctgcgcGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGCGGC..CTCGATCCGGTCAGCCTCGTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTgg                                                                     <  2:5025538/97‑1 (MQ=2)
                                    gccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                  <  1:1402640/101‑1 (MQ=255)
                                     ccgcGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                 <  1:4599193/101‑1 (MQ=255)
                                       catccatCGCCCCGCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                               >  2:3500882/7‑101 (MQ=37)
                                        cgcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                              <  1:1599981/101‑1 (MQ=255)
                                         gcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                             <  1:2553952/101‑1 (MQ=255)
                                          ccgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                            <  1:83273/99‑1 (MQ=255)
                                          cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                            >  2:2065611/1‑101 (MQ=255)
                                          cAGTCGCTCCCACCGCGGGGGCGTGGATCGAAGCGCGGCAAGCACAT..CGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCTTCGATCGaaa                                                              <  2:308898/101‑1 (MQ=2)
                                          atgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGC..TGGCTC..GACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGGGGGTCGAAAcg                                                            >  1:2967903/3‑100 (MQ=2)
                                           aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTc                                                           <  1:1916348/101‑1 (MQ=255)
                                           aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTc                                                           <  2:1831864/101‑1 (MQ=255)
                                            gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCg                                                          >  2:2239026/1‑101 (MQ=255)
                                              cGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGcc                                                        >  1:1997426/1‑101 (MQ=255)
                                                cTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCt                                                      >  2:464191/1‑101 (MQ=255)
                                                 tCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTc                                                     <  1:3262248/101‑1 (MQ=255)
                                                 gccccccGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTc                                                     <  2:1687711/100‑1 (MQ=255)
                                                  ccccccGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCg                                                    >  2:1998617/1‑101 (MQ=255)
                                                   cccccGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGa                                                   <  1:4855296/101‑1 (MQ=255)
                                                       cgcgGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcg                                               <  1:2183926/101‑1 (MQ=255)
                                                        gcgGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgc                                              >  1:2095565/1‑101 (MQ=255)
                                                        gcgGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgc                                              >  1:2723576/1‑101 (MQ=255)
                                                         cgGGGGAGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgcg                                             >  2:2904946/1‑101 (MQ=255)
                                                           ggggCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCt                                           >  1:4152917/1‑101 (MQ=255)
                                                           ggggCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCt                                           >  2:4420384/1‑101 (MQ=255)
                                                            gggCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCtt                                          >  2:3262068/1‑101 (MQ=255)
                                                              gCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGc                                        <  1:1013914/101‑1 (MQ=255)
                                                                gTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCGCGCTTGccc                                      >  2:4924543/1‑100 (MQ=255)
                                                                gTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCt                                      >  2:121496/1‑101 (MQ=255)
                                                                gTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCt                                      >  1:4872491/1‑101 (MQ=255)
                                                                gTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCACTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGACt                                      >  2:4328012/1‑101 (MQ=255)
                                                                 tGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                                     <  2:2869657/101‑1 (MQ=255)
                                                                 tGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                                     >  2:834981/1‑101 (MQ=255)
                                                                 tGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCCCGCTTCCCTc                                     >  2:556709/1‑101 (MQ=255)
                                                                  ggATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCa                                    >  1:150255/1‑101 (MQ=255)
                                                                   gATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCAc                                   >  2:3207779/1‑101 (MQ=255)
                                                                    aTCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACg                                  <  2:4126614/101‑1 (MQ=255)
                                                                     tCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGc                                 <  1:1041999/101‑1 (MQ=255)
                                                                      cGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGcc                                >  1:987302/1‑101 (MQ=255)
                                                                       gAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                               <  2:4235324/101‑1 (MQ=255)
                                                                       gAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                               >  2:3455376/1‑101 (MQ=255)
                                                                        aaaGGATGC..GTTCAT..GGCCTACGGCGGGCAGTGGACTGGTCGCTCCCCCCGTGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTg                              <  1:4275467/101‑1 (MQ=255)
                                                                          aCGATGG..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGcc                            <  1:719728/101‑1 (MQ=255)
                                                                          aCGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGcc                            <  2:2518832/101‑1 (MQ=255)
                                                                           cGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCcg                           <  1:4381923/101‑1 (MQ=255)
                                                                            gATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCcgc                          >  1:736143/1‑101 (MQ=255)
                                                                              tGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGt                        <  1:228689/101‑1 (MQ=255)
                                                                                   gTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTTCCGCGTCGc                     >  2:2833145/1‑101 (MQ=255)
                                                                                   gTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGc                     >  2:2594215/1‑101 (MQ=255)
                                                                                     tCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTc                   >  2:3013636/1‑101 (MQ=255)
                                                                                     tCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTc                   >  2:4467914/1‑101 (MQ=255)
                                                                                     tCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTc                   >  1:1413128/1‑101 (MQ=255)
                                                                                      cAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcc                  >  1:475488/1‑101 (MQ=255)
                                                                                       aT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATAGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTccc                 >  2:22764/1‑101 (MQ=255)
                                                                                        t..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccc                >  1:3404837/1‑101 (MQ=255)
                                                                                        t..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccc                <  1:1766122/101‑1 (MQ=255)
                                                                                           ggCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTccccc               >  1:2956853/1‑101 (MQ=255)
                                                                                            gCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccccc              >  1:2371862/1‑101 (MQ=255)
                                                                                              cTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCcgc            <  1:974746/101‑1 (MQ=255)
                                                                                               tACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCcgcg           >  1:661131/1‑101 (MQ=255)
                                                                                                aCGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCCCCCCCCGCgg          <  1:4367038/101‑1 (MQ=255)
                                                                                                    cGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGc      >  2:5036187/1‑101 (MQ=255)
                                                                                                     gggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCg     >  2:388777/1‑101 (MQ=255)
                                                                                                        cAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTgg  <  1:3092989/101‑1 (MQ=255)
                                                                                                        cAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTgg  <  1:2065759/101‑1 (MQ=255)
                                                                                                        |                                                                                                    
GATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCGTGG  >  NC_014034/1395855‑1396055

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