Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,069 G→C 27.1% T196T (ACC→ACG RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0690G→C27.1% 0.0 ‑0.8 ‑0.0690T196T (ACC→ACGRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (89/99):  ref base (253/246):  total (342/348)
Fisher's exact test strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.72e-01

Only 100 of 978 total aligned reads displayed.

CCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACA  >  NC_014034/221969‑222169
                                                                                                    |                                                                                                    
ccgttggctcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                                                                                      <  1:2531574/101‑1 (MQ=11)
  gtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGACCCGTGGGGTCGGTCGGATCCGTGGGCTCGGGCGGAGCCGTGGGGtcggg                                                                                                    >  1:4875470/1‑100 (MQ=11)
  gtgggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCGGTCGGGTCGGCGGGAtcggt                                                                                                    >  2:838623/1‑101 (MQ=14)
  gtgggagcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTTGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                                                    <  2:4433639/101‑1 (MQ=11)
   tgggatcggtcgggtccttggtttcggttggGTCCGTGGGATCGGGCGGGTCTGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc                                                                                                   <  1:1315716/101‑1 (MQ=11)
   tgggatcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGAACGTCGGATCGGCCGGGATCGTGGGAtcggtc                                                                                                   >  1:2949978/1‑101 (MQ=11)
   tgggatcggtcggatccgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGGGGAtcgggc                                                                                                   >  1:3952069/1‑99 (MQ=11)
   cgggatcggtcggttccgtggtatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGAtccgtg                                                                                                   <  2:2334551/100‑2 (MQ=11)
    gggatcgggcggatccgggggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGAtcggtcg                                                                                                  <  2:4821668/101‑1 (MQ=14)
     ggatcggtcgggtccgtggggtgcgtaggAGCCGGGGGATCGGTCGGATCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                                 <  1:715557/101‑1 (MQ=11)
     ggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCGGTGGGGtcggtagg                                                                                                 >  2:3854849/1‑101 (MQ=11)
     ggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCGGTCGGGTCCCTCGGGtcggtcgg                                                                                                 >  2:1665930/1‑101 (MQ=14)
     ggagtggacgggtccgtggactggtccggGTCCGAGGGATCGGTCGGGTCCGTGGGATAGGTCGGATCCGTGCGATAGGTCGGGTCCGTGGGGtcggtggg                                                                                                 >  2:180883/9‑101 (MQ=2)
       gtcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggt                                                                                               <  2:4136321/100‑1 (MQ=11)
         cggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGCCCGTGGGATCGGTCGGGTCCGTGGGATCGGCCGGGGCCGTGGGATCGgtcgggtcc                                                                                             >  1:1125334/1‑101 (MQ=11)
           gtcgggtccgtggggtccgtcggATCCGTGGGCTCGGTTGTATCCGTGGGATCTGTCGGTTCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatccgt                                                                                           <  2:3751862/101‑1 (MQ=14)
           gtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgg                                                                                           >  2:9115/1‑100 (MQ=14)
           gtcggggccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgt                                                                                           <  1:950525/101‑1 (MQ=11)
            tcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtc                                                                                          <  2:718233/101‑1 (MQ=11)
             cgggtccgtgggggccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatccgtgg                                                                                         <  1:4104282/101‑3 (MQ=11)
             cgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCTTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtgg                                                                                         <  1:3384091/101‑3 (MQ=11)
             cggatccgttggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcg                                                                                         <  1:4301436/101‑1 (MQ=11)
             cggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcg                                                                                         >  2:2687683/1‑101 (MQ=255)
              ggtctcgtgggggccgtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtggg                                                                                        <  1:1865327/96‑1 (MQ=11)
              gggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGGGGGGGCGGCCGGGTGCGTAGGAggggtcggatgggtggg                                                                                        >  2:5132804/1‑92 (MQ=2)
              gggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtggg                                                                                        >  2:3283505/1‑101 (MQ=11)
                atccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcggat                                                                                      <  2:5719106/100‑3 (MQ=255)
                atccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGgtcggatccgtggggt                                                                                      >  1:3228632/2‑101 (MQ=11)
                  ccgtgggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggtcc                                                                                    <  1:100380/101‑1 (MQ=14)
                   cgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGCTCGgtcgggtccgtgggagatc                                                                                   >  2:1865217/1‑96 (MQ=11)
                    gtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTcgggtccgtcggatcggt                                                                                  >  1:4631241/1‑101 (MQ=11)
                    gtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTggggtccgtcgggtccgt                                                                                  >  1:4224916/4‑101 (MQ=14)
                     tgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgtc                                                                                 <  2:772882/101‑2 (MQ=11)
                     gagggtcggtcggGTCCGTTGGGTCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTg                                                                                 <  1:2885480/96‑1 (MQ=14)
                     aaggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGAATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTcgggtccgtcggatcggtc                                                                                 <  2:3943213/99‑2 (MQ=11)
                      cggttcggttggGTTCGTGGGATCGGTCTGGTCCGTCGGCTCGGGGGGATCCGTGGGATCTGTCGGGTCAGTGGGATCGgtcgggtccgtcggatcgggcg                                                                                <  1:3304972/97‑6 (MQ=2)
                      aggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGGGCCGTCGGGTCCGTgg                                                                                >  2:4635132/2‑101 (MQ=14)
                      aggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGtggggtccgtcggggccggtg                                                                                >  2:2033215/2‑98 (MQ=11)
                       ggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTggg                                                                               <  2:2924732/101‑1 (MQ=14)
                       ggatcggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGCGGATCGGTCGGCTAGGGGGGCTCGGTCGGGTGCGTGGGGTGCGTcgg                                                                               >  2:5555418/1‑101 (MQ=2)
                        gatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGa                                                                              <  1:2119136/101‑1 (MQ=14)
                         gtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGac                                                                             >  1:5756475/5‑100 (MQ=11)
                         gtccgggggATCCGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAt                                                                             <  1:310716/94‑1 (MQ=17)
                          tcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGtc                                                                            <  1:3505051/101‑1 (MQ=11)
                           ccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcg                                                                           <  2:254742/96‑1 (MQ=14)
                           ccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcg                                                                           >  2:309852/6‑101 (MQ=14)
                             gtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggt                                                                         <  2:2204645/101‑1 (MQ=14)
                               gggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggtcg                                                                       <  2:4319356/100‑1 (MQ=17)
                               cggATCCGGGTGAGCGGACTGGTGGGTCGGATCCGTGGGGTCCGTCGGGGCCGTGGGATCGGTGGGGGCCGTCGGATCGGTCGGGTACGTGGGAtcggtcg                                                                       >  2:2805094/1‑101 (MQ=12)
                                gggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGGTCCGTCGGAtccgtcgg                                                                      >  1:2770202/4‑101 (MQ=11)
                                gggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcgg                                                                      >  1:2194280/4‑101 (MQ=14)
                                ggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGAtccgtggg                                                                      <  2:5404597/101‑1 (MQ=1)
                                 ggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcggg                                                                     >  2:5340265/3‑100 (MQ=14)
                                 gATCCGTGGGGGCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGa                                                                     <  1:1731677/101‑1 (MQ=11)
                                  gtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcgggt                                                                    <  2:1584611/100‑3 (MQ=11)
                                  gtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcgggg                                                                    >  1:2378664/2‑99 (MQ=14)
                                  aTCCGTGGGATCGGTCGTATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAt                                                                    <  2:1273984/101‑1 (MQ=11)
                                   tCCGTGGGAGCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATc                                                                   <  1:5027443/101‑1 (MQ=11)
                                     cGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCg                                                                 <  2:564263/101‑1 (MQ=14)
                                      gtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGt                                                                <  1:5481442/98‑1 (MQ=18)
                                      gtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGt                                                                >  2:4257304/4‑101 (MQ=17)
                                       ggggATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTg                                                               <  1:4541385/100‑1 (MQ=14)
                                        gggATCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                              <  2:5182631/101‑1 (MQ=14)
                                           aTCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAGCGGTCGGATCCGTGGGAt                                                           >  2:212497/1‑101 (MQ=18)
                                            tccgtggggTCCGGGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                          <  2:611920/95‑1 (MQ=17)
                                              cgtggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                                        <  2:1166784/97‑1 (MQ=255)
                                                 ggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGGCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                                     >  1:4957986/2‑101 (MQ=255)
                                                 cggGTCCGTGGGAACGGTCGGATCCGTGGGATCGGGCGGGTCGGGGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                     <  1:2009816/101‑1 (MQ=17)
                                                  gggTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                    >  1:1247627/1‑101 (MQ=17)
                                                    gggtgttggATCAGTCGGATCCGCGGGATGGGGCGGGTCCGTCGGATCCGGGGGGTCCGCCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGt                                                  <  1:3820812/94‑1 (MQ=9)
                                                    gTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGt                                                  >  1:3059930/1‑101 (MQ=17)
                                                      ccGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                >  1:5651628/1‑101 (MQ=255)
                                                      ccGGGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                <  1:1965922/101‑1 (MQ=17)
                                                       cGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGGTCGGTCGGGTCCGGCGGATCGGTCGGGTCCg                                               >  2:1078006/1‑101 (MQ=17)
                                                        gtcggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGt                                              <  1:852359/98‑1 (MQ=17)
                                                         tcggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTg                                             <  1:3682676/99‑1 (MQ=17)
                                                         tGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGta                                             >  1:5403358/1‑100 (MQ=17)
                                                          gggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgg                                            >  1:995589/1‑101 (MQ=255)
                                                          cggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgg                                            >  2:193980/2‑101 (MQ=17)
                                                           gtcgcggtcggGCCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg                                           <  2:3548701/97‑1 (MQ=255)
                                                            gATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggg                                          >  2:4035425/1‑101 (MQ=17)
                                                            gATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTggga                                          >  2:3194633/1‑100 (MQ=17)
                                                               ctgtcggGTCGGTGGGATCGGTCGGGTCCGTGGGATCCGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcg                                       <  1:5218902/99‑5 (MQ=14)
                                                                ggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcgg                                      >  1:2681783/1‑96 (MQ=14)
                                                                ggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcgg                                      >  2:3050703/1‑96 (MQ=14)
                                                                ggtcgcGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcgg                                      <  2:2452779/101‑6 (MQ=14)
                                                                 gtcggGTCCGTGGGATCGGTCGGGTCCGGAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtggt                                     <  2:1324330/101‑1 (MQ=255)
                                                                 gtcggGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTcgggtccgt                                     <  1:4000841/101‑6 (MQ=14)
                                                                 gtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTaggatcggt                                     >  1:1362469/1‑95 (MQ=14)
                                                                    ggATCCGTGGGATCGGTCGGGGCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggtcgg                                  <  2:4518888/101‑10 (MQ=11)
                                                                       tCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtc                               <  2:1480128/101‑1 (MQ=255)
                                                                       tCCGTCGGAACGGTCGGGTCGGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtc                               <  1:536096/101‑1 (MQ=255)
                                                                        gcGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcg                              <  1:1522243/100‑1 (MQ=255)
                                                                          gTGCGATCGGTCGGGTGCTTAGGATCGGTCGCGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggt                            <  2:3415400/101‑1 (MQ=21)
                                                                                 cggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCg                     >  1:1627923/1‑101 (MQ=255)
                                                                                  ggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGGCGGACCGc                    >  2:5280953/1‑101 (MQ=255)
                                                                                   gtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCt                   >  2:4980564/1‑101 (MQ=255)
                                                                                   gtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCt                   >  1:622402/1‑101 (MQ=255)
                                                                                         tttgTAGGACCTGTGGGTTCCGGCGGGTCCGTGGGATCGGTCTGATCCGGGGGATCTGTCTGGTCCGTGGGGCTGGTTGGATCGGTCGGACCGCTGTcgcc             <  1:4764048/98‑1 (MQ=18)
                                                                                                gATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGccgccgccg      <  1:4459439/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
CCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACA  >  NC_014034/221969‑222169

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 40 ≤  ATCG