Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,964 G→C 2.1% intergenic (+776/‑155) 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,9640G→C 2.1% 3.0 ‑1.7 ‑0.1237intergenic (+776/‑155)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (1/4):  ref base (127/105):  total (128/109)
Fisher's exact test strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.41e-01

Only 100 of 274 total aligned reads displayed.

TGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA  >  NC_014034/1395864‑1396064
                                                                                                      |                                                                                                       
tGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTg                                                                                                         <  2:3488435/101‑1 (MQ=255)
   gccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTc                                                                                                      <  2:4388267/101‑1 (MQ=255)
   cccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTc                                                                                                      <  2:3688301/100‑1 (MQ=255)
    ccgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGGCg                                                                                                     >  2:2622865/1‑101 (MQ=255)
     cgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGc                                                                                                    <  1:21109/101‑1 (MQ=255)
     cgccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGc                                                                                                    <  1:2010137/101‑1 (MQ=255)
      gccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCt                                                                                                   >  2:3103112/1‑101 (MQ=255)
      cccTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCt                                                                                                   <  2:376160/100‑1 (MQ=255)
        cTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcc                                                                                                 <  1:2857855/101‑1 (MQ=255)
          aTCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccc                                                                                               <  1:2928649/101‑1 (MQ=255)
          aTCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccc                                                                                               >  1:2425144/1‑101 (MQ=255)
          aTCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccc                                                                                               <  2:2878121/101‑1 (MQ=255)
           tCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTccccc                                                                                              <  1:3218710/101‑1 (MQ=255)
            cGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTcccccc                                                                                             <  1:2775648/101‑1 (MQ=255)
                 ccGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggg                                                                                        >  2:1565971/1‑101 (MQ=255)
                 ccGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggg                                                                                        >  2:719639/1‑101 (MQ=255)
                  cGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGTCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCgggg                                                                                       >  1:3065895/1‑101 (MQ=255)
                   ggTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCggggg                                                                                      >  1:56790/1‑101 (MQ=255)
                    gTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGc                                                                                     <  1:96089/101‑1 (MQ=255)
                    gTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGc                                                                                     <  2:1841730/101‑1 (MQ=255)
                    gTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCCCGGGGGc                                                                                     >  1:4415115/1‑101 (MQ=255)
                      cGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCgg                                                                                   >  2:2854719/1‑100 (MQ=255)
                      cGAGCGCCGCGCAGTCGCTCCCCCAGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGt                                                                                   <  2:2971172/101‑1 (MQ=255)
                        aGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGGgg                                                                                 >  1:1607112/1‑101 (MQ=255)
                          gaacaggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGC..TGGCTCGACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                               <  1:4626638/92‑1 (MQ=2)
                           gccgcGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                              <  2:2784360/101‑1 (MQ=255)
                           gccgcGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                              <  1:479992/101‑1 (MQ=255)
                              gcgcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           <  2:3443258/101‑1 (MQ=255)
                               ggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGC..TGGCTCGACCGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGGTCGaaa                                                                          >  2:328716/5‑101 (MQ=2)
                                gcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                         <  1:4398402/101‑1 (MQ=255)
                                 cAGTCGCTCCCACCGCGGGGGCGTGGATCGAAGCGCGGCAAGCACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCTTCGATCGaaa                                                                          <  2:308898/101‑1 (MQ=2)
                                  aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTc                                                                       <  1:1916348/101‑1 (MQ=255)
                                   gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCg                                                                      >  1:712807/1‑101 (MQ=255)
                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGc                                                                     <  2:2185279/101‑1 (MQ=255)
                                       cTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCt                                                                  >  2:464191/1‑101 (MQ=255)
                                        gccccccGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTc                                                                 <  2:1687711/100‑1 (MQ=255)
                                         ccccccGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGCTCGAAACTCGCCCTCg                                                                >  1:837314/1‑101 (MQ=255)
                                         ccccccGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCg                                                                <  1:2219949/101‑1 (MQ=255)
                                         ccccccGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCg                                                                >  2:1998617/1‑101 (MQ=255)
                                          cccccGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGa                                                               <  1:4855296/101‑1 (MQ=255)
                                               gcgGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCcgc                                                          >  2:1775563/1‑101 (MQ=255)
                                                  ggggCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGAAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCt                                                       >  2:4063597/1‑101 (MQ=255)
                                                     gCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGc                                                    <  1:1013914/101‑1 (MQ=255)
                                                       gTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCt                                                  >  1:732087/1‑101 (MQ=255)
                                                       gTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCt                                                  >  2:121496/1‑101 (MQ=255)
                                                        tGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                                                 >  1:3727840/1‑101 (MQ=255)
                                                        tGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCCCGCTTCCCTc                                                 >  2:556709/1‑101 (MQ=255)
                                                        tGGATCGAAACGAGGC..GTTCATGGCCTACGTCGGGCAGTGCACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                                                 <  2:4146794/101‑1 (MQ=255)
                                                         ggATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCa                                                >  1:150255/1‑101 (MQ=255)
                                                          gATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCAc                                               >  2:3207779/1‑101 (MQ=255)
                                                           aTCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACg                                              >  2:4415350/1‑101 (MQ=255)
                                                            tCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGc                                             <  1:1041999/101‑1 (MQ=255)
                                                              gAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                                           <  2:4235324/101‑1 (MQ=255)
                                                               aaaGGATGC..GTTCATGGCCTACGGCGGGCAGTGGACTGGTCGCTCCCCCCGTGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTg                                          <  1:4275467/101‑1 (MQ=255)
                                                               aaaCGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTg                                          >  1:1252000/1‑101 (MQ=255)
                                                                 aCGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGc.c                                       <  2:2518832/101‑1 (MQ=255)
                                                                    aTGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.cgct                                    >  2:4951606/1‑100 (MQ=255)
                                                                       c..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCg                               <  2:4836104/101‑1 (MQ=255)
                                                                          gTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGc                              >  1:4878176/1‑101 (MQ=255)
                                                                          gTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGc                              >  2:816159/1‑101 (MQ=255)
                                                                            tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTc                            >  1:1413128/1‑101 (MQ=255)
                                                                            tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTc                            >  1:637735/1‑101 (MQ=255)
                                                                            tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTc                            >  2:3013636/1‑101 (MQ=255)
                                                                             cATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcc                           >  1:1375475/1‑101 (MQ=255)
                                                                              aTGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTccc                          <  2:731998/101‑1 (MQ=255)
                                                                              aTGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATAGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTccc                          >  2:22764/1‑101 (MQ=255)
                                                                               tGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcccc                         <  1:1852969/101‑1 (MQ=255)
                                                                               tGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcccc                         >  1:3404837/1‑101 (MQ=255)
                                                                                 gCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcccccc                       >  1:2371862/1‑101 (MQ=255)
                                                                                  ccTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCcg                      >  2:2674060/1‑101 (MQ=255)
                                                                                  ccTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCcg                      <  1:4007528/101‑1 (MQ=255)
                                                                                   cTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCcgc                     >  2:4433299/1‑101 (MQ=255)
                                                                                     aCGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCCCCCCCCGCgg                   <  1:4367038/101‑1 (MQ=255)
                                                                                       gTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCgggg                 >  1:2618753/1‑101 (MQ=255)
                                                                                       gTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCgggg                 >  2:4363284/1‑101 (MQ=255)
                                                                                         cGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGGGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGc               <  1:4501849/101‑1 (MQ=255)
                                                                                         cGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGc               >  2:5036187/1‑101 (MQ=255)
                                                                                          gggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCg              >  1:1197335/1‑101 (MQ=255)
                                                                                          gggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCg              >  2:388777/1‑101 (MQ=255)
                                                                                          gggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCg              <  2:2166966/101‑1 (MQ=255)
                                                                                            gCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTg            >  1:1632106/1‑101 (MQ=255)
                                                                                            gCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTg            <  1:38598/101‑1 (MQ=255)
                                                                                             cAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTgg           <  1:2065759/101‑1 (MQ=255)
                                                                                             cAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGGgg           >  1:2957750/1‑101 (MQ=255)
                                                                                              aGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTggg          >  1:3026338/1‑100 (MQ=255)
                                                                                              aGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGa          >  1:363901/1‑101 (MQ=255)
                                                                                              aGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGa          >  1:1735434/1‑101 (MQ=255)
                                                                                                tGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATc        <  2:3671012/101‑1 (MQ=255)
                                                                                                tGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATc        >  1:3096508/1‑101 (MQ=255)
                                                                                                tGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATc        <  1:3501152/101‑1 (MQ=255)
                                                                                                tGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATc        <  1:1566101/101‑1 (MQ=255)
                                                                                                 tgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCTGCCTCGATCCGTTC.AGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATc        <  1:5002561/96‑1 (MQ=2)
                                                                                                 ggACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCATCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCg       >  2:4480835/1‑101 (MQ=255)
                                                                                                  gcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCg       >  2:1056044/5‑101 (MQ=2)
                                                                                                  gcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGCGGATCg       >  1:1170307/5‑101 (MQ=2)
                                                                                                  gcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTAAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCg       >  2:2752502/5‑101 (MQ=9)
                                                                                                    gcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTGAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGag     >  1:982807/3‑100 (MQ=2)
                                                                                                    gcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa     <  2:1367452/99‑1 (MQ=2)
                                                                                                    cTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGGTCGgaa    >  1:177566/1‑101 (MQ=255)
                                                                                                      ggTCCCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa  <  2:895575/101‑1 (MQ=255)
                                                                                                      |                                                                                                       
TGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGC..GTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA  >  NC_014034/1395864‑1396064

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