Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,497 C→G 17.4% intergenic (+309/‑622) 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,4970C→G17.4% 125.6 ‑1.8 ‑0.1287intergenic (+309/‑622)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (19/31):  ref base (117/120):  total (136/151)
Fisher's exact test strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.35e-01

Only 100 of 307 total aligned reads displayed.

TGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAA  >  NC_014034/1395397‑1395597
                                                                                                    |                                                                                                     
 gcccccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.g                                                                                                     <  1:638658/98‑1 (MQ=255)
 gCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.g                                                                                                     >  1:3891483/1‑101 (MQ=255)
    cccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.gtgt                                                                                                  <  1:3228443/101‑1 (MQ=255)
      cccTGTCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCg                                                                                                <  2:2532249/99‑1 (MQ=255)
       gCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGc                                                                                               <  2:1376077/101‑1 (MQ=255)
         tGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTc                                                                                             >  2:3934627/1‑101 (MQ=255)
         gtgcacTGCCGGTGATCGGTGATCCGCCGCCGCCCCCCCCGCGGGGGGGTGGATCGAAAGATCGCCGCTTTGGCGGTGGGCAGGCTGGTAAC.GTGTCGCTc                                                                                             <  2:1588092/96‑1 (MQ=255)
          ggCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTcc                                                                                            >  2:4438311/1‑101 (MQ=255)
          ggCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTcc                                                                                            <  1:275460/101‑1 (MQ=255)
              cTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTcccccc                                                                                        <  2:3854401/101‑1 (MQ=255)
               tGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCcg                                                                                       <  1:4857107/101‑1 (MQ=255)
                     tgatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCggggg                                                                                 >  2:642205/1‑101 (MQ=255)
                      gatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGc                                                                                <  2:3056649/101‑1 (MQ=255)
                      gatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGc                                                                                >  2:1018746/1‑101 (MQ=255)
                         cggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTg                                                                             <  1:820749/101‑1 (MQ=255)
                          ggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTgg                                                                            <  1:1635762/101‑1 (MQ=255)
                           gctgggcgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                          <  1:2367564/92‑1 (MQ=0)
                           gctggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                          <  2:1700424/92‑1 (MQ=2)
                            tgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                          <  1:1902821/101‑1 (MQ=255)
                            tgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                          <  1:1504465/101‑1 (MQ=255)
                            ctggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                         <  1:4009993/93‑1 (MQ=0)
                            ctggacgtggTCGCTCCCACCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                         >  2:3685569/9‑101 (MQ=0)
                             tggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                        >  1:219245/8‑101 (MQ=2)
                             tggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                        >  2:2869411/8‑101 (MQ=2)
                              ggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                       <  2:2327372/95‑1 (MQ=2)
                              ggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                       >  2:4931009/7‑101 (MQ=2)
                              ggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                       >  1:3085772/7‑101 (MQ=2)
                              ggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGCGGATCGa                                                                       >  2:4506418/7‑101 (MQ=2)
                               tcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                       >  2:4267358/1‑101 (MQ=255)
                               gacgtggTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                      <  2:579059/96‑1 (MQ=2)
                                acgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                     >  1:505163/5‑101 (MQ=2)
                                 cgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                    <  1:2721238/98‑1 (MQ=2)
                                 cggggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                    <  1:2283641/98‑1 (MQ=2)
                                  gtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                                   >  1:493847/3‑101 (MQ=2)
                                  gtggTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                                   <  2:3582553/99‑1 (MQ=2)
                                  gCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                    <  1:1736651/101‑1 (MQ=255)
                                    ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACagc                                                                 >  2:2210255/1‑99 (MQ=2)
                                    ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                                                  <  1:1218691/101‑1 (MQ=255)
                                    ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                                                  <  1:2059217/101‑1 (MQ=255)
                                     gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtc                                                                 <  2:3659068/101‑1 (MQ=255)
                                       cccaccccccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCTCTCCCCCCGCGGGGGCGTGGATCGAAACAtctc                                                               <  2:4154969/97‑1 (MQ=255)
                                       cGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgcttt                                                              >  1:2560976/1‑99 (MQ=2)
                                       cGCCCCCCCCGCGGGGGCGTGGATCGAAACTACCCCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgcttt                                                              <  1:3674651/101‑3 (MQ=2)
                                        gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg                                                              <  1:2899007/101‑1 (MQ=255)
                                        gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg                                                              >  1:34010/1‑101 (MQ=255)
                                         cTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCgg                                                             >  1:3615483/1‑101 (MQ=255)
                                          tCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGt                                                            <  2:2125428/101‑1 (MQ=255)
                                          tCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGt                                                            <  2:3428540/101‑1 (MQ=255)
                                           ccccccGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACG.GAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaat                                                          <  2:322082/101‑7 (MQ=0)
                                             ccccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGcc                                                         <  2:4350643/101‑1 (MQ=255)
                                              cccGCGGGGGCGTGGATCGAAACTACGCCGGGGCTCTGGTC.AC.GGTTTGTGACGGTGTCGCACCCCCCGCGGGGGCGTGGACCGAAACAgctttcaatggg                                                       <  1:2471888/101‑10 (MQ=2)
                                                cgcgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTc                                                      <  2:1569565/101‑1 (MQ=255)
                                                 gcgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGGGCCGTCg                                                     >  2:419310/1‑101 (MQ=255)
                                                 gcgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCCAAACATCTCGGTGCCGTCg                                                     >  2:1506879/1‑101 (MQ=255)
                                                    ggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTc                                                  <  2:1107056/101‑1 (MQ=255)
                                                    ggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTc                                                  >  2:3117919/1‑101 (MQ=255)
                                                    ggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTc                                                  >  1:222972/1‑101 (MQ=255)
                                                    ggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAA.GTGTTGGTCCCCCCGCGGGGGGGGGGGTCGAAACACCTCGGTGCCCTCGCTc                                                  >  1:454161/1‑101 (MQ=255)
                                                      ggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGt                                                <  1:739848/101‑1 (MQ=255)
                                                       gCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAATATCTCGGTGCCGTCGCTCGTc                                               <  1:3727499/101‑1 (MQ=255)
                                                        cGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTcc                                              <  2:247596/101‑1 (MQ=255)
                                                         gTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCa                                             <  2:1526942/101‑1 (MQ=255)
                                                           ggATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGc                                           >  2:2635514/1‑101 (MQ=255)
                                                           ggATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCCCGTccctc                                           >  2:4713014/1‑98 (MQ=255)
                                                               cGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGc                                       >  2:1941249/1‑101 (MQ=255)
                                                               cGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGc                                       >  2:3641671/1‑101 (MQ=255)
                                                               cGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGc                                       >  2:680381/1‑101 (MQ=255)
                                                                  aaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTg                                    <  1:4768524/101‑1 (MQ=255)
                                                                    cATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCt                                  <  1:3444660/101‑1 (MQ=255)
                                                                    cATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCt                                  >  2:2980976/1‑101 (MQ=255)
                                                                          cgcTTTGGCGGTCCACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATTGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGc                            <  1:1554760/101‑1 (MQ=255)
                                                                             tttGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTcc                         <  2:4860298/101‑1 (MQ=255)
                                                                              ttGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTccc                        >  2:2384027/1‑101 (MQ=255)
                                                                              ttGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTccc                        <  1:2999629/101‑1 (MQ=255)
                                                                              ttGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTccc                        <  2:1457832/101‑1 (MQ=255)
                                                                              ttGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTccc                        <  2:787429/101‑1 (MQ=255)
                                                                              ttGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTccc                        >  2:1795288/1‑101 (MQ=255)
                                                                              ttGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCCCTccc                        >  1:4972627/1‑101 (MQ=255)
                                                                               tGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTcccc                       >  1:1677328/1‑101 (MQ=255)
                                                                               tGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTcccc                       <  2:2082466/101‑1 (MQ=255)
                                                                                 gCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTcccccc                     <  1:1388049/101‑1 (MQ=255)
                                                                                 gCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTcccccc                     <  2:2148512/101‑1 (MQ=255)
                                                                                  gtgtggataGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCcg                    <  2:401108/93‑1 (MQ=255)
                                                                                   ggTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCcgc                   >  2:2762534/1‑101 (MQ=255)
                                                                                        aCAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCggggg              >  1:351413/1‑101 (MQ=255)
                                                                                            gCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTg          <  1:4611117/101‑1 (MQ=255)
                                                                                             cTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTgg         >  1:107815/1‑101 (MQ=255)
                                                                                               tgTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGAt       <  2:1078131/100‑1 (MQ=255)
                                                                                               ggTAAC.GTGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGAt       <  1:1795421/101‑1 (MQ=255)
                                                                                                gTAAC.GTTTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATc      >  2:2247891/1‑101 (MQ=255)
                                                                                                gTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATc      <  2:5048612/101‑1 (MQ=255)
                                                                                                gTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATc      <  2:1671381/101‑1 (MQ=255)
                                                                                                gTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATc      <  1:4092647/101‑1 (MQ=255)
                                                                                                 tAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGGGGGTCg     >  2:929179/1‑101 (MQ=255)
                                                                                                 tAAC.GTGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCg     >  2:3358588/1‑101 (MQ=255)
                                                                                                  aaC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGa    >  1:3977120/1‑101 (MQ=255)
                                                                                                  aaC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGa    <  1:3342490/101‑1 (MQ=255)
                                                                                                   aC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa   >  1:433048/1‑101 (MQ=255)
                                                                                                   aC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa   <  1:2293623/101‑1 (MQ=255)
                                                                                                   aC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa   >  1:1038922/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
TGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAAC.GTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAA  >  NC_014034/1395397‑1395597

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG