Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,090 G→C 8.8% P189P (CCC→CCG RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0900G→C 8.8% 62.0 0.0 ‑0.0359P189P (CCC→CCGRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (16/16):  ref base (160/167):  total (176/183)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.95e-01

Only 100 of 489 total aligned reads displayed.

TGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAGGCCGGCCAGACCGG  >  NC_014034/221990‑222190
                                                                                                    |                                                                                                    
tgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATGGGTCGGGTCCGTGGGATCCGTGGGGTCCTTGGGGTCGGTg                                                                                                      >  2:4529966/1‑101 (MQ=11)
tgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTg                                                                                                      >  1:4106899/1‑101 (MQ=14)
 gggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTgg                                                                                                     <  1:2238764/101‑1 (MQ=11)
 aggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTGGGATCCGTCGGGTGCGTCGGGCCCGGgg                                                                                                     >  2:4203360/2‑101 (MQ=11)
  gggtccgtgggatcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGGCCGTggg                                                                                                    <  2:1801691‑M1/74‑1 (MQ=255)
  gggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTggg                                                                                                    <  2:2549729/95‑1 (MQ=11)
  gggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGGCCGTGGGATCGGTCGGGGCCGTCGGGTCGGTggg                                                                                                    >  1:684943/7‑101 (MQ=14)
  ggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTcgg                                                                                                    <  2:3524235/101‑1 (MQ=11)
   ggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGa                                                                                                   >  2:4139590/6‑101 (MQ=11)
   gatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGa                                                                                                   >  2:1885251/1‑101 (MQ=14)
    gtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGCCGGGCCCGTGGGATCGGTGGGGGCCGGCGGATCCGGggggt                                                                                                  >  2:396353/5‑99 (MQ=2)
       ggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAGCCGGGGGGTCCGTCGGGTCCGGGGGAtcgg                                                                                               >  2:2374996/1‑101 (MQ=11)
        gtgggatcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGgtcggg                                                                                              <  2:2120756‑M1/80‑2 (MQ=255)
        gtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGAtccgg                                                                                              >  1:3697852/1‑98 (MQ=11)
        gtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccgt                                                                                              <  1:2750333/101‑1 (MQ=11)
         tgggatcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGGGGGAtgggtc                                                                                             <  1:2175895‑M1/81‑1 (MQ=255)
         tgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGGCCGGCGGGCCCGTGGGAccgggc                                                                                             >  2:3373023/3‑99 (MQ=11)
         tcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGGCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtc                                                                                             <  1:4643608/101‑1 (MQ=14)
         tcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgtg                                                                                             >  2:2628732/1‑100 (MQ=11)
         tcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgtg                                                                                             >  2:1534830/1‑100 (MQ=11)
          cggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcg                                                                                            <  1:3628423/101‑1 (MQ=14)
          cggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccgtgg                                                                                            >  1:3312242/1‑99 (MQ=11)
          cggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGGCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGGTCCGTGGGGtccgtcg                                                                                            >  1:1041686/1‑101 (MQ=11)
           gggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCCGTGGGGTCCGTCGGGTCCGTGGGGtcggtggg                                                                                           >  2:4073682/4‑101 (MQ=11)
            ggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcgggcggg                                                                                          >  2:1046563/3‑100 (MQ=14)
            gatcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGa                                                                                          <  1:3097827‑M1/84‑1 (MQ=255)
            gATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTGGGATCCGTCGGGTCCGTGGGCtcggtgggc                                                                                          >  2:3664127/1‑100 (MQ=11)
            gATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgtgggc                                                                                          >  2:2052796/1‑100 (MQ=11)
             atcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGAt                                                                                         <  2:674141‑M1/85‑1 (MQ=255)
             aTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGACCCGTGGGAtccggggggc                                                                                         >  2:4330189/1‑93 (MQ=11)
              tcggttggATCCGTCGAGGCCGTCAGATCGGTCGGATCCGTGGGATCTGTCGGGTCCGTGGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTc                                                                                        <  2:759521/95‑1 (MQ=14)
              tcggtcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATc                                                                                        <  1:2330835‑M1/86‑1 (MQ=255)
              cccGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATc                                                                                        <  2:2017560/100‑1 (MQ=11)
                cGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTGGGGTCCg                                                                                      >  2:2178612/1‑101 (MQ=14)
                  tcggatccgtgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATccggg                                                                                    <  1:1710406‑M1/90‑3 (MQ=255)
                    ggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGGCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTcgg                                                                                  >  2:1952719/1‑101 (MQ=14)
                     gATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTAGGATCCGTGGGGTCCGTGGGATCGGTCGGGGCAGTGGGa                                                                                 >  2:3443123/1‑101 (MQ=14)
                        cggtcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcg                                                                              <  2:2719412/101‑1 (MQ=14)
                         ggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGGTCCGTGGGGTCCGTCGGGTCCGGGGGAtggg                                                                             >  2:10076/1‑101 (MQ=11)
                          gtcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCCTGGGATCGGTCGGGTCCGTCGGAtcggt                                                                            >  1:366191/1‑101 (MQ=14)
                           tgGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                           <  1:3390565‑M1/99‑1 (MQ=255)
                           tcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtc                                                                           <  2:543712/101‑1 (MQ=14)
                           ncggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGGCGGGGCCGTGGGATCGGGCGGGTCCGTCGGGtcgggc                                                                           >  1:94156/2‑99 (MQ=14)
                            cggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGCCCGTCGGATCCGTGGGGCCCGTCGGGTCCGTGGGATCGGGCGGGTGCGTCGGAgcggtgg                                                                          >  2:4612280/1‑99 (MQ=18)
                               gTCCGTGGCATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGt                                                                       <  2:4685993/101‑1 (MQ=255)
                                tCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTc                                                                      <  2:2905214/101‑1 (MQ=17)
                                ctcttccgATCTGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc                                                                      <  1:848243/94‑1 (MQ=11)
                                   gtcggATCGGTCGGATGCGTGGGATGGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTGCGTGGGAGCGGTCGGGTCCGTCGGATCGGTCGGGTCCGt                                                                   <  1:2099285/98‑1 (MQ=14)
                                   gTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTTCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGCCGGGTCCgg                                                                   >  2:3134924/1‑100 (MQ=255)
                                    tcggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCCGTGGGGTCCGTCGGGTCGGTGGGATCGGTCGGGTCCGTCGGTGCGGTGGGGTCCGTg                                                                  >  1:3747356/3‑101 (MQ=17)
                                    gcggctcgggcggATCGGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTg                                                                  <  2:2253517/96‑1 (MQ=14)
                                     cggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGGCCGTgg                                                                 >  2:1860473/2‑101 (MQ=17)
                                      ggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTggg                                                                >  1:2927499/1‑101 (MQ=18)
                                        aTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggat                                                              <  1:4330871/101‑3 (MQ=14)
                                          cggtcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtg                                                            >  1:1330672/1‑101 (MQ=255)
                                          cggtcggATCCGTGGGATCGGGCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGGCGGGTGCGTCGGAACGGTGGGGTGCATgggatcg                                                            >  2:1588049/1‑97 (MQ=9)
                                            gtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggt                                                          >  1:2016394/1‑95 (MQ=11)
                                             gcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtggtg                                                         <  1:1863376/100‑1 (MQ=255)
                                               gggtCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTggg                                                       >  2:3336104/4‑101 (MQ=255)
                                                ggtCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGa                                                      >  1:1000596/3‑101 (MQ=255)
                                                 gtCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAt                                                     <  1:4597964/100‑1 (MQ=255)
                                                  tCCTTGGGATCCGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtc                                                    <  2:4436193/101‑1 (MQ=255)
                                                    cGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGGCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcgg                                                  >  2:32439/1‑101 (MQ=255)
                                                      tGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGGCCGTGGGGGTGGGGGGAtcggtc                                                >  1:3570542/1‑101 (MQ=255)
                                                       gggATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGCGGGAtcggtcg                                               >  1:3959900/1‑101 (MQ=255)
                                                          aTCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGCCGGAc                                            >  2:4754525/1‑101 (MQ=255)
                                                               tcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTg                                       >  1:4360516/1‑101 (MQ=255)
                                                               gggggTCCGTGGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTg                                       <  1:1138405/99‑1 (MQ=255)
                                                                cggGTCCGTAGGATCGGTCGGATCCGTCGGATCCGTCGGATCGGTAGGATCCGTGGGATCGGTCGGGTCTGTCGGGCTCGTTGGTTCGGTCGGAcccttgg                                      >  2:2973956/1‑96 (MQ=18)
                                                                 gggTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGAGCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTc                                     <  1:2723061/101‑1 (MQ=255)
                                                                  ggTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGCCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACAGCTGTcg                                    >  2:1344377/1‑101 (MQ=255)
                                                                      cGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgccgc                                >  1:189225/1‑101 (MQ=255)
                                                                       gTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgccgcc                               >  1:2636801/1‑101 (MQ=255)
                                                                         aGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgccgccgc                             >  2:2466495/1‑101 (MQ=255)
                                                                          ggATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCgccgccgcc                            >  2:3802092/1‑101 (MQ=255)
                                                                          ggATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCgccgccgcc                            >  1:2955389/1‑101 (MQ=255)
                                                                            aTCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGCCGCCGCCGCCgg                          >  2:1464031/1‑101 (MQ=255)
                                                                               ggtcggGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACa                       >  1:3821443/1‑101 (MQ=255)
                                                                               ggtcggGTCCGGCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACa                       <  2:1235896/101‑1 (MQ=255)
                                                                                gtggggtcggtcgggtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAg                      <  1:1238091/98‑1 (MQ=255)
                                                                                     gtccgtcgggtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCa                 >  1:3248954/1‑101 (MQ=255)
                                                                                          tcgggtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAggccg            <  2:835388/101‑1 (MQ=255)
                                                                                          tcggatccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtg            <  1:1465864‑M2/101‑53 (MQ=255)
                                                                                          tcggatccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtg            >  1:2052730‑M2/1‑49 (MQ=255)
                                                                                            ggatccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtggg          >  2:574633‑M2/1‑47 (MQ=255)
                                                                                             ggtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAggccggcc         <  2:1146103/101‑1 (MQ=255)
                                                                                             ggtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAggccggcc         >  1:2438802/1‑101 (MQ=255)
                                                                                              gtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAGGCCGGCCa        <  1:1545208/101‑1 (MQ=255)
                                                                                               tcggtGGGGTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAGGCCGGCCAg       <  1:3053704/98‑1 (MQ=255)
                                                                                               tccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAGGCCGGCCAg       >  2:2749482/1‑101 (MQ=255)
                                                                                                ccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcg      >  1:4506249‑M2/1‑43 (MQ=255)
                                                                                                ccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAGGCCGGCCAGa      <  2:1103814/101‑1 (MQ=255)
                                                                                                  gtGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGatcggtcgggtccgtggggtccgtcgggtccgtgggatcggtcggatccgtgggatcggt    <  2:3723545‑M2/101‑61 (MQ=255)
                                                                                                  gtGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcgggtccgtgggatgggc    <  2:3068441‑M2/101‑61 (MQ=255)
                                                                                                  gtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatgggc    <  2:2329012‑M2/101‑61 (MQ=255)
                                                                                                  gtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggt    >  1:351102‑M2/1‑41 (MQ=255)
                                                                                                  gtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAGGCCGGCCAGAcc    <  1:3464984/101‑1 (MQ=255)
                                                                                                  gtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAGGCCGGCCAGAcc    <  1:4754781/101‑1 (MQ=255)
                                                                                                    gggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcg  >  2:3860452‑M2/1‑39 (MQ=255)
                                                                                                    gggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcg  >  1:1966022‑M2/1‑39 (MQ=255)
                                                                                                    |                                                                                                    
TGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAGGCCGGCCAGACCGG  >  NC_014034/221990‑222190

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