Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,979 G→C 8.6% T226T (ACC→ACG RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9790G→C 8.6% 103.7 ‑0.6 ‑0.2615T226T (ACC→ACGRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (31/23):  ref base (294/267):  total (325/290)
Fisher's exact test strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.93e-01

Only 100 of 1648 total aligned reads displayed.

CCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTC  >  NC_014034/221879‑222081
                                                                                                    |                                                                                                      
 cgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcgg                                                                                                       <  1:3289917‑M1/25‑1 (MQ=255)
 cgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcggtcgggtccgtgggatcggtcgggtccgtgggatcgg                                                                                                       <  2:4076970‑M1/25‑1 (MQ=255)
 cGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcgg                                                                                                       <  1:1551524/101‑1 (MQ=255)
  gtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggt                                                                                                      <  1:3493373‑M1/26‑1 (MQ=255)
   aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                                                                     <  1:467706/101‑1 (MQ=255)
    ggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcg                                                                                                    <  1:1762360‑M1/28‑1 (MQ=255)
    cggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcg                                                                                                    <  1:142061‑M1/28‑1 (MQ=255)
       atcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggat                                                                                                 <  2:480605‑M1/31‑3 (MQ=255)
           gtgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgt                                                                                             >  2:2286988/1‑101 (MQ=255)
            tcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggagcggtcgggtccgtg                                                                                            <  2:4008105‑M1/36‑1 (MQ=255)
               ggtccgtcggatcggtgggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtggga                                                                                         >  1:4280082‑M1/63‑101 (MQ=255)
                gtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggat                                                                                        >  2:1009617‑M1/62‑101 (MQ=255)
                  tggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                                                                      <  1:820189/101‑1 (MQ=255)
                    gtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgt                                                                                    <  1:1688325‑M1/44‑7 (MQ=255)
                     tCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCgggtccgtgggatcggtcgggtccgtgggatcggtc                                                                                   <  1:3347697/101‑1 (MQ=255)
                      gggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcg                                                                                  >  2:2933530‑M1/56‑101 (MQ=255)
                         cGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                                                                               >  1:3577085/1‑101 (MQ=255)
                          gTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGGCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATAGGTCGGATc                                                                              >  1:2043107/1‑101 (MQ=255)
                           tggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                                             >  2:2215950‑M1/51‑101 (MQ=255)
                              ccGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                                                          >  1:4486190/1‑101 (MQ=255)
                                  gtcggtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAt                                                                      >  1:2748263‑M1/44‑101 (MQ=255)
                                  gtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGgaa                                                                      <  2:1102170‑M1/58‑2 (MQ=255)
                                    ccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                                                    <  2:2842967‑M1/60‑1 (MQ=255)
                                     gcgcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                                                   >  2:2639057/1‑101 (MQ=255)
                                      gtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGatcggc                                                                  <  2:2768152‑M1/62‑2 (MQ=255)
                                        gggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                <  1:40914‑M1/64‑1 (MQ=255)
                                          cATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGCCCGTGGGATCGGTCGGCTCCGTGGGATCGGTCggg                                                              >  2:4271390/1‑101 (MQ=255)
                                            tcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGGTc                                                            <  1:1239786‑M1/68‑1 (MQ=255)
                                            gcggtcgggtccgtgggatcggtcggatccgtgGGATCTGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc                                                            >  1:617669‑M1/34‑101 (MQ=255)
                                             cggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                                                           >  2:3236482‑M1/33‑101 (MQ=255)
                                              ggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTGGGGTCCg                                                          <  2:1286208‑M1/70‑1 (MQ=255)
                                              ggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                                          >  1:2547417‑M1/32‑101 (MQ=255)
                                                tcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTg                                                        <  1:3014105‑M1/72‑1 (MQ=255)
                                                tcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTg                                                        <  2:4502122‑M1/72‑1 (MQ=255)
                                                tcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                                        <  1:1643044‑M1/72‑1 (MQ=255)
                                                 cgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                       <  2:4585845‑M1/73‑1 (MQ=255)
                                                  gggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGtggt                                                      <  1:3088124‑M1/74‑2 (MQ=255)
                                                   ggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGa                                                     <  1:3673870‑M1/75‑1 (MQ=255)
                                                    gtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt                                                    >  1:2287150‑M1/26‑101 (MQ=255)
                                                     tccgggggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                   >  1:2452675‑M1/25‑101 (MQ=255)
                                                      gCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                                  >  1:2225158/1‑101 (MQ=255)
                                                         tCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                               >  1:3822269/1‑101 (MQ=255)
                                                          gggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                              <  1:1009700‑M1/82‑1 (MQ=255)
                                                             atcggtcggatccgtgGTATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGgatcggtcgggc                                           <  1:4418626‑M1/85‑3 (MQ=255)
                                                              tcggtcggtgccgtgGGATCGGTCGGGTCCGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                          >  2:4351057‑M1/16‑101 (MQ=255)
                                                              tcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATc                                          >  1:743508‑M1/16‑101 (MQ=255)
                                                              tccGTGGGGTCCGGCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATc                                          <  2:825592/99‑1 (MQ=14)
                                                                ggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGGCGGATCCg                                        <  1:187516‑M1/88‑1 (MQ=255)
                                                                ggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCg                                        <  1:2143941‑M1/88‑1 (MQ=255)
                                                                cGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg                                        <  1:2683422/101‑1 (MQ=255)
                                                                   cggatccgtgGGACCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                     >  1:4009373‑M1/11‑101 (MQ=255)
                                                                    ggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTggg                                    <  1:2828134/92‑1 (MQ=11)
                                                                     ggtccgggggatcggtcgggtccgtggggtcggtcgggtccgtgGGATCGGCCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                   <  1:3948674/93‑1 (MQ=11)
                                                                       tcggacggATCGGTCGGGTCGGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                 >  2:152885/1‑101 (MQ=255)
                                                                       tccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                                 >  2:3675182‑M1/7‑101 (MQ=255)
                                                                        ccgtgTGATCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATccgtggggtcc                                >  2:3511286‑M1/6‑97 (MQ=255)
                                                                        ccgtgGGATCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                >  2:1688214‑M1/6‑101 (MQ=255)
                                                                        ccgggggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcg                                <  2:3103215/96‑1 (MQ=11)
                                                                         cgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcgGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGGGGGAtcgg                               >  2:1994833/11‑101 (MQ=11)
                                                                          gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggt                              <  2:270790/101‑1 (MQ=14)
                                                                          gggGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                              >  1:522741‑M1/4‑101 (MQ=255)
                                                                           tgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGGTCCGTGGGATTGGTCGGGCCCGAGGGAtcttcc                             >  1:3428774/3‑97 (MQ=11)
                                                                           tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc                             <  1:2356379‑M1/99‑1 (MQ=255)
                                                                            gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                            <  1:1177001‑M1/100‑1 (MQ=255)
                                                                            cggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTACACTCGGTCGGGGCCGTGGGATGGGTCGGGTCAGTTGGAtcggtcg                            >  2:746382/1‑101 (MQ=11)
                                                                            cggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcg                            <  1:1327296/101‑1 (MQ=11)
                                                                             ggctccgtcgggtccgtcggatccgtcgggtccgtgggttcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                           <  2:3329374/95‑1 (MQ=11)
                                                                             ggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGGTCGGTCGGGGCCGTGGGATCGGTCGGCTCCGTGGGGtcggatgg                           >  2:1196290/1‑97 (MQ=11)
                                                                                tctgtggggtccgtgggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATc                        <  1:3880178/95‑1 (MQ=11)
                                                                                tcggtcgggtccgtgggatcggtcgggtccgtggggTCCGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATGGGTCGGATCCGTGGGGTCGGTCGGGTc                        >  1:457593/1‑101 (MQ=14)
                                                                                 cggtcgggtccgtgggatcggtcgggtccgtgggatCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGGcc                       >  1:565034/1‑101 (MQ=11)
                                                                                 cggtcggatccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggaaggatcc                       >  2:2380117/1‑93 (MQ=11)
                                                                                  gggcggatccgtgggatccgtggggtccgtgggatcggtCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCg                      <  2:2050128/98‑1 (MQ=11)
                                                                                   gtgggatccgggggatcggtcgggtccgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                     <  2:3759784/95‑1 (MQ=11)
                                                                                    tggggtccgcgggctcggtcgggtctgtgggttcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGtg                    <  2:1419354/99‑2 (MQ=14)
                                                                                    tcgggtccgtgggatcggtcgggtccgtcggatcggTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTg                    <  1:1993888‑M2/101‑2 (MQ=255)
                                                                                    tcgggtccgtcggatcggtcgggtccgtgggatcggTCGGGTCCGTCGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtg                    >  2:660527/1‑100 (MQ=11)
                                                                                      gggtccgtggggtccgtcggatccgtgggatcggtcGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGgg                  >  2:1868266/1‑101 (MQ=11)
                                                                                      gggtccgtgggatcggtcgggtccgtggggtccgtcGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg                  <  1:3736979‑M2/101‑4 (MQ=255)
                                                                                       ggtccgtgggatcggtcgggtccgtgggatcggtggGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTggga                 >  2:2452348‑M2/1‑97 (MQ=255)
                                                                                        atccgtgggatcggtcgggtccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgggat                <  1:4367514‑M2/100‑6 (MQ=255)
                                                                                        atccgtaggatcggtcgggtccgtaggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTGgggg                >  2:2120937/2‑99 (MQ=11)
                                                                                         tccgtgggatcgttcggctccgtgggatcgttcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc               <  2:3935923/101‑1 (MQ=11)
                                                                                         tccgtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatc               >  1:2922186‑M2/1‑95 (MQ=255)
                                                                                         gccgtgggatcgggcgggtccgtggggtccgtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc               <  1:2623748/100‑1 (MQ=11)
                                                                                          ccgtgggatcggtcggatccgtgggatcggtcggGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGTTCGTTGGGAGCGCTCAGATCCGTGGGAtcg              >  2:3845194/1‑101 (MQ=11)
                                                                                          ccgtgggagcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg              <  2:3583621/101‑1 (MQ=14)
                                                                                           cgtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGGGATCGGGCGGGTCCGTgggatcgg             <  2:302227‑M2/101‑9 (MQ=255)
                                                                                            gtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAGCCGTGGGGTCGGTCGGGTCGGTGGGAtgggg            >  1:4157606/1‑100 (MQ=11)
                                                                                             tgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc           >  1:924531/1‑101 (MQ=14)
                                                                                              gggatcggtcggttgcgggggctcggtcggTCCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTTCGTGGGATCGGTCGGATCCGTGGGAtcggtcg          <  1:320308/101‑1 (MQ=11)
                                                                                              gggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGGTCGGTCGGATCGGTGGGATCGGTCGGGTCGGTGGGGtcggtgg          >  2:3415217/1‑99 (MQ=11)
                                                                                              gggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTCGGTTCGGTCGGATCCGTGGGATCGGTCGGGTCCGGGGGAtcgggcg          >  2:4150676/1‑101 (MQ=14)
                                                                                               ggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGGTCGGTCGGATCCGTGGGAtcggtcgg         >  2:4400714/1‑101 (MQ=14)
                                                                                                gatcggtcgggaccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTGGGATCAGCCGGGTCCGTTGGGtcgggctcg        >  2:175871/1‑96 (MQ=11)
                                                                                                 atcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcgggt       <  2:1599610‑M2/101‑15 (MQ=255)
                                                                                                  tcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatc      <  2:2572114/101‑1 (MQ=11)
                                                                                                    ggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTGGGCTCGgtcgggtcgg    >  2:908137/1‑99 (MQ=11)
                                                                                                    ggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGCCGGGTCCGTGGGGTCGgtccgatccg    >  2:554600/1‑101 (MQ=11)
                                                                                                    ggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccg    <  1:1084279‑M2/101‑18 (MQ=255)
                                                                                                    |                                                                                                      
CCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTC  >  NC_014034/221879‑222081

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