Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,493,971 Δ1 bp 2.6% intergenic (+19/‑223) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,493,9710G→. 2.6% 7.4 ‑0.5 ‑0.1189intergenic (+19/‑223)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/3):  ref base (46/137):  total (49/140)
Fisher's exact test strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.29e-01

Only 100 of 214 total aligned reads displayed.

TGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA.GGGGGG.GCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGAGC  >  NC_014034/2493876‑2494064
                                                                                                 |                                                                                              
tGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCg                                                                                            <  1:151217/101‑1 (MQ=39)
 gATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGc                                                                                           >  1:5492755/1‑101 (MQ=39)
  aTGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCa                                                                                          <  1:3569915/101‑1 (MQ=39)
   tGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GGGCCGCCCCCCCTCCGcca                                                                                         >  2:4738769/1‑99 (MQ=255)
   tGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCaa                                                                                         <  2:3592959/101‑1 (MQ=38)
    ggCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAg                                                                                        <  1:1760994/101‑1 (MQ=38)
    ggCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAg                                                                                        >  2:341033/1‑101 (MQ=38)
     gCCGAGCTGCCCCCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGc                                                                                       <  2:4012335/101‑1 (MQ=38)
       cGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCgg                                                                                     >  1:5414882/1‑101 (MQ=38)
           cTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGAACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTc                                                                                 <  1:2134887/101‑1 (MQ=37)
            tGCCCCCCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTcc                                                                                <  2:587110/101‑1 (MQ=37)
            tGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCCGCCCCCCGTCCGCCACCGGGCTcc                                                                                >  2:753861/1‑101 (MQ=255)
               ccccccGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTccccc                                                                             <  2:564029/98‑1 (MQ=25)
                cgcccGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTcccccc                                                                            <  2:5745835/99‑1 (MQ=25)
                cACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTcccccc                                                                            >  1:5246225/1‑101 (MQ=37)
                cACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTcccccc                                                                            >  2:4430323/1‑101 (MQ=37)
                cACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTcccccc                                                                            <  1:3665497/101‑1 (MQ=37)
                cACCCGCCCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTcccccc                                                                            <  2:2145222/101‑1 (MQ=25)
                  cccGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCgg                                                                          >  1:431406/1‑101 (MQ=25)
                   ccGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCggg                                                                         <  2:3248960/101‑1 (MQ=25)
                   ccGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCggg                                                                         >  1:3910379/1‑101 (MQ=25)
                    cgcgcgCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGa                                                                        <  2:1585252/101‑1 (MQ=25)
                     gcgcgcCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGat                                                                       <  2:566809/101‑1 (MQ=25)
                       gcgcCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGatat                                                                     <  1:753929/101‑1 (MQ=25)
                         gccgcGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATAttt                                                                   >  1:3597065/1‑101 (MQ=25)
                          ccgcGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTg                                                                  <  2:3028916/101‑1 (MQ=25)
                          ccgcGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTg                                                                  <  2:1728269/101‑1 (MQ=25)
                           cgcgTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgc                                                                 >  2:3152893/1‑101 (MQ=25)
                            gcgTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgca                                                                <  1:3074085/101‑1 (MQ=25)
                            gcgTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgca                                                                <  2:2392124/101‑1 (MQ=25)
                             ggTCTGGCCCAGCGCGGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcag                                                               <  2:4669248/100‑1 (MQ=21)
                              gTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GGGCTGCCCCCCGTCCGCAAGCGGGCCCCCCCCCGGGTATTTgcagc                                                              >  1:4743777/1‑101 (MQ=25)
                               tCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcagca                                                             <  1:1097514/101‑1 (MQ=25)
                               tCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcagca                                                             <  1:1295757/101‑1 (MQ=25)
                               tCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcagca                                                             <  2:4712387/101‑1 (MQ=25)
                                cTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCaa                                                            <  2:1850605/101‑1 (MQ=25)
                                 tGGCCCGGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAg                                                           <  1:4722034/101‑1 (MQ=21)
                                 tGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCCTCCGCAAGCGGACTCCCCCCGGGAAATTTGCAGGaaa                                                           >  2:3875878/1‑100 (MQ=25)
                                  ggCCCAGCGCCGCGGCGGGGGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGa                                                          <  1:5686735/101‑1 (MQ=21)
                                  ggCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GGGCGGGCCCCCCGCCGCCCGCGGACCCCCCCCCGGGGATTTTCCGGAAGa                                                          >  1:3747665/1‑101 (MQ=255)
                                   gCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGAt                                                         <  2:4055277/101‑1 (MQ=21)
                                    cccCGCGCCGCGGCGGGCGGGCGGCAGGACAGGCCGAGTAA.GGGGGGCGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATg                                                        <  1:2682086/101‑1 (MQ=18)
                                     ccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGa                                                       <  2:4612123/101‑1 (MQ=21)
                                     ccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGa                                                       >  1:1935271/1‑101 (MQ=21)
                                     ccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGa                                                       <  2:2800768/101‑1 (MQ=21)
                                        gcgcCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAc                                                    <  2:2526510/101‑1 (MQ=21)
                                         cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                                   <  1:2709425/101‑1 (MQ=21)
                                          gccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCg                                                  <  1:2046516/101‑1 (MQ=21)
                                           ccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCgg                                                 >  1:5308800/1‑101 (MQ=21)
                                             gcggcgGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.At                                              <  2:5533574/101‑1 (MQ=21)
                                              cggcggGCGCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATg                                             <  2:3747437/101‑1 (MQ=21)
                                               ggcgggcgCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGc                                            <  1:1229830/101‑1 (MQ=21)
                                                 cgggcgCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCt                                          >  1:3114725/1‑101 (MQ=21)
                                                 cgggcgCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCt                                          <  1:5606574/101‑1 (MQ=21)
                                                  gggcgCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCtt                                         <  2:5331538/101‑1 (MQ=21)
                                                  gggcgCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCtt                                         <  2:4272241/101‑1 (MQ=21)
                                                  gggcgCGCGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACCCCCCCCCGGACATTTGCAGCCAGATGAAAACGG.AGGccgt                                         >  1:5716720/1‑99 (MQ=18)
                                                    gcgcgcGGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTc                                       <  1:3257067/101‑1 (MQ=21)
                                                      gcgcgGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCAt                                     <  2:5376790/101‑1 (MQ=21)
                                                        gcgGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCt                                   <  2:3431272/101‑1 (MQ=21)
                                                         cgGCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCtt                                  <  2:4126816/101‑1 (MQ=21)
                                                          ggCACGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTg                                 <  1:4043740/101‑1 (MQ=21)
                                                            cgcGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGcc                               <  1:334478/99‑1 (MQ=18)
                                                             aCGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGccc                              <  1:513639/101‑1 (MQ=21)
                                                             aCGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGccc                              >  1:2371653/1‑101 (MQ=21)
                                                             aCGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGccc                              <  1:2862297/101‑1 (MQ=21)
                                                             aCGACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGGCCGCCAGCGGACTCCCCCCCGGATATTTGCCGCAAGAAGAAAACGGAAGGCCTTTCTTCtgccc                               >  1:1041465/1‑97 (MQ=9)
                                                               gACAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCCG.ATGCCTTTCATCTTGgccca                            >  1:2620391/1‑99 (MQ=18)
                                                                aCAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGaa                           <  2:5258250/101‑1 (MQ=21)
                                                                aCAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGaa                           <  1:1126193/101‑1 (MQ=21)
                                                                 cAGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGaaa                          <  2:3984427/101‑1 (MQ=21)
                                                                  aGGCCGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAat                         <  2:4629544/101‑1 (MQ=21)
                                                                  aGGCCGAGTAA.GGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAata                        <  1:4922381/101‑1 (MQ=255)
                                                                    gccgggtagggggggg.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAatat                       <  1:290385/91‑1 (MQ=18)
                                                                     ccGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATc                      <  2:1002550/101‑1 (MQ=18)
                                                                      cGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATcc                     >  1:1610409/1‑101 (MQ=18)
                                                                      cGAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATcc                     <  1:1768512/101‑1 (MQ=21)
                                                                       gAGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATccc                    >  2:5357442/1‑101 (MQ=18)
                                                                        ggaagggggggg.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCg                   <  2:3377528/95‑1 (MQ=18)
                                                                        aGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCCAGCGGACTCCCACCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCg                   <  2:2268504/101‑1 (MQ=18)
                                                                        aGTAAGGGGGGG.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCg                   <  2:1405575/101‑1 (MQ=18)
                                                                         gtaaggggggg.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCgg                  <  2:4477155/96‑1 (MQ=255)
                                                                         gtaaggggggg.gCGCGGCCCACCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCgg                  <  2:1457348/96‑1 (MQ=18)
                                                                           aaggggggg.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCgggg                >  2:190980/4‑101 (MQ=255)
                                                                            aggggggg.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCggggg               <  1:2298560/99‑1 (MQ=255)
                                                                            aggggggg.gCGCTGCCCCCCGTCCGCAAGCGGACCCCCCCCGGGGTTTTTGCCGCCAGCTGACACCGG.ATGCCCTTCATCTTCCCCTAAATAATCCggggg               >  2:2288534/3‑101 (MQ=12)
                                                                             ggggggg.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCgggggg              <  1:4207758/100‑1 (MQ=255)
                                                                               ggggg.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg            <  1:5728365/101‑1 (MQ=255)
                                                                                gggg.gCGCTGCCCCCCGGCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa           <  2:5200174/101‑1 (MQ=255)
                                                                                gggc.gCGCTCCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa           <  2:5703803/101‑1 (MQ=18)
                                                                                 ggg.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGaa          >  2:2477211/1‑101 (MQ=255)
                                                                                  gg.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAAt         <  1:1688307/101‑1 (MQ=255)
                                                                                  gg.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAAt         <  1:679154/101‑1 (MQ=255)
                                                                                  gg.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAAt         <  1:1563632/101‑1 (MQ=255)
                                                                                     gcgcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAATgg       >  1:1959831/1‑101 (MQ=255)
                                                                                        cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTTCCGCAAGATGAAACCGG.ATGCCTTTCATCTGGCCCGAAATATCCCGGGGGGTGAATGGcaa    >  2:1679067/1‑99 (MQ=21)
                                                                                         tGCCCCCCGTCCGCAAGCGGGCTCCCCCCGGGATATTTGCAGCAAGATGCAACCGG.GTGCCTTTCCTCTTTTCCCAAATAACCCGGGGGGGGacaggggaa   >  2:4961732/1‑93 (MQ=12)
                                                                                          gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGGCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAATGGCGaac  >  1:5799507/1‑99 (MQ=18)
                                                                                          gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGaag  <  2:4908975/101‑3 (MQ=255)
                                                                                          gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGaag  <  1:4336343/101‑3 (MQ=255)
                                                                                                 |                                                                                              
TGATGGCCGAGCTGCCCACCCGCGCGCCGCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA.GGGGGG.GCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGG.ATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGAGC  >  NC_014034/2493876‑2494064

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