Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,494 T→G 18.4% intergenic (+306/‑625) 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,4940T→G18.4% 190.5 ‑4.0 ‑1.7394intergenic (+306/‑625)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (23/50):  ref base (150/169):  total (174/220)
Fisher's exact test strand distribution p-value = 1.85e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.06e-02

Only 100 of 418 total aligned reads displayed.

AACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCG  >  NC_014034/1395394‑1395594
                                                                                                       |                                                                                                    
aaCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGt                                                                                                      <  2:4806765/101‑1 (MQ=255)
 aCTGCCCCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTa                                                                                                     <  2:2127539/101‑1 (MQ=255)
 aCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTa                                                                                                     <  1:4148418/101‑1 (MQ=255)
 aCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTa                                                                                                     >  1:2972067/1‑101 (MQ=255)
  cTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTaa                                                                                                    <  1:6450452/101‑1 (MQ=255)
   tGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAAc                                                                                                   >  2:5087139/1‑101 (MQ=255)
      aCCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGCCAGGCTGGTAACgcg                                                                                                >  2:1351670/1‑99 (MQ=255)
       cccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACgtgt                                                                                               <  2:1455811/101‑1 (MQ=255)
        ccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTc                                                                                              <  2:1641019/101‑1 (MQ=255)
         cGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCg                                                                                             <  1:6771000/101‑1 (MQ=255)
         cGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCg                                                                                             >  2:5344640/1‑101 (MQ=255)
          gCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGc                                                                                            <  2:5474588/101‑1 (MQ=255)
          gCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGGAACGTGTCGc                                                                                            >  1:2661362/1‑101 (MQ=255)
              gCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTccc                                                                                        <  1:4675789/101‑1 (MQ=255)
              gCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTccc                                                                                        <  1:100025/101‑1 (MQ=255)
                ccTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTccccc                                                                                      <  2:3536974/101‑1 (MQ=255)
                 cTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccccc                                                                                     >  1:3415308/1‑101 (MQ=255)
                  tGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCcg                                                                                    <  2:4954754/101‑1 (MQ=255)
                   gcgggtgatcggtgctcCGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCcgc                                                                                   <  2:2023382/98‑1 (MQ=255)
                   gCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCcgc                                                                                   >  1:4737716/1‑101 (MQ=255)
                       gtgatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCgggg                                                                               <  1:557146/101‑1 (MQ=255)
                         gatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGc                                                                             >  2:2289707/1‑101 (MQ=255)
                            cggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTg                                                                          <  2:3739982/101‑1 (MQ=255)
                             ggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTgg                                                                         <  1:3338880/101‑1 (MQ=255)
                             ggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTgg                                                                         <  2:1847939/101‑1 (MQ=255)
                              gtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGa                                                                        >  2:670421/1‑101 (MQ=255)
                               tgccgctgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGG.CGA.ATACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGa                                                                        <  1:5367378/92‑1 (MQ=2)
                               ctggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..ACACCGGG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                      <  1:1907630/93‑1 (MQ=0)
                                tggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..ACACCGGG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                     <  2:6572241/94‑1 (MQ=2)
                                gatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                      <  2:1203036/101‑1 (MQ=255)
                                 atcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                     <  1:2376017/101‑1 (MQ=255)
                                 atcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                     >  2:3250141/1‑101 (MQ=255)
                                 atcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                     <  2:2847743/101‑1 (MQ=255)
                                 atcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGGGGATCg                                                                     >  1:2076237/1‑101 (MQ=255)
                                  tcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCTCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                    <  2:5585230/101‑1 (MQ=255)
                                   cAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                   <  1:1985748/101‑1 (MQ=255)
                                    aGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                  <  1:3143625/101‑1 (MQ=255)
                                     gtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..ACACCGGG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGCAACa                                                                >  2:305981/3‑101 (MQ=2)
                                     gtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..ACACCGGG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                                <  2:4274681/99‑1 (MQ=2)
                                     gCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                 <  1:6965863/101‑1 (MQ=255)
                                      cGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                                >  2:1067385/1‑101 (MQ=255)
                                       ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGGTCCCCCCGCGGGGGCGGGGATCGAAACAt                                                               >  1:1488535/1‑101 (MQ=255)
                                        gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..ACACCGGG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGGGGGTCGAAACAcct                                                             >  2:7043087/1‑101 (MQ=0)
                                        gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtc                                                              >  1:3642701/1‑101 (MQ=255)
                                         tCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..ACACCGGG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctt                                                            <  1:2848035/101‑2 (MQ=2)
                                         tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtct                                                             >  1:5659688/1‑101 (MQ=255)
                                          cGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtctc                                                            >  2:19991/1‑101 (MQ=255)
                                          cGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtctc                                                            <  1:6217824/101‑1 (MQ=255)
                                           gCTCCCCCCGCGGGGGCGTGGATCGAAACT..ACACCGGG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttc                                                          <  1:740277/101‑4 (MQ=2)
                                           gCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg                                                           >  1:1161415/1‑101 (MQ=255)
                                           gCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg                                                           >  1:3839374/1‑101 (MQ=255)
                                            cgccccgcGCGGGGGCGTGGATCGCACCA..TCGCCGCT.TTGGCCGTCGACAGGCTCGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCgg                                                          <  2:65015/99‑1 (MQ=37)
                                            cTCCCCCCGCGGGGGCGTGGATCGAAACT..ACACCGGG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttca                                                         <  1:5822231/101‑5 (MQ=2)
                                            cTCCCCCCGCGGGGGCGTGGATCGAAACT..ACACCGGG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttca                                                         <  1:3269440/101‑5 (MQ=2)
                                            cTCCCCCCGCGGGGGCGTGGATCGAAACT..ACACCGGG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttca                                                         <  2:301044/101‑5 (MQ=2)
                                             gccccccGCGGGGGCGTGGATCGAAACT..ACACCGGG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaa                                                        <  1:1825736/100‑6 (MQ=0)
                                              ccccccGCGGGGGCGTGGATCGAAACT..ACACCGGG.GCTCTGGTC.ACGGTTTGGCACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaat                                                       <  1:861952/101‑7 (MQ=0)
                                               cccccGCGGGGGCGTGGATCGAAACT..ACACCGGG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaatg                                                      <  2:335451/101‑8 (MQ=2)
                                                ccccGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGcc                                                      <  2:71059/101‑1 (MQ=255)
                                                     cgGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGc                                                 >  1:4364972/1‑101 (MQ=255)
                                                     cgGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGc                                                 >  2:4425763/1‑101 (MQ=255)
                                                      gggggCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCt                                                <  1:1727795/101‑1 (MQ=255)
                                                            gTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCa                                          >  2:6078652/1‑101 (MQ=255)
                                                             tGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTTCCATCGCTCGTCCAg                                         >  2:591367/1‑101 (MQ=255)
                                                              ggATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGCCGCTCGTCCAGc                                        >  2:3910154/1‑101 (MQ=255)
                                                                   gAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGcc                                   <  1:6728613/101‑1 (MQ=255)
                                                                    aaaCA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCt                                  <  1:1241700/101‑1 (MQ=255)
                                                                    aaaCA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCt                                  <  1:1576947/101‑1 (MQ=255)
                                                                     aaCA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTg                                 >  1:6098403/1‑101 (MQ=255)
                                                                     aaCA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTg                                 >  2:5682041/1‑101 (MQ=255)
                                                                     aaCA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTg                                 <  1:2240420/101‑1 (MQ=255)
                                                                     aaCA..GCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTg                                 <  1:559098/101‑1 (MQ=255)
                                                                        a..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtc                              <  2:1891547/101‑1 (MQ=255)
                                                                          agcgccgcT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtc                              <  2:4670891/99‑1 (MQ=255)
                                                                            cgccgcT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgt                            <  1:1409203/101‑1 (MQ=255)
                                                                             gccgcT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATTGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgtc                           <  1:3918751/101‑1 (MQ=255)
                                                                             gccgcT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGGtcgtc                           >  2:2564652/1‑101 (MQ=255)
                                                                             gccgcT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGCGCCGTCGCTCGTCCAGCCGGCCTGCtcgtc                           >  1:483669/1‑101 (MQ=255)
                                                                                gcT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGcc                        >  2:3440429/1‑100 (MQ=255)
                                                                                gcT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCt                        <  2:6306192/101‑1 (MQ=255)
                                                                                gcT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCt                        <  2:6254381/101‑1 (MQ=255)
                                                                                gcT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCt                        >  1:2473712/1‑101 (MQ=255)
                                                                                 cT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTc                       >  1:6594888/1‑101 (MQ=255)
                                                                                 cT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTc                       <  2:2070953/101‑1 (MQ=255)
                                                                                  t.ttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTAGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTcc                      >  1:5100944/1‑101 (MQ=255)
                                                                                      ggCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTccccc                   >  2:2457311/1‑101 (MQ=255)
                                                                                         ggTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCcgc                <  1:1780708/101‑1 (MQ=255)
                                                                                          gTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTGCCGCCGGCCTGCTCGTCGCTCCCCCcccg               >  2:6163043/1‑101 (MQ=255)
                                                                                            cGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGGTCCCCCCGCggg             >  1:7080735/1‑101 (MQ=255)
                                                                                               cAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGc          <  2:1753228/101‑1 (MQ=255)
                                                                                                 ggCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGt        >  1:4558014/1‑101 (MQ=255)
                                                                                                  gCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTg       >  1:6382314/1‑101 (MQ=255)
                                                                                                  gCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTg       >  1:4622481/1‑101 (MQ=255)
                                                                                                   cTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTgg      >  2:6997094/1‑101 (MQ=255)
                                                                                                    tGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGa     <  1:2766275/101‑1 (MQ=255)
                                                                                                    tGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGa     >  2:4895176/1‑101 (MQ=255)
                                                                                                    tGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCGGCTCGTCGCTCCCCCCGCGGGGGCGTGGa     >  2:54058/1‑101 (MQ=255)
                                                                                                     ggTACCGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGAt    <  2:4060107/101‑1 (MQ=255)
                                                                                                     ggTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGAt    <  2:3156785/101‑1 (MQ=255)
                                                                                                      gTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATc   >  2:2687919/1‑101 (MQ=255)
                                                                                                       |                                                                                                    
AACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..TCGCCGCT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCG  >  NC_014034/1395394‑1395594

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