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% 87.5 191intergenic (+25/‑148)RCAP_rcc01687/RCAP_rcc01688HK97 family phage major capsid protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (108/55):  ref base (0/28):  total (108/83)
Fisher's exact test strand distribution p-value = 3.37e-12
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.98e-01

Only 100 of 224 total aligned reads displayed.

CTATGCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.C.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGCGCCCG  >  NC_014034/1823640‑1823837
                                                                                                        |                                                                                                    
ctATGCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.a.a                                                                                                       >  2:1384787/1‑101 (MQ=255)
 tATGCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.c                                                                                                     >  2:4788886/1‑101 (MQ=255)
    gCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.c.c.cc                                                                                                <  2:2388593/101‑2 (MQ=255)
    gCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.c.c.cc                                                                                                <  2:4630518/101‑2 (MQ=255)
      gAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.c.c.ccgg                                                                                              <  1:1767085/101‑4 (MQ=255)
      gAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.c.c.ccgg                                                                                              <  1:34206/101‑4 (MQ=255)
      gAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.c.c.ccgg                                                                                              >  1:4208149/1‑98 (MQ=255)
       aGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.c.c.ccggc                                                                                             <  2:2265635/101‑5 (MQ=255)
          aaGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTgc                                                                                          <  1:1683964/101‑1 (MQ=255)
           aGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTGGTAACGGGCGGTGAAGGGGCC.GGGGGG.A.ACC.C.C.GGCTgcc                                                                                         >  1:1900821/1‑101 (MQ=255)
              gggTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTgccgcc                                                                                      >  1:1201357/1‑101 (MQ=255)
              gggTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTgccgcc                                                                                      >  1:4180810/1‑101 (MQ=255)
                 tGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTTTTGAAGTTTGCCGCCTCGTAACGGGCGGGGAAGGGGCC.GGGGGC.A.C.C.CGC.GGCTGCCGCCGGt                                                                                   <  2:3632919/101‑1 (MQ=255)
                 tGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGt                                                                                   >  1:3896654/1‑101 (MQ=255)
                  gggCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGtt                                                                                  <  1:2258892/101‑1 (MQ=255)
                   ggcggcGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTc                                                                                 <  1:2622791/101‑1 (MQ=255)
                    gcggcgATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCg                                                                                <  1:3060131/101‑1 (MQ=255)
                     cggcgATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGc                                                                               <  1:5748727/101‑1 (MQ=255)
                       gcgATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCg                                                                             <  1:2054373/101‑1 (MQ=255)
                        cgATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGTGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCgg                                                                            <  2:821188/101‑1 (MQ=255)
                        cgATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCgg                                                                            <  1:2492291/101‑1 (MQ=255)
                        cgATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCgg                                                                            >  2:1733185/1‑101 (MQ=255)
                         gATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCggg                                                                           >  2:2866406/1‑100 (MQ=255)
                         gATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGt                                                                           <  1:5217829/101‑1 (MQ=255)
                          aTGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTc                                                                          >  1:657972/1‑101 (MQ=255)
                          aTGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTc                                                                          <  1:5153479/101‑1 (MQ=255)
                           tGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTcg                                                                         <  2:3189559/101‑1 (MQ=255)
                           tGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTcg                                                                         >  1:4261229/1‑101 (MQ=255)
                            gTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.CAC.C.GGCTGCCGCCGGTTCGCCGGCcgg                                                                        >  2:5162487/1‑101 (MQ=255)
                             tCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGCAA.A.C.C.C.GGCTGCCGCCGGTTCGCCGGTcggg                                                                       >  1:2724098/1‑101 (MQ=255)
                              cAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTcgggc                                                                      >  1:4852547/1‑101 (MQ=255)
                                   aTTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGggcgggcg                                                                 <  1:1527022/101‑1 (MQ=255)
                                   aTTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGggcgggcg                                                                 >  2:1842626/1‑101 (MQ=255)
                                     ttCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGgcgggcgcg                                                               >  1:4021594/1‑101 (MQ=255)
                                       cgccgcGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGgcgggcgcggg                                                             >  1:111674/1‑101 (MQ=255)
                                        gccgcGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGgcgggcgcgggc                                                            >  1:53462/1‑101 (MQ=255)
                                        gccgcGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGGCGGgcgggcgcgggc                                                            >  1:4914595/1‑101 (MQ=255)
                                          cgcgATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCtt                                                          >  1:5583936/1‑101 (MQ=255)
                                          cgcgATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCtt                                                          <  1:1014338/101‑1 (MQ=255)
                                           gcgATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.ACC.C.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTc                                                         <  1:89656/101‑1 (MQ=255)
                                           gcgATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGGGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGACGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTc                                                         >  2:5252640/1‑101 (MQ=255)
                                              aTCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCgg                                                      <  1:366544/101‑1 (MQ=255)
                                               tCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCCGc                                                     >  1:2926097/1‑101 (MQ=255)
                                                cAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGcc                                                    >  1:1309448/1‑101 (MQ=255)
                                                cAAGCTGTGGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.ACC.C.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGcc                                                    <  2:2130929/101‑1 (MQ=255)
                                                 aaGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCg                                                   >  1:5531306/1‑101 (MQ=255)
                                                  aGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGc                                                  >  2:3689037/1‑101 (MQ=255)
                                                   gCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.ACC.C.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCt                                                 <  2:1084170/101‑1 (MQ=255)
                                                    cTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGGCGGTCGGGCGGGGGCGGGCTTCCGGCCGCTg                                                >  2:4221594/1‑101 (MQ=255)
                                                    cTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTg                                                >  1:3532695/1‑101 (MQ=255)
                                                    cTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTg                                                <  2:3811765/101‑1 (MQ=255)
                                                    cTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTg                                                <  1:2181404/101‑1 (MQ=255)
                                                    cTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTg                                                <  2:5274824/101‑1 (MQ=255)
                                                     tGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGc                                               >  1:5775223/1‑101 (MQ=255)
                                                      gTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGcc                                              <  1:5274427/101‑1 (MQ=255)
                                                      gTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGcc                                              <  1:1380626/101‑1 (MQ=255)
                                                         gAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAg                                           >  2:1658147/1‑101 (MQ=255)
                                                             tttGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgctgc                                       >  2:4509190/1‑101 (MQ=255)
                                                             tttGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgctgc                                       <  1:445155/101‑1 (MQ=255)
                                                             tttGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAgctgc                                       >  1:3323550/1‑101 (MQ=255)
                                                               tGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.ACC.C.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctc                                     <  2:3509746/101‑1 (MQ=255)
                                                               tGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctc                                     >  2:1409564/1‑101 (MQ=255)
                                                               tGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctc                                     >  1:3990587/1‑101 (MQ=255)
                                                                gccgccTCGTAACGGGCGGTGAGGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctcc                                    <  2:2671344/101‑1 (MQ=255)
                                                                gccgccTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctcc                                    >  2:1233035/1‑101 (MQ=255)
                                                                 ccgccTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.ACC.C.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccc                                   <  2:2166029/101‑1 (MQ=255)
                                                                 ccgccTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccc                                   >  1:2106718/1‑101 (MQ=255)
                                                                  cgccTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccct                                  <  2:4357018/101‑1 (MQ=255)
                                                                  cgccTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGGCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctcccc                                  >  2:4984782/1‑100 (MQ=255)
                                                                   gccTCGTACCGGGGCGTGAAGGGGCC.GGGGGC.ACC.C.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccctc                                 <  1:4995857/101‑1 (MQ=255)
                                                                   gccTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccctc                                 <  2:3523590/101‑1 (MQ=255)
                                                                   gccTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGctccctc                                 >  1:3485473/1‑101 (MQ=255)
                                                                      tCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCtccgt                              >  2:557924/1‑101 (MQ=255)
                                                                       cGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCtccgtc                             >  1:2402133/1‑101 (MQ=255)
                                                                       cGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCtccgtc                             <  1:1549020/101‑1 (MQ=255)
                                                                        gTAACGGGCGGTGAAGGGGCC.GGGGGC.ACC.C.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCtccgtcc                            <  1:4379781/101‑1 (MQ=255)
                                                                        gTAACGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTGCGGCCGCTGCCTTGCTGCTCCCtccggcc                            >  2:2207804/1‑101 (MQ=255)
                                                                         taggcggCGGGGAACGGGCC.GGGGGA.A.A.C.C.CGGGCTACCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCtccgtccg                           <  1:5267956/96‑1 (MQ=255)
                                                                          aaCGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCtccgtcccc                          >  2:4884370/1‑99 (MQ=255)
                                                                          aaCGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGa                          >  2:3057834/1‑101 (MQ=255)
                                                                           aCGGGCGGTGAAGGGGCC.GGGGGC.ACC.C.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAg                         <  2:5774928/101‑1 (MQ=255)
                                                                           aCGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAg                         <  2:111654/101‑1 (MQ=255)
                                                                           aCGGGCGGGAAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAg                         <  1:2866567/101‑1 (MQ=255)
                                                                            cGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGt                        >  1:2005099/1‑101 (MQ=255)
                                                                            cGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGt                        >  1:4531412/1‑101 (MQ=255)
                                                                            cGGGCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGt                        >  1:854624/1‑101 (MQ=255)
                                                                             gggGGGGGAGGGGGCG.GGGGGC.ACC.C.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTg                       <  1:3922818/101‑1 (MQ=255)
                                                                             gggCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTg                       <  1:2463093/101‑1 (MQ=255)
                                                                              ggCGGTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTgg                      >  2:4204345/1‑101 (MQ=255)
                                                                                 ggTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCgg                   >  1:4632255/1‑101 (MQ=255)
                                                                                  gTGAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTACCTGCTCCCTCCGTCCGAGTGGCggg                  >  2:3406465/1‑100 (MQ=255)
                                                                                   gggaGGGGCCGGGGGGC.A.C.C.C.CGGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGt                  <  2:3990343/98‑1 (MQ=255)
                                                                                    gAAGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGc                >  2:669496/1‑101 (MQ=255)
                                                                                     aaGGGGCC.GGGGGC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCg               >  1:1449909/1‑101 (MQ=255)
                                                                                           cc.GGGGGC.ACC.C.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCt         <  2:3018976/101‑1 (MQ=255)
                                                                                            cggggggC.A.C.C.C.CGGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCt         <  2:4813750/99‑1 (MQ=255)
                                                                                            c.GGGGGC.ACC.C.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTg        <  1:4793833/101‑1 (MQ=255)
                                                                                             ggggggC.A.CCC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTg        <  1:3112422/100‑1 (MQ=255)
                                                                                              gggggC.A.ACC.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTgc       >  2:2295912/1‑101 (MQ=255)
                                                                                                   c.A.A.C.C.CGGGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGCGCCCg  <  2:3261154/101‑1 (MQ=255)
                                                                                                        |                                                                                                    
CTATGCGAGCAAGCGGGTGGGCGGCGATGTCAGCGATTTCGCCGCGATCAAGCTGTTGAAGTTTGCCGCCTCGTAACGGGCGGTGAAGGGGCC.GGGGGC.A.A.C.C.C.GGCTGCCGCCGGTTCGCCGGTCGGGCGGGCGCGGGCTTCCGGCCGCTGCCTAGCTGCTCCCTCCGTCCGAGTGGCGGTGCGGGGTCTGCGCCCG  >  NC_014034/1823640‑1823837

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