Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,698 Δ1 bp 6.2% intergenic (+510/‑421) 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,6980A→. 6.2% 47.5 0.0 ‑0.0304intergenic (+510/‑421)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (9/10):  ref base (135/150):  total (144/160)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.39e-01

Only 100 of 324 total aligned reads displayed.

CCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTG  >  NC_014034/1395598‑1395798
                                                                                                      |                                                                                                       
 cgccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAg                                                                                                        <  1:4537754/99‑1 (MQ=255)
        cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTcc                                                                                                 >  1:2100483/1‑101 (MQ=255)
        cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTcc                                                                                                 >  2:4567470/1‑101 (MQ=255)
          gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTcccc                                                                                               <  2:279425/101‑1 (MQ=255)
            ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTcccccc                                                                                             >  1:2191126/1‑101 (MQ=255)
             tGTCTCTTGCGGGGGCCTAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                                                            <  2:210493/101‑1 (MQ=37)
             tGTCGCTTCCGCGGGGCGAGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                                                            <  1:4318643/101‑1 (MQ=255)
             tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                                                            <  2:419156/101‑1 (MQ=255)
              gTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                                                           >  2:5311451/1‑101 (MQ=255)
              gTCACTTTCGC.GGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgcg                                                                                          >  2:398376/1‑101 (MQ=255)
               tCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgcg                                                                                          <  1:3976555/101‑1 (MQ=255)
                cccTTTGGCGGGGCGAAGCGTCGCTCCCCCGGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgg                                                                                         <  2:4745535/99‑1 (MQ=255)
                 ccTTTCGCGGGGGGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggg                                                                                        <  1:449680/100‑1 (MQ=255)
                   tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                                                      <  1:3541201/101‑1 (MQ=255)
                   tttCGCGGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGGAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                                                      <  1:2447720/101‑1 (MQ=255)
                    ttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGGAGTCGGCGG.CCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGc                                                                                     <  1:5299461/101‑1 (MQ=255)
                     tCGGGGGCCGAGGCGCCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCg                                                                                    <  2:4949166/101‑1 (MQ=38)
                     tCGCGGGGCGAAGCGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCg                                                                                    <  1:4030392/101‑1 (MQ=255)
                      cgcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCACTTCGA.GCCGCAGTCGCTCCGCCCGCGGGGGCgg                                                                                   >  1:4973124/1‑100 (MQ=255)
                       gcgGGGCGAGGCGTCGCGCCCCCCGCGGGGGCGGGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTg                                                                                  <  1:3193946/101‑1 (MQ=255)
                       gcgGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTg                                                                                  <  2:66279/101‑1 (MQ=255)
                         gtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                                                 <  1:268315/93‑1 (MQ=2)
                         gtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                                                 <  1:1855369/93‑1 (MQ=2)
                          gggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                               >  2:4215510/1‑101 (MQ=255)
                          gggCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                               <  2:3592893/101‑1 (MQ=255)
                           ggCGAGGCGTCGCTCCCCCCGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                              <  1:1019576/101‑1 (MQ=255)
                           ggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                              >  2:694992/1‑101 (MQ=255)
                            gCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                             <  2:5727585/101‑1 (MQ=255)
                            gCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                             <  1:3437958/101‑1 (MQ=255)
                             cGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                            >  2:1067220/1‑101 (MQ=255)
                             cGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                            <  2:2859818/101‑1 (MQ=255)
                              gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGGTCGaa                                                                           >  1:1210895/1‑101 (MQ=255)
                              gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           <  2:34782/101‑1 (MQ=255)
                              gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           >  2:5449150/1‑101 (MQ=255)
                              gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           >  1:2914021/1‑101 (MQ=255)
                               aaccGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                          <  2:611790/98‑1 (MQ=255)
                                tgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGGTCGaaa                                                                          >  2:190505/2‑101 (MQ=2)
                                tgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                          >  1:1659929/2‑101 (MQ=2)
                                  cgcccccccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTCCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                      <  1:5434839/94‑1 (MQ=255)
                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTa                                                                   >  1:5333049/1‑101 (MQ=255)
                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTa                                                                   >  2:4349504/1‑101 (MQ=255)
                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTa                                                                   >  2:5346040/1‑101 (MQ=255)
                                        tCCCCCGGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                <  2:3821302/101‑1 (MQ=255)
                                        cccccccGCGGGGGCGGGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                <  2:3946707/100‑1 (MQ=255)
                                          cccccGGGGGGGCGGGGATCGAAACATCTGCATGTGGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGAt                                                              <  1:5683637/101‑1 (MQ=255)
                                            tccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                            <  1:1528175/100‑1 (MQ=255)
                                            cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                            <  2:2204460/101‑1 (MQ=255)
                                            cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                            <  1:1788474/101‑1 (MQ=255)
                                             ccGCGGGGGCGTGGATCAAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTg                                                           <  1:4187352/101‑1 (MQ=255)
                                               gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGgc                                                         >  1:5070293/1‑101 (MQ=255)
                                                 gggggCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGgcgc                                                       >  2:491180/1‑101 (MQ=255)
                                                    ggCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCaaa                                                    >  2:1243562/1‑101 (MQ=255)
                                                    ggCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCaaa                                                    <  1:3776250/101‑1 (MQ=255)
                                                     gCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAAc                                                   >  1:4829436/1‑101 (MQ=255)
                                                     gCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAAc                                                   >  1:5721522/1‑101 (MQ=255)
                                                       gTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACgg                                                 >  2:3194808/1‑101 (MQ=255)
                                                        tGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..a                                              >  2:1574342/1‑101 (MQ=255)
                                                          gATCGAGACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..aaa                                            <  1:945727/101‑1 (MQ=255)
                                                          gATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..aaa                                            <  1:4138823/101‑1 (MQ=255)
                                                           aTCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAg                                           >  1:624109/1‑101 (MQ=255)
                                                           aTCAAAACCTCGGCATATAGTCGGCGG.GCATTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAg                                           <  2:3332188/101‑1 (MQ=37)
                                                            tCGAATCTTCTGCGTGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGt                                          <  2:1975368/101‑1 (MQ=255)
                                                            tCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGt                                          <  2:827303/101‑1 (MQ=255)
                                                             cgcaaCATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTc                                         >  1:3297113/4‑101 (MQ=255)
                                                             cGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTc                                         <  2:413936/101‑1 (MQ=255)
                                                              gAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCt                                        <  1:5218408/101‑1 (MQ=255)
                                                              gAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGGATGGATCTGGCGCAAACGG..AAAGTCt                                        >  2:2060028/1‑101 (MQ=255)
                                                                aaCATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgt                                      >  2:4846339/1‑101 (MQ=255)
                                                                aaCATCTGCATGTAGGCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgt                                      <  1:1067295/101‑1 (MQ=255)
                                                                    tCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcgtc                                  <  2:4495295/101‑1 (MQ=255)
                                                                    tCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcgtc                                  >  2:4044563/1‑101 (MQ=255)
                                                                      tGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACG.GGTATGGAGCTGGCGCAAACGG..AAAGTCTGTCGTCGc                                >  2:1361705/1‑101 (MQ=255)
                                                                       gCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGcc                               >  1:3102790/1‑100 (MQ=255)
                                                                            tAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTccccc                          >  1:3696244/1‑101 (MQ=255)
                                                                             ggccggcgg.GCAGTTGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTcccccc                         <  1:3030287/98‑1 (MQ=255)
                                                                             aggcggcgg.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTcccccc                         <  1:3837876/98‑1 (MQ=255)
                                                                             aGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTcccccc                         >  2:4107372/1‑101 (MQ=255)
                                                                              gTCGGGGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCcg                        <  2:2431838/101‑1 (MQ=255)
                                                                               tgggcgg.GCATTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCcgc                       <  2:5201928/99‑1 (MQ=255)
                                                                                cggcgg.GCAGTGCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCcgcg                      <  1:5530866/101‑1 (MQ=255)
                                                                                 ggcgg.CCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCgg                     <  2:5043636/101‑1 (MQ=255)
                                                                                  gcgg.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCggg                    >  1:1378958/1‑101 (MQ=255)
                                                                                    gg.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCggggg                  >  2:560231/1‑101 (MQ=255)
                                                                                     g.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGc                 >  1:5397658/1‑101 (MQ=255)
                                                                                       gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGGTCGAAAAC.GGGATGGGTCTTGGCCAAACGG..AAAGGCCGTCGGCCCCCCCCCCCCcgggggg                >  1:2013780/1‑99 (MQ=12)
                                                                                          gTTGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTgg             <  1:2551897/101‑1 (MQ=255)
                                                                                          gTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTgg             >  1:196773/1‑101 (MQ=255)
                                                                                              caagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTgg             >  2:1821220/8‑101 (MQ=2)
                                                                                               aagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGa            <  1:3569634/95‑1 (MQ=2)
                                                                                                  ccGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa      >  2:343802/1‑101 (MQ=255)
                                                                                                  attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCg         >  2:3052593/4‑101 (MQ=2)
                                                                                                  attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCg         <  2:3171379/98‑1 (MQ=2)
                                                                                                  attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCg         <  2:5237643/98‑1 (MQ=2)
                                                                                                   ttcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGa        <  2:2660514/99‑1 (MQ=2)
                                                                                                    tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       <  2:3757600/100‑1 (MQ=2)
                                                                                                    gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc    >  2:3290382/1‑101 (MQ=255)
                                                                                                     cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa      <  1:875239/101‑1 (MQ=2)
                                                                                                      aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaaa     >  1:485874/1‑100 (MQ=2)
                                                                                                      aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc     >  1:2374453/1‑101 (MQ=2)
                                                                                                      aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAACTg  >  2:4093635/1‑101 (MQ=255)
                                                                                                      |                                                                                                       
CCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTG  >  NC_014034/1395598‑1395798

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG