Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,823,740 +C 100% intergenic (+25/‑148) RCAP_rcc01687 → / → RCAP_rcc01688 HK97 family phage major capsid protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,823,7401.→C100.0% 62.7 163intergenic (+25/‑148)RCAP_rcc01687/RCAP_rcc01688HK97 family phage major capsid protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (92/50):  ref base (0/21):  total (92/71)
Fisher's exact test strand distribution p-value = 3.76e-09
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.78e-03

Only 100 of 197 total aligned reads displayed.

TATGCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.A.CCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGC  >  NC_014034/1823641‑1823832
                                                                                                     |                                                                                            
 aTGCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCC.A.cc                                                                                            >  2:2180911/1‑101 (MQ=255)
 aTGCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.A.cc                                                                                            <  1:761590/101‑1 (MQ=255)
   gCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.A.cccc                                                                                          <  1:4096377/101‑2 (MQ=255)
   gCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.A.cccc                                                                                          <  1:2284758/101‑2 (MQ=255)
    cGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.A.ccccg                                                                                         >  2:3818134/1‑99 (MQ=255)
    cGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.A.ccccg                                                                                         <  2:1903127/101‑3 (MQ=255)
     naNNANNCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.A.ccccgn                                                                                        <  1:1906454/100‑4 (MQ=255)
     gAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.A.ccccgg                                                                                        >  2:567850/1‑98 (MQ=255)
     gAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.A.ccccgg                                                                                        >  2:4305762/1‑98 (MQ=255)
      aGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.A.ccccggc                                                                                       >  2:4273405/1‑97 (MQ=255)
      aGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.A.ccccggc                                                                                       <  1:2107648/101‑5 (MQ=255)
       gCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.A.ccccggct                                                                                      >  2:2163568/1‑96 (MQ=255)
            cGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTgccgc                                                                                  >  1:1092621/1‑101 (MQ=255)
             gggTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTgccgcc                                                                                 >  1:1344084/1‑101 (MQ=255)
             gggGGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTgccgcc                                                                                 <  1:3200571/101‑1 (MQ=255)
               gTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCgg                                                                               >  1:3126773/1‑101 (MQ=255)
               gTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCgg                                                                               >  2:23358/1‑101 (MQ=255)
                 gggCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGtt                                                                             <  1:4453007/101‑1 (MQ=255)
                 gggCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGtt                                                                             <  1:322163/101‑1 (MQ=255)
                    cggcgATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGc                                                                          >  1:1059961/1‑101 (MQ=255)
                     ggcgATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGcc                                                                         >  1:906717/1‑101 (MQ=255)
                     ggcgATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGcc                                                                         >  1:486132/1‑101 (MQ=255)
                     ggcgATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGcc                                                                         >  2:3505764/1‑101 (MQ=255)
                       cgATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCgg                                                                       <  1:3818419/101‑1 (MQ=255)
                        gATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGt                                                                      >  1:296006/1‑101 (MQ=255)
                        gATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGt                                                                      >  2:3804617/1‑101 (MQ=255)
                        gATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGAA.ACCCCGGCTGCCGCCGGTTCGCCGGt                                                                      >  1:3703244/1‑101 (MQ=255)
                         aTGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTc                                                                     <  2:2606297/101‑1 (MQ=255)
                          tGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTcg                                                                    <  2:4639916/101‑1 (MQ=255)
                            tCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGGCGGGcggg                                                                  >  2:3593572/1‑101 (MQ=255)
                            tCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTcggg                                                                  <  2:4538527/101‑1 (MQ=255)
                              aGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGTCA.ACCCCGGCTGCCGCCGGTTCGCCGGTcgggcg                                                                >  1:1966890/1‑101 (MQ=255)
                              aGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTcgggcg                                                                <  2:2269787/101‑1 (MQ=255)
                               gCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTcgggcgg                                                               >  2:3149198/1‑101 (MQ=255)
                               gCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTcgggcgg                                                               <  2:4014852/101‑1 (MQ=255)
                                cGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTcgggcggg                                                              >  2:1146035/1‑101 (MQ=255)
                                cGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTcgggcggg                                                              >  1:4387853/1‑101 (MQ=255)
                                 gATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCgggcgggc                                                             <  2:3370615/101‑1 (MQ=255)
                                 gATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCgggcgggc                                                             >  1:3069600/1‑101 (MQ=255)
                                  aTTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGggcgggcg                                                            <  1:2524037/101‑1 (MQ=255)
                                   tttCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGgcgggcgc                                                           >  1:1807454/1‑101 (MQ=255)
                                   tttCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGgcgggcgc                                                           >  1:2919881/1‑101 (MQ=255)
                                    ttCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGgcgggcgcg                                                          >  2:2244217/1‑101 (MQ=255)
                                      cgccgcGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGgcgggcgcggg                                                        >  1:1423381/1‑101 (MQ=255)
                                      cgccgcGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGgcgggcgcggg                                                        >  2:4352849/1‑101 (MQ=255)
                                       gccgcGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGgcgggcgcgggc                                                       >  2:1038902/1‑101 (MQ=255)
                                         cgctATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCtt                                                     <  1:3603489/101‑1 (MQ=255)
                                         cgcgATCAAGCTGTGGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCtt                                                     <  2:3613512/101‑1 (MQ=255)
                                          gcgATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTc                                                    <  1:2250418/101‑1 (MQ=255)
                                           cgATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTcc                                                   >  1:2431666/1‑101 (MQ=255)
                                           cgATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTcc                                                   <  2:2862855/101‑1 (MQ=255)
                                            gATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCg                                                  >  2:2754416/1‑101 (MQ=255)
                                            gATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGAGCA.ACCCCGGCTGCCGCCGGTTCGCCGGCCGGGCGGGAGCCGGCGTCCg                                                  >  2:1348609/1‑101 (MQ=255)
                                              tCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGc                                                >  2:4181066/1‑101 (MQ=255)
                                              tCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGc                                                >  1:163087/1‑101 (MQ=255)
                                                aaGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCACC.CCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCg                                              <  1:2765483/101‑1 (MQ=255)
                                                aaGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCg                                              >  2:3094592/1‑101 (MQ=255)
                                                  gCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCt                                            <  2:1946393/101‑1 (MQ=255)
                                                    tGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGc                                          >  2:3702074/1‑101 (MQ=255)
                                                    tGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGc                                          >  1:4641359/1‑101 (MQ=255)
                                                     gTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGcc                                         >  1:2137507/1‑101 (MQ=255)
                                                       tGAAGTTTGCCGCCTCGTAACGGGGGGTGAAGGGGCGGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTa                                       <  1:824373/101‑1 (MQ=255)
                                                         aaGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgc                                     <  1:3779448/101‑1 (MQ=255)
                                                         aaGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgc                                     >  1:4199079/1‑101 (MQ=255)
                                                          attttGCCGCCTCGTAAGGGGCGGTGAAGGGGCCGGGGGCACC.CCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgct                                    <  1:23363/99‑1 (MQ=255)
                                                          aGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgct                                    <  2:1807328/101‑1 (MQ=255)
                                                          aGTTTCCCGCCTCGTAACGGGCGGTGAAGGGGCGGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgct                                    <  2:2925330/101‑1 (MQ=255)
                                                           gTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgctg                                   >  1:892736/1‑101 (MQ=255)
                                                            tttGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgctgc                                  >  1:3044813/1‑101 (MQ=255)
                                                            tttGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgctgc                                  >  2:238805/1‑101 (MQ=255)
                                                             ttGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgctgct                                 >  1:4580012/1‑101 (MQ=255)
                                                               gccgccTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctcc                               >  1:2534990/1‑101 (MQ=255)
                                                                 cgccTCGTAGCGGGCGGGGAGGGGGCCGGGGGCACC.CCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccct                             <  2:833655/101‑1 (MQ=255)
                                                                 cgccTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccct                             >  2:2633546/1‑101 (MQ=255)
                                                                 cgccTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccct                             <  1:620312/101‑1 (MQ=255)
                                                                 cgccTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCACCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctcccc                             >  2:661785/1‑100 (MQ=255)
                                                                   ccTCGTAGCGGGGGGTGAAGGGGCGGGGGGCACC.CCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccctcc                           <  1:983688/101‑1 (MQ=255)
                                                                   ccTCGTAACGGGCGGTGAAGGGGCCGGGGGCACC.CCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccctcc                           <  1:1146127/101‑1 (MQ=255)
                                                                   ccTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccctcc                           <  1:1509339/101‑1 (MQ=255)
                                                                    cTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCtccg                          >  1:1210347/1‑101 (MQ=255)
                                                                       gTACCGGGCGGTGAAGGGGCCGGGGGCACC.CCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCtccgtcc                       <  2:367258/101‑1 (MQ=255)
                                                                       gTAACGGGCGGTGAAGGGGCCGGGGGCACC.CCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCTGCCGCTGCCTAGCTGCTCCCtccgtcc                       <  2:1092532/101‑1 (MQ=255)
                                                                       gTAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCtccgtcc                       >  2:1341991/1‑101 (MQ=255)
                                                                        tAACGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGCCGGGCGCGCCCTTCCGGCCGCCGCCTAGCTGCCCCCtccgtccc                      >  2:1863396/1‑100 (MQ=255)
                                                                         aaCGGGCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGa                     >  2:1512051/1‑101 (MQ=255)
                                                                          aggggCGGTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAg                    <  2:1722037/99‑1 (MQ=255)
                                                                            gggCGGTGAAGGGGCCGGGGGCACC.CCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTg                  <  2:3213361/101‑1 (MQ=255)
                                                                             ggCGGTGAAGGGGCCGGGGGCACC.CCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTgg                 <  1:1165525/101‑1 (MQ=255)
                                                                                ggTGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCgg              >  1:681409/1‑101 (MQ=255)
                                                                                  tGAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTg            >  2:3416523/1‑101 (MQ=255)
                                                                                   gAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGc           >  1:4562302/1‑101 (MQ=255)
                                                                                   gAAGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGGGc           >  2:1319161/1‑101 (MQ=255)
                                                                                    aaGGGGCCGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCg          >  1:3657296/1‑101 (MQ=255)
                                                                                     aGGGGCCGGGGGCACC.CCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCgg         <  1:3473901/101‑1 (MQ=255)
                                                                                           cGGGGGCACC.CCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTg   <  2:3381246/101‑1 (MQ=255)
                                                                                           cGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTg   >  2:879192/1‑101 (MQ=255)
                                                                                           cGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTg   <  1:2181061/101‑1 (MQ=255)
                                                                                           cGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTg   >  2:3342723/1‑101 (MQ=255)
                                                                                           cGGGGGCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGGCTg   >  1:3586378/1‑101 (MQ=255)
                                                                                            gggggCA.ACCCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGGCTgc  >  1:554808/1‑101 (MQ=255)
                                                                                                     |                                                                                            
TATGCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGCA.A.CCCGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGC  >  NC_014034/1823641‑1823832

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