Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,012 A→G 4.6% D215D (GAT→GAC RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0120A→G 4.6% 277.0 ‑1.0 0.03696D215D (GAT→GACRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (80/100):  ref base (1690/1819):  total (1775/1921)
Fisher's exact test strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 4236 total aligned reads displayed.

GATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAT  >  NC_014034/221912‑222112
                                                                                                    |                                                                                                    
       gggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtcg                                                                                               >  2:7528174‑M1/38‑101 (MQ=255)
       gggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                                                                               <  1:3928875‑M1/64‑1 (MQ=255)
        ggatcggtgggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                              <  1:3765508‑M1/65‑1 (MQ=255)
            cggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                                                          <  1:1312118‑M1/69‑1 (MQ=255)
            cggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                                                                                          >  1:3000604‑M1/33‑101 (MQ=255)
            ccgtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                                                          >  2:6862854‑M1/33‑101 (MQ=255)
            ccgtcggatccgggggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                                                          >  1:5503636‑M1/33‑101 (MQ=255)
             ggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGGTCCg                                                                                         <  2:6228753‑M1/70‑1 (MQ=255)
             ggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGCTCGGTCGGATCCGTGGGATCGGTCGGGTCCg                                                                                         <  1:8416734‑M1/70‑1 (MQ=255)
                ggTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGGCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGtcg                                                                                      >  2:203404/1‑99 (MQ=255)
                   ttccgtgggctcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                                   >  2:4070565‑M1/26‑101 (MQ=255)
                   gtccgtgggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                                   >  2:42399‑M1/26‑101 (MQ=255)
                     ccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                                                                 >  1:2120434‑M1/24‑101 (MQ=255)
                     ccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                                                                 >  1:3452243‑M1/24‑101 (MQ=255)
                       gtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                                                               <  2:6065041‑M1/80‑1 (MQ=255)
                       gtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                                                               <  1:1137202‑M1/80‑1 (MQ=255)
                       gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggc                                                                               >  2:100805/1‑100 (MQ=255)
                           gatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGa                                                                           <  1:7117675‑M1/84‑1 (MQ=255)
                             gcggtcgggtccgtgGGATCGGTCCGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGggatcggtcgggtc                                                                         >  1:5967831‑M1/16‑98 (MQ=255)
                               ggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCg                                                                       >  1:7650089‑M1/14‑101 (MQ=255)
                                gtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGt                                                                      <  2:1348436‑M1/89‑1 (MQ=255)
                                  cgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTgg                                                                    <  2:3420427‑M1/91‑1 (MQ=255)
                                   gggtcgggcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTcgg                                                                   <  1:1000386/98‑1 (MQ=11)
                                     gtcggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt                                                                 <  2:5844809/100‑1 (MQ=11)
                                     gtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt                                                                 <  1:8601137‑M1/94‑1 (MQ=255)
                                     gtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                 <  2:3230299‑M1/94‑1 (MQ=255)
                                      tccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                                                                >  1:3140366‑M1/7‑101 (MQ=255)
                                      tccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                                                                >  1:4826645‑M1/7‑101 (MQ=255)
                                      tccgtgGGATCGGGCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                >  2:8345471‑M1/7‑101 (MQ=255)
                                       cggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtcg                                                               >  1:4420920/1‑101 (MQ=11)
                                       cggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcg                                                               <  1:2280142/101‑1 (MQ=11)
                                       ccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATccgtggggtcc                                                               <  2:3586874‑M1/96‑5 (MQ=255)
                                       ccgggggatcggtcggggccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcg                                                               <  1:8925823/96‑1 (MQ=11)
                                        ggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                                              >  2:5345530/1‑101 (MQ=255)
                                        cgtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                              >  1:3625964/2‑101 (MQ=11)
                                         gtcggTTCGGTGGCATCGGTGGGACCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                             <  2:7013645/101‑1 (MQ=11)
                                         gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                                             >  1:5616286/1‑101 (MQ=255)
                                         gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGAtcggg                                                             >  2:8678673/1‑100 (MQ=14)
                                          gtgcatcggtcgggtccgtggggcccgtcgggtccgtgggATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                            <  2:4806484/97‑1 (MQ=11)
                                            ggatcggtcgggtccgtggggtccgtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTGGGAtcggtcgg                                                          >  2:8517297/1‑101 (MQ=11)
                                              atcggtgggctccgtgggatcggtcgggtccgtgggATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGt                                                        <  1:6374298/101‑1 (MQ=11)
                                              atcggtcgggtccgtgggatcggtcgggtccgtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGGCGGGTCCGTCGGATCGGTCgggc                                                        >  2:3069314/1‑100 (MQ=11)
                                              atcggtcgggtccgtgggatcggtcgggtccgtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGt                                                        <  2:5540009/101‑1 (MQ=11)
                                               tcggtgggatccgtggggtcggtcgggtccgtgggaTCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc                                                       <  1:8747396/101‑1 (MQ=11)
                                               tcggtcgggtcggtgggatcggtcgggtccgtgggaTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatc                                                       >  2:8828219‑M2/1‑87 (MQ=255)
                                               gtcgtggcgtccttgggattggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                                       <  2:8194910/93‑1 (MQ=11)
                                                cgttcgggtccgtcggatcggtcgggtccgtgggatcggTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATcc                                                      <  2:2311708/98‑1 (MQ=11)
                                                cggtcgggtccgtgggatcgggcgggtccgtgggatCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                                                      <  1:8276046/101‑1 (MQ=11)
                                                cggcgcgttccgtcggatccgtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                      <  2:5207546/93‑1 (MQ=11)
                                                 ggtcgggtccgtgggatcggtcgggtccgtggggtcCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTAGGGTcag                                                     >  2:8058766/1‑99 (MQ=14)
                                                  gtcgggtccgtgggatcggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCgg                                                    >  2:564286/1‑100 (MQ=11)
                                                  gtcgggtccgtcggatccgtcgggtccgtcggatctGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                                    <  2:2172443/101‑1 (MQ=11)
                                                  gtcggatccgtggggtccgtcgggtccgtgggatcgGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                                    >  2:7709590/1‑101 (MQ=11)
                                                  gtcggatccgtgggatcggtcggatccgtgggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGt                                                    <  1:2793140/101‑1 (MQ=11)
                                                   tcgggtccgtgggatcggtcgggtccgtgggatcggTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTg                                                   >  2:6691069‑M2/1‑100 (MQ=255)
                                                    cgggtccgtgggatcggtcgggtccgtgggatcggtCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgg                                                  <  2:3909478‑M2/101‑3 (MQ=255)
                                                     gggtccgtgggatcggtcgggtccgtggggtccgtcGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg                                                 <  1:4094674‑M2/101‑4 (MQ=255)
                                                     gggtccgtgggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTGGGGTCGGTCGGGTCAGTGgg                                                 >  2:2821654/1‑101 (MQ=11)
                                                     ggctccgggggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg                                                 >  1:4715170‑M2/1‑98 (MQ=255)
                                                      ggtccgtgggatcggtcgggtccgtggggtccgtcgGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                                >  2:4712492/1‑101 (MQ=11)
                                                      ggtccgtgggatcggtcgggtccgtggggtccgtcgGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggga                                                <  1:7464906‑M2/101‑5 (MQ=255)
                                                        ttcgtgggatcggtcgggtccgtggggtcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                                              <  2:7151378/99‑1 (MQ=14)
                                                        tccgtggggtcggtcgggtccggcgggtccgtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                              <  1:5207596/101‑1 (MQ=14)
                                                        tccgtggggtcggtcggggccgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                              <  1:3261904/101‑1 (MQ=11)
                                                         ccgtggggtccgtcgggtccgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAtcg                                             >  1:1459407/1‑101 (MQ=11)
                                                           gtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggt                                           <  1:6963822‑M2/101‑10 (MQ=255)
                                                           gtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcggtcgggtccgt                                           <  2:1181185‑M2/101‑27 (MQ=255)
                                                            ggggatgggtcgggtgcgtgggatcggtcggGGCCGTGGGAGGGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTTGGATCGGTCGGGTCCGTGGGAtcggtc                                          <  1:4259949/100‑1 (MQ=11)
                                                                  cggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatcc                                    >  1:261796‑M2/1‑85 (MQ=255)
                                                                  cggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcggatcc                                    <  2:2475111‑M2/101‑17 (MQ=255)
                                                                  cggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtcc                                    <  1:6182360/101‑1 (MQ=11)
                                                                    gtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgg                                  >  2:4973522/1‑100 (MQ=11)
                                                                      cgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgg                                <  1:175062‑M2/101‑21 (MQ=255)
                                                                       gggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtggg                               >  2:1682376/1‑101 (MQ=11)
                                                                       ggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgg                               >  1:6816354/4‑101 (MQ=255)
                                                                        ggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatccgtgggg                              >  2:7314862/1‑101 (MQ=11)
                                                                         gtccgtgggatcggtcggGTCCGTGGGATCGGTCGGCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggat                             >  2:1535300‑M2/1‑78 (MQ=255)
                                                                         gtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggat                             <  2:1847765‑M2/101‑24 (MQ=255)
                                                                               gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcg                       <  1:5348658‑M2/101‑30 (MQ=255)
                                                                               gggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTgg                       >  2:2637222/1‑101 (MQ=11)
                                                                                 gatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTgggg                     >  2:7862810/1‑100 (MQ=11)
                                                                                   tcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtc                   <  2:4403120/101‑1 (MQ=14)
                                                                                   tccgtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtc                   <  1:1592576/98‑1 (MQ=11)
                                                                                    gggtcgggTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtcc                  >  1:2455537‑M2/2‑67 (MQ=255)
                                                                                    cggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcggtccc                  <  1:7801453‑M2/101‑35 (MQ=255)
                                                                                    cggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcg                  <  1:6674795/101‑1 (MQ=14)
                                                                                    cggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtcc                  <  1:9107317/101‑2 (MQ=11)
                                                                                      gtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgt                <  1:4571575‑M2/101‑37 (MQ=255)
                                                                                      gtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgg                >  2:4531828/1‑98 (MQ=11)
                                                                                       tcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtggggtcggtcgggtccggg               <  2:6550259‑M2/101‑38 (MQ=255)
                                                                                       tcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtg               <  1:5000737‑M2/101‑38 (MQ=255)
                                                                                        cggATCCGTGGGCTCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgtgg              >  1:8151586/1‑99 (MQ=11)
                                                                                         gggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcgg             >  1:8004283/4‑101 (MQ=14)
                                                                                         ggatgggTGGGCTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggtcgg             <  2:9037695/95‑1 (MQ=18)
                                                                                         ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccgtggg             >  2:6411984/1‑101 (MQ=11)
                                                                                           aTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAt           >  2:3735402/1‑101 (MQ=11)
                                                                                            tCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATc          <  1:7267523/101‑1 (MQ=11)
                                                                                             cggtgGGGTCGGTCGGGTCCGTGGGATCGGTGGGATCCGTGgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtaggatcggtcgggtcc         >  1:1941522‑M2/3‑41 (MQ=255)
                                                                                             ccGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtaggatcggtcgggtcc         >  1:580904‑M2/1‑41 (MQ=255)
                                                                                                  ggATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtaggatcggtcgggtccgtcgg    <  1:6503378‑M2/101‑66 (MQ=255)
                                                                                                    |                                                                                                    
GATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAT  >  NC_014034/221912‑222112

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