Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,953 G→C 2.6% intergenic (+765/‑166) 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,9530G→C 2.6% 5.6 ‑0.8 ‑0.1235intergenic (+765/‑166)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/4):  ref base (120/107):  total (123/112)
Fisher's exact test strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.58e-01

Only 100 of 249 total aligned reads displayed.

TGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCG  >  NC_014034/1395853‑1396052
                                                                                                      |                                                                                                   
tGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCgg                                                                                                     >  2:4391882/1‑101 (MQ=255)
 ggATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCggg                                                                                                    >  1:2185434/1‑101 (MQ=255)
     cGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGt                                                                                                <  1:3518232/101‑1 (MQ=255)
      gAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTg                                                                                               >  1:3671290/1‑101 (MQ=255)
       aaaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTgg                                                                                              <  2:398564/101‑1 (MQ=255)
       aaaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCGCCGCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTgg                                                                                              >  1:793342/1‑101 (MQ=255)
        aaCTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGa                                                                                             <  2:1020171/101‑1 (MQ=255)
           tGCGCCGCCTAACGCATCCGGTCGAGCGCCGCGCAGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTg                                                                                          <  2:793240/101‑1 (MQ=255)
            gcgcCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTgg                                                                                         <  2:2585300/101‑1 (MQ=255)
             cgccgcCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGt                                                                                        <  2:3101028/101‑1 (MQ=255)
              gccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTc                                                                                       <  2:4388267/101‑1 (MQ=255)
              cccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTc                                                                                       <  2:3688301/100‑1 (MQ=255)
               ccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCg                                                                                      >  1:4976927/1‑101 (MQ=255)
               ccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGGCg                                                                                      >  2:2622865/1‑101 (MQ=255)
                cgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGc                                                                                     <  1:2010137/101‑1 (MQ=255)
                   cTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTcc                                                                                  <  1:2857855/101‑1 (MQ=255)
                   cTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTcc                                                                                  >  1:5029718/1‑101 (MQ=255)
                     aTCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTcccc                                                                                <  2:2878121/101‑1 (MQ=255)
                     aTCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTcccc                                                                                >  2:31217/1‑101 (MQ=255)
                     aTCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTcccc                                                                                <  1:2928649/101‑1 (MQ=255)
                      tCGCATCCGGTCGAGCGCCGCGCAGTCGCTACCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGGCGCTccccc                                                                               >  1:4492324/1‑101 (MQ=255)
                       cGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTcccccc                                                                              >  1:2317295/1‑101 (MQ=255)
                         cATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCcgc                                                                            <  1:976344/101‑1 (MQ=255)
                            ccGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggg                                                                         >  1:3830483/1‑101 (MQ=255)
                            ccGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggg                                                                         >  2:719639/1‑101 (MQ=255)
                               gTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGc                                                                      >  1:1290091/1‑101 (MQ=255)
                               gTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCCCGGGGGc                                                                      >  1:4415115/1‑101 (MQ=255)
                                tCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCg                                                                     <  1:3987524/101‑1 (MQ=255)
                                  gAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGGGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCCGTGGCCTGGTCGCCCCCCCCGCGGGGGCggg                                                                   >  1:1194560/1‑99 (MQ=38)
                                  gAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTg                                                                   >  1:1687840/1‑101 (MQ=255)
                                   aGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTgg                                                                  >  2:3683644/1‑101 (MQ=255)
                                     gtctgcgcGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGCGGCCTCGATCCGGTCAGCCTCGTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTgg                                                                  <  2:5025538/97‑1 (MQ=2)
                                     cgccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                <  2:2997020/101‑1 (MQ=255)
                                      gccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                               <  2:1840048/101‑1 (MQ=255)
                                      gccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                               <  1:1140522/101‑1 (MQ=255)
                                      gccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                               <  2:2784360/101‑1 (MQ=255)
                                       ccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                              <  1:2848561/101‑1 (MQ=255)
                                       ccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                              <  1:2070468/101‑1 (MQ=255)
                                       ccgcGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                              <  1:4599193/101‑1 (MQ=255)
                                       acaggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGCTC..GACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                              <  1:2597771/94‑1 (MQ=2)
                                         gcgcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                            <  2:3443258/101‑1 (MQ=255)
                                         gcgcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                            <  2:29672/101‑1 (MQ=255)
                                          ggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGCTC..GACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGGTCGaaa                                                           >  2:328716/5‑101 (MQ=2)
                                          cgcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                           <  1:1599981/101‑1 (MQ=255)
                                           gcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                          >  2:2841088/1‑101 (MQ=255)
                                            cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                         <  2:168745/101‑1 (MQ=255)
                                            aggTCGCGCCACCCGGGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                         <  2:2258586/99‑1 (MQ=255)
                                             aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTc                                                        <  1:1916348/101‑1 (MQ=255)
                                              gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCg                                                       >  1:712807/1‑101 (MQ=255)
                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGcc                                                     >  1:1620817/1‑101 (MQ=255)
                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCt                                                   >  2:464191/1‑101 (MQ=255)
                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCCCGGGGCCG.GGGTTGGAACCCGCCCTc                                                  >  2:1350029/1‑101 (MQ=37)
                                                   tCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTc                                                  <  1:3262248/101‑1 (MQ=255)
                                                    ccccccGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGCTCGAAACTCGCCCTCg                                                 >  1:837314/1‑101 (MQ=255)
                                                     cccccGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGa                                                <  1:4855296/101‑1 (MQ=255)
                                                          gcgGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgc                                           >  1:2723576/1‑101 (MQ=255)
                                                          gcgGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgc                                           >  2:1775563/1‑101 (MQ=255)
                                                           cgGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgcg                                          >  2:1852836/1‑101 (MQ=255)
                                                             ggggCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCt                                        >  1:4152917/1‑101 (MQ=255)
                                                             ggggCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCt                                        >  2:4420384/1‑101 (MQ=255)
                                                             ggggCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGAAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCt                                        >  2:4063597/1‑101 (MQ=255)
                                                              gggCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCtt                                       <  1:2182988/101‑1 (MQ=255)
                                                                gCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGc                                     <  1:1013914/101‑1 (MQ=255)
                                                                 cGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGcc                                    >  2:2286639/1‑101 (MQ=255)
                                                                  gTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCGCGCTTGccc                                   >  2:4924543/1‑100 (MQ=255)
                                                                  gTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCt                                   >  1:732087/1‑101 (MQ=255)
                                                                  gTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCt                                   >  1:4872491/1‑101 (MQ=255)
                                                                   tGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                                  >  2:834981/1‑101 (MQ=255)
                                                                   tGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCCCGCTTCCCTc                                  >  2:556709/1‑101 (MQ=255)
                                                                   tGGATCGAAACGAGGCGTTCAT..GGCCTACGTCGGGCAGTGCACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                                  <  2:4146794/101‑1 (MQ=255)
                                                                       tCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGc                              <  1:1041999/101‑1 (MQ=255)
                                                                         gAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                            >  1:3080522/1‑101 (MQ=255)
                                                                          aaaGGATGCGTTCAT..GGCCTACGGCGGGCAGTGGACTGGTCGCTCCCCCCGTGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTg                           <  1:4275467/101‑1 (MQ=255)
                                                                             cGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCcg                        >  2:3073563/1‑101 (MQ=255)
                                                                             cGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCcg                        <  1:4381923/101‑1 (MQ=255)
                                                                                 gCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTc                    >  1:1197788/1‑101 (MQ=255)
                                                                                 gCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTc                    <  2:1446199/101‑1 (MQ=255)
                                                                                  cGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCg                   >  1:3027118/1‑101 (MQ=255)
                                                                                    ttCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCt                 <  2:2017453/101‑1 (MQ=255)
                                                                                     tCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTc                >  1:637735/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)
                                                                                      cAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcc               >  1:1375475/1‑101 (MQ=255)
                                                                                        t..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccc             >  1:3404837/1‑101 (MQ=255)
                                                                                        t..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccc             <  1:1852969/101‑1 (MQ=255)
                                                                                        t..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccc             >  1:1492346/1‑101 (MQ=255)
                                                                                        t..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccc             <  1:1766122/101‑1 (MQ=255)
                                                                                           ggCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTccccc            <  2:4269668/101‑1 (MQ=255)
                                                                                           ggCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTccccc            >  1:2956853/1‑101 (MQ=255)
                                                                                            gCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccccc           >  1:2371862/1‑101 (MQ=255)
                                                                                            gCCTACGTCGGGCAGGGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTcccccc           <  2:3459970/101‑1 (MQ=255)
                                                                                             ccTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCcg          >  2:3259117/1‑101 (MQ=255)
                                                                                             ccTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCcg          >  2:2674060/1‑101 (MQ=255)
                                                                                                aCGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCCCCCCCCGCgg       <  1:4367038/101‑1 (MQ=255)
                                                                                                 cGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCggg      >  2:2747184/1‑101 (MQ=255)
                                                                                                  gTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCgggg     >  2:4363284/1‑101 (MQ=255)
                                                                                                    cGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGc   >  2:3828170/1‑101 (MQ=255)
                                                                                                    cGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGc   >  2:5036187/1‑101 (MQ=255)
                                                                                                     gggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCg  <  2:2166966/101‑1 (MQ=255)
                                                                                                      |                                                                                                   
TGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCAT..GGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGCCGCGTCGCTCCCCCCGCGGGGGCG  >  NC_014034/1395853‑1396052

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