Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,009 G→A 6.2% P216P (CCC→CCT RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0090G→A 6.2% 0.0 ‑1.9 ‑0.73684P216P (CCC→CCTRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (117/103):  ref base (1588/1728):  total (1780/1904)
Fisher's exact test strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.82e-01

Only 100 of 4222 total aligned reads displayed.

CCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGG  >  NC_014034/221909‑222109
                                                                                                    |                                                                                                    
tcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTg                                                                                                      <  1:4722601‑M1/54‑1 (MQ=255)
   gatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGa                                                                                                   <  1:5998284‑M1/57‑1 (MQ=255)
        gtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                                              >  1:6858928‑M1/40‑101 (MQ=255)
         ggggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                                             >  2:1369860‑M1/39‑101 (MQ=255)
            ggtccgtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCggg                                                                                          <  2:6419023‑M1/66‑1 (MQ=255)
               cggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                                                       <  1:1312118‑M1/69‑1 (MQ=255)
               cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                                                       <  1:8162246/101‑1 (MQ=255)
                    ggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                                                                                  >  1:5633998‑M1/28‑101 (MQ=255)
                    ggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGGccgtcgg                                                                                  <  2:4076522‑M1/74‑4 (MQ=255)
                      gtccgtgggatcggtcgggtgcgtgGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                                >  1:3295990‑M1/26‑101 (MQ=255)
                      gtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt                                                                                >  1:5013910‑M1/26‑101 (MQ=255)
                       tccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                                                                               <  2:5674911‑M1/77‑1 (MQ=255)
                        ccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                                                              >  1:6697488‑M1/24‑101 (MQ=255)
                         cgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                                                             >  1:6205123‑M1/23‑101 (MQ=255)
                          gtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                                                            <  1:1137202‑M1/80‑1 (MQ=255)
                           tgcgATGCGTCGCCTCCCTGGCATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                           <  2:8684869/98‑1 (MQ=11)
                           ggggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                           >  1:6394995‑M1/21‑101 (MQ=255)
                             gtatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtccgtcgg                                                                         >  2:5227116‑M1/19‑101 (MQ=255)
                                tGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                                                      <  2:5889387/101‑1 (MQ=255)
                                 cggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                                                                     <  1:1459990‑M1/87‑1 (MQ=255)
                                  ggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGtgggatcggtcggcttcg                                                                    <  2:3748121‑M1/88‑6 (MQ=255)
                                   gTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                                                   <  2:6704780/101‑1 (MQ=11)
                                     gggATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg                                                                 >  1:317019/1‑99 (MQ=1)
                                      gggtcgggcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTcgg                                                                <  1:1000386/98‑1 (MQ=11)
                                       ggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtgggg                                                               <  2:1867307‑M1/93‑2 (MQ=255)
                                        gtcggtgggGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                                              <  2:6598010/100‑1 (MQ=11)
                                        gtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCgtggggt                                                              >  2:1552771‑M1/8‑99 (MQ=255)
                                         tccgtgGGATCGGGCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                             >  2:8345471‑M1/7‑101 (MQ=255)
                                          ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGACCCGTGGGATCGGTCGGATCCGTGGGAacg                                                            >  1:4278904/6‑101 (MQ=11)
                                            gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGCTCGGGTCCGCCGGAtcgcg                                                          >  1:3497310/1‑99 (MQ=14)
                                               ggatcggtcgggtccgtgggatcggtcgggtccgtcGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGCGGAtcggtcgg                                                       >  2:7017572/1‑101 (MQ=11)
                                                 atcggtcggatccgtgggatcggtcggatccgtgggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggggggac                                                     >  2:687022/1‑95 (MQ=11)
                                                  tcggtcgggtcggtgggatcggtcgggtccgtgggaTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatc                                                    >  2:8828219‑M2/1‑87 (MQ=255)
                                                  tcggtcggatccgtgggatcggtcgggtccgtgggaTCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATc                                                    <  1:7277098/101‑1 (MQ=11)
                                                  tcggccgtgtccgtgggcgcggtcgggtccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                                    <  2:2534088/101‑1 (MQ=14)
                                                   cggtcgggtccgtgggatcggtcgggtccgtggggtCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                   >  1:7670097/1‑101 (MQ=11)
                                                   cggtcgggtccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATcc                                                   <  2:1023221/101‑1 (MQ=14)
                                                   cggtcgggtccgtcggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                   <  2:8068186/101‑1 (MQ=11)
                                                   cggtcggatccgtgggatcggtcggatccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGACCGGTCGGATcg                                                   >  2:1196867/1‑100 (MQ=11)
                                                   ccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGCGCCGTCGGGCCGGTCGAGTcc                                                   >  2:3440298/3‑101 (MQ=11)
                                                    ggtcgggtccgtgggatcggtcgggtccgtggggtcCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGCTCGGATCCGTGGGATCGGTCGGGTCCg                                                  >  1:2811041/1‑101 (MQ=14)
                                                    ggtcgggtccgtcggatcggtcgggtccggggggtcCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                                  <  2:3596279/101‑1 (MQ=11)
                                                     gtcggatccgtgggatcggtcgggtccgtgggatcgGTCGGATCCGTGGGATGGGTCGGGTACGTGGGTTCGGTCGGGTGCGTGGGGTCGGTTGGGTCCgg                                                 >  2:2345422/1‑100 (MQ=11)
                                                     ggcggagccgtggggtcggtcggagccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGt                                                 <  1:3812411/94‑1 (MQ=11)
                                                       cgggtccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTGGGATCCGTGGGATCGGTCGGGTCGGtgg                                               >  2:4713120/1‑99 (MQ=14)
                                                       cgggtccgtcggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGCGGGATCGGTCGGGGCCGTGGGATCGGTCGGGGCCGtgg                                               >  1:4109767/1‑99 (MQ=11)
                                                        gggtccgtggggtcgttcgggtgggtgggattggtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCTGGTCCGTGGGATCGGTCGGGTCCGTGgg                                              <  1:5216718/101‑1 (MQ=11)
                                                        gggtccgtgggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTCGGGGGGGTCCGTCGGATCCGGGGGAgcggtcggagccggggg                                              >  2:2692692/1‑92 (MQ=14)
                                                        gggtccgtgggatcggtcgggtccgtcggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGACCGTGGGATCGGGCGGGTGCggggg                                              >  2:4808635/1‑97 (MQ=11)
                                                         ggtccgtgggatcggtcgggtccgtggggtccgtcgGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCgagggg                                             >  2:2202553/1‑96 (MQ=11)
                                                         ggtccgtgggatcggtcgggtccgtggggtccgtcgGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggga                                             >  1:3939459‑M2/1‑97 (MQ=255)
                                                          gtccgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggat                                            <  1:3025254‑M2/101‑6 (MQ=255)
                                                          gtccgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggat                                            >  2:2054432‑M2/1‑96 (MQ=255)
                                                          atccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGGTCGGTCGGGGCCGTGgggg                                            >  2:4549963/2‑99 (MQ=11)
                                                           tccgtgtgatcggtcgggtgcgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatc                                           >  2:2848927‑M2/1‑95 (MQ=255)
                                                           tccgtggggtcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcggtcgggtc                                           >  1:8647475‑M2/1‑78 (MQ=255)
                                                           tccgtgggatcggtctggtccgtgggatcggtcggCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatc                                           >  2:8212027‑M2/1‑95 (MQ=255)
                                                           gccgtgggatcggtcgggtccgtgggatcgntcgggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatc                                           >  1:2143356‑M2/2‑95 (MQ=255)
                                                            ccgtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcg                                          <  1:8297000‑M2/101‑8 (MQ=255)
                                                              gtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggt                                        <  2:8518788‑M2/101‑10 (MQ=255)
                                                               tgggatcggtcgggcctctgggctcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                       <  2:5370274/101‑1 (MQ=11)
                                                                ggggtcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                      <  1:2295003/101‑1 (MQ=11)
                                                                   atccgtgggatcggtcgggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggt                                   <  1:2620791/101‑1 (MQ=14)
                                                                    tcggtcgtgtgcgtggggtcggtcggATCCGGGGGACCCGTCCGCTCCGTGGGACCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtc                                  <  2:3241845/101‑1 (MQ=14)
                                                                     gggtcgggttcgggagagcggttggGTTGTTGGGATCGGGCGGGTCCGTGGGATCAGGCGGGTCCGTGGGCTTGGTCGGGTCCGTGGGATCGgtcgggtcc                                 <  1:4581209/100‑1 (MQ=2)
                                                                     cggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtcc                                 <  1:5519640‑M2/101‑17 (MQ=255)
                                                                     cggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtcc                                 >  1:1765588/1‑101 (MQ=11)
                                                                     cggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtcc                                 <  1:4846933‑M2/101‑17 (MQ=255)
                                                                     cggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtcc                                 <  1:804927/101‑1 (MQ=11)
                                                                      ggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccg                                <  1:7136710‑M2/101‑18 (MQ=255)
                                                                      ggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggagcggtcgggtccg                                <  2:4991495‑M2/101‑18 (MQ=255)
                                                                       gtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgt                               <  1:7483574‑M2/101‑19 (MQ=255)
                                                                       gtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgt                               <  2:7872126/101‑1 (MQ=11)
                                                                       gtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgt                               <  1:583194‑M2/101‑19 (MQ=255)
                                                                         cgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtcg                             <  1:914526/101‑1 (MQ=11)
                                                                          ggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtgtggtccgtcgg                            >  1:6537041/4‑101 (MQ=11)
                                                                           ggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatccgtgggg                           >  2:7314862/1‑101 (MQ=11)
                                                                               cgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgggggatcgg                       >  1:6978756/1‑92 (MQ=11)
                                                                               cgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggaccggtcgggtccgtgggatcggtcgggtccg                       <  1:8798219‑M2/101‑44 (MQ=255)
                                                                                gtggggtccgtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgt                      <  1:897693/101‑1 (MQ=11)
                                                                                  ggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTggggtccgtggggtccgtcg                    >  2:3920114/1‑99 (MQ=11)
                                                                                  gggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGGTCCTTgg                    >  2:42572/1‑101 (MQ=11)
                                                                                   gggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTcgg                   <  1:1599137/95‑1 (MQ=11)
                                                                                   ggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCGGGCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgg                   <  2:7539845‑M2/101‑31 (MQ=255)
                                                                                   ggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTggg                   <  2:3094871/101‑1 (MQ=11)
                                                                                      tggggcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtc                <  1:7823533/96‑1 (MQ=14)
                                                                                      tcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtaggatc                <  2:1724932‑M2/101‑51 (MQ=255)
                                                                                       gggtcgggTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtcc               >  1:2455537‑M2/2‑67 (MQ=255)
                                                                                       cggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcg               <  1:5355556/101‑1 (MQ=14)
                                                                                       ccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtcc               <  2:6253976/99‑2 (MQ=11)
                                                                                         gtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGGCCGTGGGAtcggg             >  2:6076951/1‑100 (MQ=255)
                                                                                         ggcggATCCGTGGGATCGGTCGTATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgt             <  2:7724807/99‑1 (MQ=11)
                                                                                           cggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGGTCGGTCGGGTGCGCCGGGtgggtcg           >  2:4412853/1‑101 (MQ=14)
                                                                                           cggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcggatccgtgggatcggtcggatccgtgg           <  1:8448370‑M2/101‑39 (MQ=255)
                                                                                            ggatgggTGGGCTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggtcgg          <  2:9037695/95‑1 (MQ=18)
                                                                                             gATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggaacgggcgggtccgtggga         <  2:5554627‑M2/101‑41 (MQ=255)
                                                                                              aTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAt        >  1:4431406/1‑101 (MQ=11)
                                                                                                ccGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcg      <  2:7154482‑M2/101‑44 (MQ=255)
                                                                                                   tGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtc   >  1:3645376‑M2/1‑55 (MQ=255)
                                                                                                   ggggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtc   >  1:1957909‑M2/2‑55 (MQ=255)
                                                                                                    |                                                                                                    
CCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGG  >  NC_014034/221909‑222109

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 41 ≤  ATCG