Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,823,743 +G 1.2% intergenic (+28/‑145) 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,7431.→G 1.2% 3.2 ‑3.2 ‑0.0337intergenic (+28/‑145)RCAP_rcc01687/RCAP_rcc01688HK97 family phage major capsid protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (184/149):  total (184/153)
Fisher's exact test strand distribution p-value = 4.16e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.97e-01

Only 100 of 419 total aligned reads displayed.

TGCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.CC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGCGCCCGAACGCC  >  NC_014034/1823643‑1823843
                                                                                                          |                                                                                                    
 gCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.cc.cc                                                                                                      <  1:8070343/101‑2 (MQ=255)
 gCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.cc.cc                                                                                                      >  2:5970410/1‑100 (MQ=255)
   gAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.cc.ccgg                                                                                                    <  2:6872439/101‑4 (MQ=255)
   gAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.cc.ccgg                                                                                                    <  1:3636696/101‑4 (MQ=255)
   gAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.cc.ccgg                                                                                                    >  2:2239442/1‑98 (MQ=255)
   gAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.cc.ccgg                                                                                                    >  2:2605540/1‑98 (MQ=255)
     gCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGNC.GGGGGCAA.A.CC.C.GGCt                                                                                                  >  2:3776827/1‑101 (MQ=255)
     gCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCCGGGGGGA.A.A.CC.C.GGcc                                                                                                  >  1:5010941/1‑100 (MQ=255)
     gCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.cc.ccggct                                                                                                  >  1:4127274/1‑96 (MQ=255)
     gCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.cc.ccggct                                                                                                  <  1:1633104/101‑6 (MQ=255)
      cAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTg                                                                                                 >  1:339584/1‑101 (MQ=255)
       aaGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTgc                                                                                                >  2:8176260/1‑101 (MQ=255)
           gggTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTgccgcc                                                                                            >  1:2247428/1‑101 (MQ=255)
             gTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCgg                                                                                          >  2:3462938/1‑101 (MQ=255)
             gTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCgg                                                                                          >  1:6855967/1‑101 (MQ=255)
             gTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGGGAAGGGGCC.GGGGGG.A.ACCC.C.GGCTGCCGCCgg                                                                                          >  1:2925390/1‑101 (MQ=255)
             gTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGGGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCgg                                                                                          >  2:4219619/1‑101 (MQ=255)
                ggcggcGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTc                                                                                       <  1:8201815/101‑1 (MQ=255)
                ggcggcGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTc                                                                                       <  2:6239720/101‑1 (MQ=255)
                ggcggcGATGTCAGCGATTTCGCCGCGATCAACCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTc                                                                                       <  1:6265243/101‑1 (MQ=255)
                 gcggcgATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCg                                                                                      <  2:8040652/101‑1 (MQ=255)
                   ggcgATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGcc                                                                                    >  2:5291575/1‑101 (MQ=255)
                    gcgATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCg                                                                                   >  1:2450330/1‑101 (MQ=255)
                     cgATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.CGCTGCGGCCGGTTCGCCgg                                                                                  >  2:4726198/1‑101 (MQ=255)
                      gATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGt                                                                                 >  1:3996171/1‑101 (MQ=255)
                       aTGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTc                                                                                >  2:4466236/1‑101 (MQ=255)
                       aTGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTc                                                                                <  2:7622479/101‑1 (MQ=255)
                       aTGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTc                                                                                >  1:1971717/1‑101 (MQ=255)
                       aTGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.CCAC.GGCTGCCGCCGGTTCGCCGGTc                                                                                >  2:509426/1‑101 (MQ=255)
                         gTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTcgg                                                                              <  2:4641889/101‑1 (MQ=255)
                               gATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCgggcgggc                                                                        >  1:2022146/1‑101 (MQ=255)
                               gATTTCGCCGCGATCAAGCTGTGGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.ACC.CC.C.GGCTGCCGCCGGTTCGCCGGTCgggcgggc                                                                        <  1:4404636/101‑1 (MQ=255)
                                aTTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGggcgggcg                                                                       >  1:9080238/1‑101 (MQ=255)
                                aTTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGggcgggcg                                                                       >  1:7794045/1‑101 (MQ=255)
                                 tttCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGgcgggcgc                                                                      >  1:7263713/1‑101 (MQ=255)
                                  ttCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGgcgggcgcg                                                                     >  2:517824/1‑101 (MQ=255)
                                     gccgcGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCCGCCGCCGGTCCGCCGGTCgggggggctgggcc                                                                  >  1:8488564/1‑95 (MQ=255)
                                        gcgATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCGGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTc                                                               >  2:6767275/1‑101 (MQ=255)
                                        gcgATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGNCGGTTCGCCGGTCGGGCGGGCGCGGGCTTc                                                               >  1:1370374/1‑101 (MQ=255)
                                         cgATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTcc                                                              >  1:4807668/1‑101 (MQ=255)
                                            tCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGc                                                           >  1:1018158/1‑101 (MQ=255)
                                             cAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGcc                                                          >  2:901980/1‑101 (MQ=255)
                                              aaGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGTCGGGCGCGGGCTTCCGGCCg                                                         >  1:5651151/1‑101 (MQ=255)
                                              aaGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCg                                                         >  1:61593/1‑101 (MQ=255)
                                              aaGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCg                                                         >  1:2696760/1‑101 (MQ=255)
                                               aGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGc                                                        >  1:8366315/1‑101 (MQ=255)
                                                gCTGTTGAGGTTTGCCGCCTCGTAGCGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCt                                                       <  2:8846938/101‑1 (MQ=255)
                                                gCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCt                                                       >  2:6577589/1‑101 (MQ=255)
                                                   gTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCgattcc                                                    >  2:2021418/1‑96 (MQ=255)
                                                   gTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGcc                                                    >  2:1829927/1‑101 (MQ=255)
                                                   gTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGcc                                                    >  1:4966485/1‑101 (MQ=255)
                                                    ttGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCt                                                   <  1:5189887/101‑1 (MQ=255)
                                                    ttGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCt                                                   <  2:2796063/101‑1 (MQ=255)
                                                     tGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.ACC.CC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTa                                                  <  1:1768938/101‑1 (MQ=255)
                                                        aGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.ACC.CC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgct                                               <  2:3103590/101‑1 (MQ=255)
                                                        aGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgct                                               <  1:598950/101‑1 (MQ=255)
                                                        aGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgct                                               <  1:6941555/101‑1 (MQ=255)
                                                        aGTTTGCCCCCTCGTAAGGGGCGGGGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgct                                               <  1:6373405/101‑1 (MQ=255)
                                                         gTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.ACC.CC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgctg                                              <  2:5543310/101‑1 (MQ=255)
                                                           ttGCGGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgctgct                                            <  2:1313870/101‑1 (MQ=255)
                                                            tGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctc                                           >  2:5317687/1‑101 (MQ=255)
                                                             gccgccTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.ACC.CC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctcc                                          <  2:791709/101‑1 (MQ=255)
                                                                gccTCGTACCGGGCGGTGAAGGGGCC.GGGGGC.AGC.CC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccctc                                       <  2:7162627/101‑1 (MQ=255)
                                                                 ccTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccctcc                                      <  2:3181524/101‑1 (MQ=255)
                                                                 ccTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccctcc                                      >  2:3649039/1‑101 (MQ=255)
                                                                    cGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCtccgtc                                   >  2:6728117/1‑101 (MQ=255)
                                                                     gTAACGGGCGGTTAAGGGGCC.GGGGGC.AGC.CC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccctagatcg                                  >  1:6529558/1‑95 (MQ=255)
                                                                     gTAACGGGCGGTGAAGGGGCC.GGGGGC.ACC.CC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCtccgtcc                                  <  2:3552494/101‑1 (MQ=255)
                                                                     gTAACGGGCGGTGAAGGGGCC.GGGGGC.ACC.CC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCtccgtcc                                  <  2:3304897/101‑1 (MQ=255)
                                                                     gTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCtccgtcc                                  >  1:3161834/1‑101 (MQ=255)
                                                                      tAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCtccgtccg                                 >  1:1216213/1‑101 (MQ=255)
                                                                      tAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCtccgtccg                                 >  1:7387074/1‑101 (MQ=255)
                                                                         cGGGCGGTGAAGGGGCC.GGGGGC.ACC.CC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGt                              <  1:4609620/101‑1 (MQ=255)
                                                                         cGGGCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAgg                              >  2:5797768/1‑100 (MQ=255)
                                                                          gggCGGTGAAGGGGCC.GGGGGC.ACC.CC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTg                             <  1:5747455/101‑1 (MQ=255)
                                                                           ggCGGTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTgg                            >  1:2138678/1‑101 (MQ=255)
                                                                              ggTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCgg                         >  2:2704423/1‑101 (MQ=255)
                                                                              ggTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCgg                         >  2:643628/1‑101 (MQ=255)
                                                                               gTGAAGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGt                        <  2:8493133/101‑1 (MQ=255)
                                                                                ggAAGGGGCC.GGGGGC.A.A.CC.CGGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTg                       <  2:4966325/100‑1 (MQ=255)
                                                                                 gAAGGGGCC.GGGGGC.ACC.CC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGc                      <  1:1438973/101‑1 (MQ=255)
                                                                                 gAAGGGGCC.GGGGGC.A.A.CC.CGGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGc                      <  2:6983581/101‑1 (MQ=255)
                                                                                  aaGGGGCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCg                     >  1:8267531/1‑101 (MQ=255)
                                                                                     gggCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCgggg                  >  2:1051076/1‑101 (MQ=255)
                                                                                       gCC.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCgggggc                >  2:2437949/1‑99 (MQ=255)
                                                                                        cc.GGGGGC.ACC.CC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCt               <  2:8618808/101‑1 (MQ=255)
                                                                                        cc.GGGGGC.A.ACCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCt               <  2:2553255/101‑1 (MQ=255)
                                                                                          ggggggC.A.C.CCAC.GGTAGCCGCCGTGTCGCCGGTCGGGCGGGCGCGGGGTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTg              <  2:9095271/100‑1 (MQ=255)
                                                                                           gggggC.A.CCCC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTgc             <  1:7991148/101‑1 (MQ=255)
                                                                                              ggggcacccc.c.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTgcg            <  1:5106631/93‑1 (MQ=255)
                                                                                                gccacccc.c.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGCGcc          <  2:74286/95‑1 (MQ=255)
                                                                                                  cacccc.c.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCTCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGCGCCCg        <  2:2187425/97‑1 (MQ=255)
                                                                                                  cacccc.c.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGCGCCCg        <  1:2833738/97‑1 (MQ=255)
                                                                                                   aaccc.c.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGCGCCCGa       >  1:5490627/4‑101 (MQ=255)
                                                                                                   aaccc.c.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGCGGCGGTGCGGGGCCTGCGCCCGa       >  1:644474/4‑101 (MQ=255)
                                                                                                     ccc.c.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGCGCCCGAAc     <  2:6283113/100‑1 (MQ=255)
                                                                                                     ccc.c.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGCGCCCGAAc     <  2:1452284/100‑1 (MQ=255)
                                                                                                       c.c.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGCGCCCGAACGc   <  1:86601/101‑1 (MQ=255)
                                                                                                       c.c.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGCGCCCGAACGc   <  1:8750709/101‑1 (MQ=255)
                                                                                                         c.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGCGCCCGAACGcc  <  2:1664829/101‑1 (MQ=255)
                                                                                                          |                                                                                                    
TGCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.CC.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGCGCCCGAACGCC  >  NC_014034/1823643‑1823843

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