Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,688 T→A 5.4% intergenic (+500/‑431) 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,6880T→A 5.4% 23.9 ‑0.9 ‑0.0242intergenic (+500/‑431)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (5/8):  ref base (122/107):  total (127/115)
Fisher's exact test strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.83e-01

Only 100 of 274 total aligned reads displayed.

TGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGA  >  NC_014034/1395588‑1395788
                                                                                                     |                                                                                                     
tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGt                                                                                                       <  1:1674427/101‑1 (MQ=255)
  gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTc                                                                                                     <  1:2775524/101‑1 (MQ=255)
  gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTc                                                                                                     <  1:618053/101‑1 (MQ=255)
  gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTc                                                                                                     <  2:2920813/101‑1 (MQ=255)
  gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGCGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTc                                                                                                     >  2:4722053/1‑101 (MQ=255)
    tCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGa                                                                                                   >  2:3511641/1‑101 (MQ=255)
    tCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGa                                                                                                   >  2:1788551/1‑101 (MQ=255)
    tCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGa                                                                                                   >  2:1739953/1‑101 (MQ=255)
    tCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGa                                                                                                   >  1:1291487/1‑101 (MQ=255)
    tCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGa                                                                                                   <  1:2190555/101‑1 (MQ=255)
     cGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.g                                                                                                 >  2:3277359/1‑101 (MQ=255)
      gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.Gc                                                                                                >  2:1618242/1‑101 (MQ=255)
       aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.Gcc                                                                                               >  2:412842/1‑101 (MQ=255)
         accaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGc                                                                                             >  1:1049710/1‑101 (MQ=255)
           caccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAg                                                                                           >  2:2103009/1‑101 (MQ=255)
             ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTc                                                                                         <  1:1999324/101‑1 (MQ=255)
               ttGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGc                                                                                       >  1:3052107/1‑101 (MQ=255)
                 gCGGAGTGTGCCTTCCGCGGGGCGAAGCGGCGCTCCCCCCGCGGGGGCGTGGATCTAAACAGCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTc                                                                                     <  1:4610784/101‑1 (MQ=37)
                  cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTcc                                                                                    <  2:2206333/101‑1 (MQ=255)
                  cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTcc                                                                                    >  1:1895722/1‑101 (MQ=255)
                   ggATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTccc                                                                                   >  2:2345546/1‑101 (MQ=255)
                    gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTTGCTcccc                                                                                  >  2:269110/1‑101 (MQ=255)
                        gTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                                              <  2:909710/101‑1 (MQ=255)
                        gTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                                              <  1:1694552/101‑1 (MQ=255)
                            cTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgggg                                                                          <  2:3895058/101‑1 (MQ=255)
                            cTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgggg                                                                          >  1:2143582/1‑101 (MQ=255)
                             tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                                         >  1:2277616/1‑101 (MQ=255)
                                cgcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGt                                                                      >  1:6269/1‑101 (MQ=255)
                                cgcgGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGt                                                                      <  2:4102/101‑1 (MQ=255)
                                 tggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGt                                                                      <  2:370319/91‑1 (MQ=2)
                                 tggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGt                                                                      <  1:100425/91‑1 (MQ=2)
                                   ggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGa                                                                   <  1:2400245/101‑1 (MQ=255)
                                   ggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGa                                                                   <  1:412886/101‑1 (MQ=255)
                                   ggggCGAAGCGCCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGa                                                                   <  2:4575061/101‑1 (MQ=255)
                                    tgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                   >  2:3658365/8‑101 (MQ=0)
                                      gCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                >  1:2881889/1‑101 (MQ=255)
                                        gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                              <  1:3899687/101‑1 (MQ=255)
                                        gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                              >  2:1451597/1‑101 (MQ=255)
                                        gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCCCGGGGGCGTGGATCGaa                                                              >  2:3617230/1‑101 (MQ=255)
                                         aaGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                             >  1:3741065/1‑101 (MQ=255)
                                          tgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAGCCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTAGATCGaaa                                                             >  2:2423521/2‑101 (MQ=2)
                                           gcgcccccccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGGGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                           <  1:921802/93‑1 (MQ=255)
                                            cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                                          <  1:745632/101‑1 (MQ=255)
                                             gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgg                                                         >  1:3262063/1‑101 (MQ=255)
                                              tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGt                                                        <  1:1221657/101‑1 (MQ=255)
                                               cGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTa                                                       <  1:4538773/101‑1 (MQ=255)
                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATg                                                     >  2:1562879/1‑101 (MQ=255)
                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATg                                                     <  2:3342316/101‑1 (MQ=255)
                                                  tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                                    >  1:2558229/1‑101 (MQ=255)
                                                   ccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGa                                                   <  1:2692802/101‑1 (MQ=255)
                                                   ccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGa                                                   <  2:2277483/101‑1 (MQ=255)
                                                     ccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATc                                                 <  1:2307920/101‑1 (MQ=255)
                                                        cgcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTgg                                              <  2:4440136/101‑1 (MQ=255)
                                                         gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                                             <  2:2848268/101‑1 (MQ=255)
                                                          cgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGGTCCCCCCGCGGGTGCGTGGATCGAAACTGGTATGGATCTGgcg                                            >  1:303260/1‑101 (MQ=255)
                                                          cgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcg                                            >  1:1083018/1‑101 (MQ=255)
                                                          cgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcg                                            >  2:11217/1‑101 (MQ=255)
                                                           gggggCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcgc                                           >  2:3981513/1‑101 (MQ=255)
                                                                cGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACg                                      >  1:1103903/1‑101 (MQ=255)
                                                                 gTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                                     >  1:3340734/1‑101 (MQ=255)
                                                                 gTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                                     >  1:1524427/1‑101 (MQ=255)
                                                                  tGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                                    >  1:3713010/1‑101 (MQ=255)
                                                                   ggATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa                                   >  2:4600258/1‑101 (MQ=255)
                                                                      tCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGt                                >  1:4380702/1‑101 (MQ=255)
                                                                       ccaaaCCTCTGCGTGCAGGCCGCGG.GGAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc                               <  2:3103709/99‑1 (MQ=25)
                                                                       cGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc                               >  2:859857/1‑101 (MQ=255)
                                                                       cGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCACTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc                               >  2:2971104/1‑101 (MQ=255)
                                                                        gAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt                              >  2:945078/1‑101 (MQ=255)
                                                                        gAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt                              >  2:1822788/1‑101 (MQ=255)
                                                                        gAAACATCTGCACGTAGGCGGCGG.GCAGTGCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt                              <  2:1484410/101‑1 (MQ=255)
                                                                         aaccATCCGCATGTGGTCGGCGG.GCAGTTGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg                             <  2:1906405/98‑1 (MQ=37)
                                                                         aaaCATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg                             >  1:3525496/1‑101 (MQ=255)
                                                                         aaaCATCTGCATGTAGGCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg                             <  1:874158/101‑1 (MQ=255)
                                                                          aaCATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt                            <  1:4085980/101‑1 (MQ=255)
                                                                          aaCATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt                            >  2:2466824/1‑101 (MQ=255)
                                                                           ctaTCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtc                           <  2:713545/99‑1 (MQ=255)
                                                                           aCATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtc                           <  1:2466967/101‑1 (MQ=255)
                                                                             aTCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgt                         >  1:2441259/1‑101 (MQ=255)
                                                                              tgtccATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc                        <  2:1049631/97‑1 (MQ=255)
                                                                               cTGCATGTAGGCGGCGG.GCAGTCCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGtcgtcg                       <  2:1982267/101‑1 (MQ=255)
                                                                                tGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGc                      >  1:4704638/1‑101 (MQ=255)
                                                                                 gCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCt                     >  1:824685/1‑101 (MQ=255)
                                                                                  cATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTc                    >  2:1305446/1‑101 (MQ=255)
                                                                                     gTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcccc                 >  2:3359987/1‑101 (MQ=255)
                                                                                        gTCGGCGG.GCATTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCcg              <  1:1788660/101‑1 (MQ=255)
                                                                                        gTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCcg              <  2:2384527/101‑1 (MQ=255)
                                                                                        gTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCcg              >  1:4006729/1‑101 (MQ=255)
                                                                                          cggcgg.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCcgcg            >  1:1113122/1‑101 (MQ=255)
                                                                                          cggcgg.GCAGTGCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCcgcg            <  1:945155/101‑1 (MQ=255)
                                                                                           ggcgg.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCgg           >  1:1634383/1‑101 (MQ=255)
                                                                                           ggcgg.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCgg           <  2:2652036/101‑1 (MQ=255)
                                                                                             cgg.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCgggg         >  1:3352971/1‑101 (MQ=255)
                                                                                              gg.gCTGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCggggg        <  2:4460973/101‑1 (MQ=255)
                                                                                              gg.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCggggg        >  2:1097349/1‑101 (MQ=255)
                                                                                               g.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGc       <  1:4470103/101‑1 (MQ=255)
                                                                                               g.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGc       <  1:2168271/101‑1 (MQ=255)
                                                                                                 gCATTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCg      <  1:4289933/101‑1 (MQ=255)
                                                                                                   aGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCcgagatcggaag    >  1:3640038/1‑91 (MQ=255)
                                                                                                    gTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTgg   <  2:6264/101‑1 (MQ=255)
                                                                                                     ttCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGa  <  1:3178278/101‑1 (MQ=255)
                                                                                                     |                                                                                                     
TGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGA  >  NC_014034/1395588‑1395788

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