Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0300A→G 7.6% 258.1 ‑6.0 ‑1.41672D209D (GAT→GACRCAP_rcc00185hypothetical protein
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (43/88):  ref base (716/821):  total (762/910)
Fisher's exact test strand distribution p-value = 2.49e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.73e-02

Only 100 of 2085 total aligned reads displayed.

TGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGC  >  NC_014034/221930‑222130
                                                                                                    |                                                                                                    
 gtccgtgggatcggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                                                     <  1:2363085‑M1/70‑1 (MQ=255)
 atccgtgggatcggtcgggtccgtgggatcggtcggaTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                                                     <  2:4357354‑M1/64‑1 (MQ=255)
    cgtgggatcggtcggatccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                                                  <  1:751097‑M1/67‑1 (MQ=255)
      tgggatcggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                                                >  1:3766092‑M1/27‑101 (MQ=255)
       gggatcggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                                               >  2:2195759‑M1/26‑101 (MQ=255)
           tccgtgggatcggtcgggtcggtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGGTCCGTGggatcggtcgggtc                                                                                           >  1:3618009‑M1/79‑98 (MQ=255)
           tGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                                                                           >  1:4119293/1‑101 (MQ=255)
            cggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                                                                          >  1:918019‑M1/21‑101 (MQ=255)
            ccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTcc                                                                                          <  1:1294808‑M1/24‑1 (MQ=255)
               tGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                                                                       >  2:2887691/1‑101 (MQ=255)
                cgggtccgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                                                      >  2:1997079‑M1/23‑101 (MQ=255)
                     cggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                                                                 <  1:3339707/101‑1 (MQ=255)
                     ccgttggatcggtggggtccctgggctcggtcgggtccgtaGGATCTGTCGGGTCCGTCGGATCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                                                                 <  2:3512555/96‑1 (MQ=11)
                       gtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGCGGATCGGTCGGGGCCGGGGGCTCGGCCGGGCCCGGGGGCtcggc                                                                               >  1:2081689/1‑100 (MQ=2)
                        tgtgatctttcgggtccgtgggatcggtcgggtccgtgggatcggTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                              <  2:2969352/92‑1 (MQ=14)
                        tggggtcggtcggaTCCCTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                              >  2:3887110‑M1/15‑101 (MQ=255)
                        tgggatgagtccgggccgtgggatccgtcggatccgtgggatcgGTCGGATCCGTGGGTTCGGTCGTGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                              <  2:1310893/93‑1 (MQ=11)
                        tgggagcgggcgggtccgtgggatcggtggggtccgtgGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                              <  1:1359842/99‑1 (MQ=11)
                         gggggccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAtcggtcg                                                                             <  2:247712/94‑1 (MQ=14)
                          ggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                            >  2:1220751‑M1/13‑101 (MQ=255)
                            gtcggcggggtcagggggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGt                                                                          <  1:2454539/94‑1 (MQ=11)
                            atcggtcgggtccgtgggatcggtcgggtccgtgggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGat                                                                          >  2:2579453‑M2/1‑99 (MQ=255)
                            atcggtcggatccgtgggatcggtcggatccgtaggATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGt                                                                          <  1:2010410/101‑1 (MQ=11)
                            atcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGt                                                                          <  1:4397927‑M1/91‑1 (MQ=255)
                             tcggtcgggtccgtgggatcggtcgggtccgtgggaTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGatc                                                                         >  1:3131679‑M2/1‑98 (MQ=255)
                             tcggtcgggtccgtgggatcggtcgggtccgtgggaTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGatc                                                                         >  1:4188267‑M2/1‑98 (MQ=255)
                              cggtcggatccgtgggatcggtcggatccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATcc                                                                        <  1:4739143/101‑1 (MQ=11)
                              cggccgggtccgtgggatcggtcgggtccgtcggatCGGTTGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTcc                                                                        <  1:2306226/101‑1 (MQ=11)
                               ggtgggatccgtgggatcggtcgggtccgtgggatcGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                                                       >  1:167116/1‑101 (MQ=11)
                               ggtcgggtccgtgggatcggtcgggtccgtgggatcGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccg                                                                       <  2:1561911‑M2/101‑6 (MQ=255)
                               ggtcgggtccgtgggatcggtcgggtccgtgggatcGGTCGGGACCGTGGGATCGGTCGGGTCAGTGGGGGCGGCCGGGTCAGCGGGTTAGGTCGGATCCg                                                                       >  1:4659337/1‑101 (MQ=11)
                               gggggtgtccgtgggatcggtcgggtccgttggatccgtcggGTCTGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                                                       <  1:4556089/95‑1 (MQ=11)
                                gtgggatcgggcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                                                      >  1:3015762‑M1/58‑101 (MQ=255)
                                gtcggatccgtgggatcggtcgggtccgtgggatcgGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgt                                                                      <  1:1132289‑M2/101‑7 (MQ=255)
                                 tcgggtccgtgggatcggtcgggtccgtggggtccgTCGGATCCGTGGGCTCGGCCGGCTCCGCGGGATCGGTCGGACCCGTCGGGTCGGTCGGGTCCgcg                                                                     >  1:2619500/1‑99 (MQ=14)
                                  gggttccctggcatcggtcgtgtccgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGtgg                                                                    <  2:2298554/97‑3 (MQ=11)
                                   gggtggctgggctcggtggtgtgcgctggatgggtcggGCGCGTGGGGTCGGGCGGGTGCCTCGGACCCGTGGGATCCGTGGGATCGGTCGGGTCCGTGgg                                                                   <  1:3692799/94‑1 (MQ=2)
                                      tccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatc                                                                >  2:1697116‑M2/1‑89 (MQ=255)
                                      tccgtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGatccgtgggatc                                                                >  2:2070499‑M2/1‑89 (MQ=255)
                                      tccgtgggatcggtcgggtccgtgggatcggacggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCTGatccatgggata                                                                <  1:82816‑M2/101‑13 (MQ=255)
                                       ccgtgggatcggtgggttccgggggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                                               <  1:1002194/101‑1 (MQ=11)
                                       ccgtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTGGGCTCCGTGGGATCGGTCGGATCCgcgggttcg                                                               >  1:3876439/1‑93 (MQ=11)
                                        ggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGCCCGTAGGAtcgg                                                              <  2:3605937‑M1/52‑1 (MQ=255)
                                        cgtgggatcggtcggatccgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcgg                                                              >  1:3959537/1‑101 (MQ=11)
                                         gtgggctcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcggtcgggtccgt                                                             >  2:1465770‑M2/1‑75 (MQ=255)
                                         gtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGGTCGGTCGGGTCCGTGGGAtcggg                                                             >  2:3669909/1‑100 (MQ=14)
                                         ggtgcgtcggtcgtgtccgtgggtcctgtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                             <  2:3347065/95‑1 (MQ=11)
                                          tgggatcggtcggatcggtgggatcggtgggGTCCGTGGGATCGGTGGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGAtcggtc                                                            <  2:3072929/101‑1 (MQ=11)
                                           gggatcggtcgggtccgtggggtcggtcggATCGGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                           <  2:1780134/101‑1 (MQ=17)
                                           gggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                           >  2:1937901/1‑101 (MQ=14)
                                           gggatcggtcggatccgtgggatcggtcggGTCAGAGGGATCGGTCGGGTCCGTGGGATAGGACGGGACCGTGGGAGCGGGCGGAGCCGTGGGAtcgagcg                                                           >  2:344489/1‑97 (MQ=2)
                                            ggatcggtcgggtccgggggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtcgg                                                          <  1:222187/101‑1 (MQ=11)
                                            ggatccgtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                          <  2:722141/101‑1 (MQ=11)
                                              gtcggtggggtccgtaggttcggtcggGTCCGTGGGGTCGGTTTGGTCCGTAGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggt                                                        <  1:2141614/100‑1 (MQ=11)
                                              atggggcgggtccgggggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggt                                                        <  1:671760/95‑1 (MQ=11)
                                                 ggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGACCGTGGGATCGGTCGGGTCGGTTTGATCGGTCGGGTCCGTGGGATCGgtcgggtccg                                                     >  1:4056312/1‑101 (MQ=11)
                                                 ggtcggacccggggggtcgggaggTTCCGTCGGATCCGTCGGGCCGGTGGGATGGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatccg                                                     <  2:2662078/101‑1 (MQ=2)
                                                  gtcggggccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgt                                                    >  1:284640‑M2/1‑77 (MQ=255)
                                                   tcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtg                                                   <  1:4267873‑M2/101‑26 (MQ=255)
                                                    cgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtgg                                                  >  1:1741346/1‑99 (MQ=11)
                                                      ggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAGCGGTCGGGTCCGTGGGATCGggcggctccgtggga                                                >  2:1390906/1‑100 (MQ=11)
                                                      ggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtggga                                                >  2:4157262‑M2/1‑73 (MQ=255)
                                                         ccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcggtcgggtccgtgggatcggccgggtcc                                             <  2:632378‑M2/101‑43 (MQ=255)
                                                          cgtggggtccgtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTCGGGTCGgtggggtccgtggggtcgg                                            >  2:3761990/1‑99 (MQ=11)
                                                           gtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGtcggggccgtggggtgcgg                                           >  2:917942/1‑100 (MQ=11)
                                                           gtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGGCCGTGGGATCGGTcgggtccgtggggtccgt                                           >  1:1192172/1‑101 (MQ=14)
                                                           gtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggg                                           <  2:2124443‑M2/101‑34 (MQ=255)
                                                            tgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtg                                          <  2:3255969‑M2/101‑35 (MQ=255)
                                                            tgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATGGGTCGGGTCCGTGGGATCCGTGGGGTCCTTGGGGTCGGTg                                          >  2:4529966/1‑101 (MQ=11)
                                                             gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcg                                         >  2:1881419‑M2/1‑66 (MQ=255)
                                                                  ctgttgcaTCGGTCGGATCGTTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcg                                    <  2:4304411/94‑1 (MQ=14)
                                                                  cggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgagatc                                    <  1:2543599‑M2/101‑41 (MQ=255)
                                                                   ggtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcggatccg                                   >  1:3232752‑M2/1‑60 (MQ=255)
                                                                      cggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcg                                <  2:553441/101‑1 (MQ=14)
                                                                      cggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgg                                >  2:4438724‑M2/1‑57 (MQ=255)
                                                                      cggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgtgg                                >  1:1449835/1‑99 (MQ=11)
                                                                       gggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCCGTGGGGTCCGTCGGGTCCGTGGGGtcggtggg                               >  2:4073682/4‑101 (MQ=11)
                                                                       ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcggatccgtggg                               >  2:3990420‑M2/1‑56 (MQ=255)
                                                                       ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgtggg                               >  1:2719566/1‑101 (MQ=11)
                                                                         gtcCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggat                             >  1:573563‑M2/2‑54 (MQ=255)
                                                                         gtCCGTGGGATCGGGCGAATCACTGGCATCGGTTGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAt                             <  2:4650475/100‑1 (MQ=11)
                                                                         atcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGAt                             <  1:3435277‑M1/85‑1 (MQ=255)
                                                                          tCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatc                            >  2:1138464‑M2/1‑53 (MQ=255)
                                                                          tCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtaggatcggtcgggtc                            >  1:1859743‑M2/1‑42 (MQ=255)
                                                                           ccGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtaggatcggtcgggtcc                           <  1:1553634‑M2/101‑61 (MQ=255)
                                                                             gTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggt                         <  1:4489754‑M2/101‑52 (MQ=255)
                                                                              ttggAGCGGTTGGGTGCGTGGGTTCGGTCGGCTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTg                        <  1:3501039/99‑1 (MQ=14)
                                                                              taggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGGTCGGGCGGGTCCGtc                        >  2:4752290/3‑100 (MQ=14)
                                                                                   tccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtc                   >  1:1053660‑M2/1‑44 (MQ=255)
                                                                                   gcggtcggGTCCGTCGGATCGGTCGGATCGGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtc                   <  1:1946236/100‑1 (MQ=14)
                                                                                    cggtcgtGTCCGTCGGATGGGTCGGATCCGGGGGAGCGGTCGGGTCCGGCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGGCGGGTCCGTCGGAtcg                  <  2:3172964/101‑1 (MQ=18)
                                                                                    ccgtcggatCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtcc                  <  2:796147‑M2/99‑59 (MQ=255)
                                                                                     ggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGGTCCGTGGGGTCCGTCGGGTCCGGGGGAtggg                 >  2:10076/1‑101 (MQ=11)
                                                                                       tcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgta               >  2:2752792‑M2/1‑40 (MQ=255)
                                                                                            tccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatc          >  2:2411225‑M2/1‑35 (MQ=255)
                                                                                             ccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcg         <  2:2339964‑M2/101‑68 (MQ=255)
                                                                                                 cggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTAGGTCGGGTCCGTgg     >  2:106404/2‑101 (MQ=17)
                                                                                                  ggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg    <  2:4200329/101‑1 (MQ=255)
                                                                                                  ggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTggg    >  1:2927499/1‑101 (MQ=18)
                                                                                                    aTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggat  <  1:4330871/101‑3 (MQ=14)
                                                                                                    |                                                                                                    
TGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGC  >  NC_014034/221930‑222130

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