Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,023 C→A 6.8% intergenic (+71/‑171) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0230C→A 6.8% 15.7 ‑1.2 ‑0.2133intergenic (+71/‑171)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (6/3):  ref base (55/69):  total (61/72)
Fisher's exact test strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.90e-01

Only 100 of 145 total aligned reads displayed.

GCGGGCGCGCGGCACGACAGGCCGAGTAA.GGGGG.G.GCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AG.CAGAAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGG  >  NC_014034/2493924‑2494122
                                                                                                      |                                                                                                         
 cgggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCt                                                                                                          <  1:576105/101‑1 (MQ=21)
  gggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                                         <  1:4577252/101‑1 (MQ=21)
  gggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                                         <  2:3907789/101‑1 (MQ=21)
  gggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                                         <  2:3691524/101‑1 (MQ=21)
   ggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCGCGTCCGCAAGCGGACTCCCCCCGGGGTATTTGCAGCAAGATGAAACCGGATGCCttt                                                                                                        <  1:1624367/101‑1 (MQ=21)
   ggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                                                                                        <  2:399985/101‑1 (MQ=21)
    gcgcgcGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                                                                                       <  2:363510/101‑1 (MQ=21)
    gcgcgcGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                                                                                       <  1:4452450/101‑1 (MQ=21)
         cgGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                                  <  1:2097046/101‑1 (MQ=21)
         cgGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                                  <  2:1020681/101‑1 (MQ=21)
         cgGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                                  <  1:1966740/101‑1 (MQ=21)
          ggCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTg                                                                                                 >  2:1364262/1‑101 (MQ=21)
            cACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGcc                                                                                               >  2:1311081/1‑101 (MQ=21)
              cGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                             <  1:2581620/101‑1 (MQ=21)
              cGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                             <  1:2796169/101‑1 (MQ=21)
              cGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                             <  2:4362146/101‑1 (MQ=21)
              cGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                             <  1:1254986/101‑1 (MQ=21)
               gACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGa                                                                                            <  2:110926/101‑1 (MQ=21)
               gACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGa                                                                                            <  2:2833658/101‑1 (MQ=21)
               gACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGACGCCTTTCATCTTTGCCGa                                                                                            >  2:1674435/1‑101 (MQ=17)
                aCAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaa                                                                                           <  2:4427230/101‑1 (MQ=21)
                 cAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                          <  2:1981220/101‑1 (MQ=21)
                  aGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat                                                                                         <  1:3816435/101‑1 (MQ=21)
                  aGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat                                                                                         <  2:647773/101‑1 (MQ=21)
                   ggCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                                        <  2:3559894/101‑1 (MQ=21)
                   ggCCGAGTAA.GGGGG.GCGGTCGGCCCCGCGCCCGCAAGGGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCCTCTTGCCCGAAata                                                                                        <  2:2306986/101‑1 (MQ=17)
                    gCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat                                                                                       <  1:3778516/101‑1 (MQ=21)
                     tcGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc                                                                                      <  1:4054669/100‑1 (MQ=18)
                      cGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                                     <  1:1686562/101‑1 (MQ=18)
                      cGAGTAAGGGGGG.G.GCCCTGCCCCCCGTCCGCAAGCCGACTCCCCCCGGCATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                                     <  2:3151033/101‑1 (MQ=18)
                      agAGAAA.GGGGGAG.TCGCTGCCCCCCGTCCACAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                                     <  1:530034/100‑1 (MQ=18)
                       gAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                                    <  1:4192533/101‑1 (MQ=18)
                       gAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                                    >  1:3418496/1‑101 (MQ=18)
                       gAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                                    <  2:3429093/101‑1 (MQ=18)
                        aGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCg                                                                                   <  1:1833890/101‑1 (MQ=18)
                        aGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCg                                                                                   >  1:1000152/1‑101 (MQ=18)
                        aGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAAGCATTTAATCTTGCCAGAAATATTCCg                                                                                   >  2:187457/1‑101 (MQ=17)
                         gtaagggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGTATTTGCAGCAAGATGAAAACGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                                                  >  1:906249/6‑101 (MQ=17)
                         gtaagggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                                                  <  1:1350477/96‑1 (MQ=255)
                         gtaagggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                                                  >  1:1486152/6‑101 (MQ=255)
                           aagggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                                >  1:1628290/4‑101 (MQ=255)
                              ggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggggg                                                                             >  2:3028533/1‑100 (MQ=255)
                              ggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGt                                                                             <  2:960816/101‑1 (MQ=255)
                               gggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggggg                                                                            >  2:1621830/1‑99 (MQ=255)
                               gggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                                            <  1:4511415/101‑1 (MQ=255)
                                ggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                                           <  2:2208328/101‑1 (MQ=255)
                                    g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                                       <  1:4340028/101‑1 (MQ=255)
                                    g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCAGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                                       <  1:3232344/101‑1 (MQ=255)
                                        gcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..g                                                                  >  2:2159991/1‑101 (MQ=255)
                                         cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGAAGCAAGATGAAACCGGATGCCTTTCATTTTGACCGAAATATCCCGGGGGGTGAA.TGGC..G.a                                                                >  2:2703254/1‑101 (MQ=18)
                                          tGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aa                                                               >  1:14024/1‑100 (MQ=255)
                                           gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aag                                                              <  1:704842/101‑3 (MQ=255)
                                           gCCCCCCGTCCGCAAGCGGAATCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.GGGC..G.aac                                                              >  2:2183415/1‑99 (MQ=18)
                                            cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcc                                                            <  1:434851/101‑5 (MQ=255)
                                            cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcc                                                            <  2:1568902/101‑5 (MQ=255)
                                               cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..G.aagccag                                                          >  1:2985216/5‑95 (MQ=18)
                                                ccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccaga                                                         <  2:432097/98‑8 (MQ=255)
                                                 cgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagag                                                        <  1:2418397/99‑9 (MQ=255)
                                                  gtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagagg                                                       >  2:4129599/2‑92 (MQ=255)
                                                        aaGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCa                                                 <  1:814009/101‑1 (MQ=18)
                                                         aGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCag                                                <  2:1264520/101‑2 (MQ=18)
                                                              aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGG..GAAG.CCAGAGGGGGGC.AGCG.ccc                                         >  2:4255159/1‑101 (MQ=17)
                                                              aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.ccc                                         >  2:2256315/1‑101 (MQ=18)
                                                               cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccc                                        >  2:486759/1‑101 (MQ=18)
                                                                tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCG.ccccc                                       <  2:1714876/101‑1 (MQ=17)
                                                                tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.ccccc                                       >  1:2173295/1‑101 (MQ=18)
                                                                  cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccccct                                     <  1:2160141/101‑3 (MQ=17)
                                                                   ccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccccctg                                    <  2:2984987/101‑4 (MQ=17)
                                                                   ccccGGGATATCTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccccctg                                    >  1:2729047/1‑98 (MQ=17)
                                                                       gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCAGCAC.CCCCCTGcccc                                 >  2:344576/1‑101 (MQ=17)
                                                                        ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcg                                >  2:4133702/1‑101 (MQ=2)
                                                                        ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcg                                >  2:1903843/1‑101 (MQ=0)
                                                                        ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcg                                <  1:2338829/101‑1 (MQ=0)
                                                                        ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcg                                >  1:1282178/1‑101 (MQ=2)
                                                                        ggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcg                                >  1:991533/1‑101 (MQ=18)
                                                                        ggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcg                                <  1:1726617/101‑1 (MQ=17)
                                                                         gATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgc                               >  2:4290240/1‑100 (MQ=2)
                                                                          atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgcc                              <  1:1127028/101‑3 (MQ=2)
                                                                          atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgcc                              >  1:3869648/1‑99 (MQ=0)
                                                                          atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgcc                              <  2:2444885/101‑3 (MQ=2)
                                                                           tatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgccg                             <  2:2091493/101‑1 (MQ=2)
                                                                            atTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgccgg                            >  1:4479250/1‑101 (MQ=2)
                                                                              ttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCAGCGG.CCCCCGGCCCCGGCGGTc                          >  2:714330/1‑101 (MQ=18)
                                                                              ttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTc                          <  1:1467624/101‑1 (MQ=17)
                                                                               tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTcc                         >  1:3926664/1‑101 (MQ=18)
                                                                                gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCg                        >  1:1658069/1‑101 (MQ=17)
                                                                                gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCACAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCg                        >  1:4037011/1‑101 (MQ=18)
                                                                                 cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGCGAA.TGGC..GAAG.CCACAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGt                       >  1:2256617/1‑101 (MQ=18)
                                                                                   gcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAG.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTgg                     >  2:294835/1‑101 (MQ=18)
                                                                                    caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggc                    >  2:2049065/1‑101 (MQ=18)
                                                                                    caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGAGAA.TGGC..GAAG.CCAGAGGGGGGCCGACC.CCCCCTGCCCCGGCGGTCCGCggg                    >  1:677081/1‑100 (MQ=17)
                                                                                         tGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGgcgc               >  1:4604971/1‑101 (MQ=17)
                                                                                           aaaCCGGATGCCTTTCATCTTCCCCGAAATGTCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCt             <  2:3182635/101‑1 (MQ=17)
                                                                                             ccgggctGCCTTTAATCTTG.CCAAAACTTGCCGGGGGGGGAAGGGGC..GAAG.CGGGGGGGGGCCACCCC.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg           <  2:2728849/95‑1 (MQ=255)
                                                                                                ggATGCCTTTCATCTTTCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcg        <  1:3085165/101‑1 (MQ=17)
                                                                                                 gATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgg       >  1:1199947/1‑101 (MQ=17)
                                                                                                   tGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggc     >  2:3668466/1‑101 (MQ=18)
                                                                                                    gCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcg    >  1:3532805/1‑101 (MQ=18)
                                                                                                      cTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcggg  <  2:3699520/101‑1 (MQ=17)
                                                                                                      cTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcggg  >  1:4153737/1‑101 (MQ=18)
                                                                                                      |                                                                                                         
GCGGGCGCGCGGCACGACAGGCCGAGTAA.GGGGG.G.GCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AG.CAGAAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGG  >  NC_014034/2493924‑2494122

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG