Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0480A→G 2.1% 78.7 ‑5.2 ‑1.32423D203D (GAT→GACRCAP_rcc00185hypothetical protein
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (13/42):  ref base (1009/1354):  total (1026/1397)
Fisher's exact test strand distribution p-value = 5.30e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.59e-02

Only 100 of 3060 total aligned reads displayed.

GCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGAC  >  NC_014034/221948‑222148
                                                                                                    |                                                                                                    
 gttcgcgggctcgggcgggtccgcgggatcggtcggggccgggggaTCGGGCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                                                     <  1:5325521/91‑1 (MQ=11)
 gtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt                                                                                                     >  1:4350406‑M1/8‑101 (MQ=255)
 gtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAt                                                                                                     >  1:5936640‑M1/8‑101 (MQ=255)
  tcggtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                                    >  1:4211853‑M1/7‑101 (MQ=255)
   cggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcg                                                                                                   >  2:332993/1‑101 (MQ=11)
    ggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcgg                                                                                                  >  1:86374/1‑101 (MQ=11)
    ggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTGGGGtcgg                                                                                                  >  1:4080019/1‑101 (MQ=11)
     gtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAttggt                                                                                                 <  1:3270207‑M1/98‑1 (MQ=255)
      tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGatcggtg                                                                                                <  2:3596404‑M1/99‑2 (MQ=255)
      tcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                                                <  2:2216719/101‑1 (MQ=255)
      ggGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                                                >  1:6193994‑M1/3‑101 (MQ=255)
       cggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGAtcgggcg                                                                                               >  2:4029492/1‑101 (MQ=14)
        gggtcggtcgggtccgtgggatcggtcgggtccgtggggTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                              <  1:2202433/98‑1 (MQ=11)
          gtcggttgcgtgcgtgggatcggtcgggtccgtgggaTCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcggat                                                                                            <  2:708691/100‑3 (MQ=11)
           tcggtcgggtccgtgggatcggtcgggtccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCGGTGGGGTCGGGCGGGTGGGAGGGGtcggtccgatc                                                                                           >  2:5217935/1‑96 (MQ=11)
            cggtcgggtccgtggggtccgtcgggtccgtgggatCGGTCGGCGCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                                                          <  2:5896953/101‑1 (MQ=11)
            cggtcggggccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATcc                                                                                          <  1:138494/101‑1 (MQ=14)
            cggtcggatccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcggagcc                                                                                          >  2:5758059/1‑97 (MQ=11)
            cggtcggatccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAGCGGTGGGATCCGTGGGCccgggcggagcg                                                                                          >  2:5498895/1‑97 (MQ=11)
            cggtcgagtccgtgggagcggtcgggtccgtgggatCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                                                                                          <  2:6801435/101‑1 (MQ=11)
              gtcggatccgtgggatcggtcgggtccgtgggatcgGTCGGATCCGTGGGGTCGGTCGGATCCGTGGGCTCCGTGGGATCGGTTGGGtcgggcggcacggg                                                                                        >  2:4505129/1‑95 (MQ=11)
              gtcggatccgtgggatcggtcggatccgtaggatcgGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                                                                                        <  2:2655443/101‑1 (MQ=11)
                tgggtccgttgggtcggtcgggtccgtgggatctgtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtgg                                                                                      <  1:2666806/100‑3 (MQ=11)
                cgggtccgtgggatcggtcgggtccgtggggtccgtCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgg                                                                                      <  1:4506248‑M2/101‑3 (MQ=255)
                cgggtccgtgggatcggtcgggtccgtcggatcggtGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTCGGATCGGTCGGGGCCGtgg                                                                                      >  2:5610760/1‑99 (MQ=11)
                 gggtctgtgggatcggtcgggtccgtcggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGgg                                                                                     >  1:6869015/1‑101 (MQ=11)
                 gggtccgtgggatcggtcgggtccgtcggatcggtcGGGTCCGTGGGATCGGGCGGGTCCGTGGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCgcggg                                                                                     >  2:1850570/1‑97 (MQ=11)
                  ggtacgtgggatcggtcgggtccgtgggatcgggcgGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggga                                                                                    >  2:1114783‑M2/1‑97 (MQ=255)
                  gatccgtgagatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTAGGATCCGTGGGATCGGTCGGATCCGTGGGa                                                                                    <  1:1271283/99‑1 (MQ=11)
                   gtccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCGGTgggga                                                                                   <  2:2868349‑M2/101‑6 (MQ=255)
                    tttttgggatctgtcgggtccgtgggatcggtcgggtccGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatc                                                                                  >  2:6664832‑M2/5‑95 (MQ=255)
                    tccgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatc                                                                                  >  2:207029‑M2/1‑95 (MQ=255)
                    tccgtgggatcggtcgggtccggggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGGTCGGCCGGATCCGTGGGATCGGTCGGGGCCCTGGGAtc                                                                                  >  2:4636388/1‑101 (MQ=14)
                     ccgtcggatcggtcgcgtctgtggagtccgtcggATGCGAGGGATCGGGTGTATCCGTGGGACCGGCCAGATCCGTGGGATCGGTCGGGCCCGTGGGAtcg                                                                                 <  2:2753883/101‑1 (MQ=2)
                      cgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGgggcggg                                                                                >  1:1000254/1‑96 (MQ=11)
                      cgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggggtccg                                                                                >  1:6533844‑M2/1‑93 (MQ=255)
                        gtgggtcggttggatccgtgggatcgggcgtATCCGCGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGAtcggtc                                                                              <  1:517652/96‑1 (MQ=11)
                         gggatcggtcgggtccgtggggtccgtcggATCCGTGGGCTCGGTCGGATCCGTGGGCTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtccgccg                                                                             >  1:4768888/1‑96 (MQ=11)
                         gggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTgggatcgggcg                                                                             <  2:6942314‑M2/101‑12 (MQ=255)
                          ggatcggtcgggtccgtgggatcgggcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                                                                            <  1:3159439/101‑1 (MQ=17)
                           gatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTGGGGTCGGGGGGATCGgacggg                                                                           >  2:4933297/1‑101 (MQ=11)
                             tcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtc                                                                         <  2:6399337/101‑1 (MQ=255)
                              cggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCAGTGGGGTCGGTGGGGTCCGTAGGTTCGGTGGGGGCAGTGCGATCGgaggggtcc                                                                        >  1:2894289/1‑101 (MQ=2)
                              cggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcggatcc                                                                        >  1:6656237‑M2/1‑85 (MQ=255)
                                gtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgt                                                                      >  1:6768093‑M2/1‑83 (MQ=255)
                                 tcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtc                                                                     <  2:6260898/101‑1 (MQ=255)
                                  cgggtccgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgg                                                                    <  1:567999‑M2/101‑21 (MQ=255)
                                   gggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtggg                                                                   <  2:4163367‑M2/101‑22 (MQ=255)
                                   ggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtggg                                                                   <  1:2438714‑M2/101‑22 (MQ=255)
                                   ggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATGCGTGGGATCGGACGGGTTCGTGGGATCGgtcgggtccgtggg                                                                   >  2:892953/4‑101 (MQ=14)
                                      tccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggtc                                                                <  1:6374816/101‑1 (MQ=11)
                                      gccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggtc                                                                <  2:6135123/100‑1 (MQ=11)
                                         gtgggatccgtgtgATGCGTGGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTcgggtccgtcgggtccgt                                                             <  2:5068766/101‑1 (MQ=17)
                                           gggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTgg                                                           >  1:1106838/1‑101 (MQ=11)
                                           cgggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTgg                                                           >  2:1954236/11‑101 (MQ=11)
                                           cgggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTgg                                                           >  2:5979086/11‑101 (MQ=11)
                                            ggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgg                                                          <  1:3498262‑M2/101‑31 (MQ=255)
                                              atcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggt                                                        <  2:6545441‑M2/101‑33 (MQ=255)
                                               tcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtggggtcggtggggtc                                                       <  2:2538963‑M2/101‑34 (MQ=255)
                                                cggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCCGGGGGGTCCGTCGGGTCCGTGGGGtcg                                                      >  2:4294547/1‑101 (MQ=11)
                                                cggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtcc                                                      >  2:3233006/1‑100 (MQ=11)
                                                   gcggGTCCGTGGGATCGGTAGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtc                                                   <  1:1947874/100‑1 (MQ=11)
                                                       gtcggTGGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGAt                                               <  1:1579191/97‑1 (MQ=17)
                                                        tgggtcggATCGGTGGGGTACGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATc                                              <  1:4636627/95‑1 (MQ=17)
                                                          cGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtggggtccgtcgggtccgtgggatcgg                                            <  1:3751902‑M2/101‑45 (MQ=255)
                                                          cGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcgg                                            <  2:3921555‑M2/101‑45 (MQ=255)
                                                           gtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCgc                                           >  2:2408038/4‑100 (MQ=17)
                                                             gggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcg                                         <  1:5448302‑M2/101‑48 (MQ=255)
                                                             gggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcg                                         <  1:2793234‑M2/101‑48 (MQ=255)
                                                              ggATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgg                                        >  2:3815666‑M2/1‑53 (MQ=255)
                                                              ggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                                        <  1:4931227/101‑1 (MQ=255)
                                                               ggtcggtcggGTCCGTGGGATCGGTCGGAGCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                       <  1:2063368/99‑1 (MQ=17)
                                                               ggtccgtcggatCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgga                                       >  1:5558894‑M2/3‑52 (MQ=255)
                                                               gATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                       >  1:715245/1‑101 (MQ=17)
                                                                    gtggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                  >  1:3019287/4‑101 (MQ=255)
                                                                    gtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgt                                  <  1:2338262‑M2/101‑55 (MQ=255)
                                                                     tcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggagcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtg                                 <  1:3878267‑M2/101‑56 (MQ=255)
                                                                      cggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgggggatcggtcggggccgtgg                                <  2:1391371‑M2/101‑57 (MQ=255)
                                                                      cggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGCTCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                >  1:1738973/1‑101 (MQ=17)
                                                                       gggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                               <  2:4778439/101‑1 (MQ=255)
                                                                       ggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgggggatcggtcgggtccgtggg                               <  2:173793‑M2/101‑58 (MQ=255)
                                                                            cGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                          <  2:2980327/101‑1 (MQ=255)
                                                                             gTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgt                         <  1:5797770‑M2/101‑64 (MQ=255)
                                                                               aggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATGCGGCGGGTCGGtcg                       >  2:1004871/2‑99 (MQ=14)
                                                                                  gtccggcggATCCGGGGGACCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggat                    <  2:1551144/95‑3 (MQ=14)
                                                                                  aTCGGTCGGGTCCGTGGGATCGGTCGGGTCC.TAGGATCGGTCGGGTCCGTCGGGTCCCTGGGATCGGTCGGATCCGGGGGATCGGCCGGGTCCGTggggtg                   >  2:1249489/1‑99 (MQ=18)
                                                                                     ggtgggATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtcgc                 <  1:3955719‑M2/101‑72 (MQ=255)
                                                                                     gggcggggccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccg                 >  1:3596540‑M2/9‑30 (MQ=255)
                                                                                      gtcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtgga                >  2:774032/1‑100 (MQ=255)
                                                                                      gtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggg                >  2:3693105/1‑95 (MQ=14)
                                                                                      gtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggt                >  1:1846/1‑95 (MQ=14)
                                                                                       tcggactggggggATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtggtg               <  1:6063849/91‑1 (MQ=18)
                                                                                       tcggGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtg               >  2:3489797‑M2/1‑28 (MQ=255)
                                                                                         ggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtcgg             >  1:2559619/1‑92 (MQ=14)
                                                                                             ccGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGCCTGGTGGGAtcg         <  1:5058059/101‑1 (MQ=255)
                                                                                                tGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggtc      >  1:3401494/1‑101 (MQ=255)
                                                                                                 gggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcg     <  1:7020774‑M2/101‑84 (MQ=255)
                                                                                                 gggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatccgtcgggtccgggggatcggtcg     <  2:5359193‑M2/101‑84 (MQ=255)
                                                                                                 gggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggggccgtgggatcggtcg     <  2:4385417‑M2/101‑84 (MQ=255)
                                                                                                  ggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcgggcggatccgtgggatcggtggg    <  2:3096565‑M2/101‑85 (MQ=255)
                                                                                                    |                                                                                                    
GCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGAC  >  NC_014034/221948‑222148

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