Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,003 A→G 5.9% intergenic (+51/‑191) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0030A→G 5.9% 15.4 ‑2.6 ‑0.3152intergenic (+51/‑191)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (6/3):  ref base (50/93):  total (56/96)
Fisher's exact test strand distribution p-value = 7.61e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.22e-01

Only 100 of 169 total aligned reads displayed.

GCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA.GGGGG.G.GCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AGCAGA.AGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTG  >  NC_014034/2493904‑2494102
                                                                                                      |                                                                                                         
gcgTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCTCGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgca                                                                                                           <  1:4044647/101‑1 (MQ=21)
   tCTGGCCCGGCGCCGCGGCGGTCGCGCGGCTCGACAGCCGGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcagca                                                                                                        <  1:4519000/101‑1 (MQ=18)
   tCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTgcagca                                                                                                        <  2:1373026/101‑1 (MQ=25)
     tGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAg                                                                                                      <  2:3276650/101‑1 (MQ=21)
     tGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAg                                                                                                      <  1:4623226/101‑1 (MQ=21)
      ggCCCAGCGCCGCTGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGa                                                                                                     <  2:4420555/101‑1 (MQ=21)
      ggCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGa                                                                                                     <  1:425134/101‑1 (MQ=21)
      ggCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCCGGGTATTTGCAGCCAGa                                                                                                     >  2:3449973/1‑101 (MQ=21)
       gCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGAt                                                                                                    <  1:2153340/101‑1 (MQ=21)
        cccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATg                                                                                                   <  1:3897437/101‑1 (MQ=21)
         ccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCTCCCGGGATATTTACAGCAAGATGa                                                                                                  <  2:3203720/101‑1 (MQ=21)
         ccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGa                                                                                                  <  1:1854239/101‑1 (MQ=21)
         ccAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGa                                                                                                  <  2:1260855/101‑1 (MQ=21)
          cAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaa                                                                                                 <  2:86021/101‑1 (MQ=21)
           aGCGCGGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaaa                                                                                                <  2:1850840/101‑1 (MQ=18)
           aGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATNTTTGCAGGaannnnnnn                                                                                                >  2:453059/1‑94 (MQ=18)
           aGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaaa                                                                                                <  1:3263804/101‑1 (MQ=21)
           aGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaaa                                                                                                <  2:1071748/101‑1 (MQ=21)
             ctccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                                                                              <  2:720224/99‑1 (MQ=18)
             cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGTATTTGCAGCAAGATGAAAcc                                                                                              >  2:3893179/1‑101 (MQ=21)
                cgcgGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGa                                                                                           <  1:2415967/101‑1 (MQ=21)
                cgcgGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGa                                                                                           <  2:3085716/101‑1 (MQ=21)
                cgcgGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGa                                                                                           <  1:3827145/101‑1 (MQ=21)
                cgcgCCGGGCGCGCGCCCCGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGa                                                                                           <  2:4138218/101‑1 (MQ=17)
                 gcggcgGGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAt                                                                                          <  2:3799038/101‑1 (MQ=21)
                  cggcggGCGCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATg                                                                                         <  1:2093370/101‑1 (MQ=21)
                   ggcgggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCACCGGGAGATTTGCAGCAAGATGAAACAGGCTGc                                                                                        >  1:4013220/1‑101 (MQ=18)
                    gcgggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGcc                                                                                       <  2:3405661/101‑1 (MQ=21)
                     cgggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCt                                                                                      <  1:576105/101‑1 (MQ=21)
                      gggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                     <  2:3691524/101‑1 (MQ=21)
                       ggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                                                                    <  2:1947490/101‑1 (MQ=21)
                       ggcgCGCGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                                                                    <  2:399985/101‑1 (MQ=21)
                        gcgcgcGGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                                                                   <  2:363510/101‑1 (MQ=21)
                            gcgGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCt                                                                               <  1:2334011/101‑1 (MQ=21)
                             cgGCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                              <  1:1966740/101‑1 (MQ=21)
                              ggCACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTg                                                                             >  2:2736730/1‑101 (MQ=21)
                                cACGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGcc                                                                           >  2:1311081/1‑101 (MQ=21)
                                 aCGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATTTTGccc                                                                          >  2:3578950/1‑101 (MQ=18)
                                  cGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                         <  1:1254986/101‑1 (MQ=21)
                                  cGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                         <  2:4362146/101‑1 (MQ=21)
                                  cGACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                         <  1:2796169/101‑1 (MQ=21)
                                   gACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGa                                                                        <  2:110926/101‑1 (MQ=21)
                                   gACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGa                                                                        <  2:2833658/101‑1 (MQ=21)
                                   gACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGa                                                                        <  1:3525447/101‑1 (MQ=21)
                                   gACAGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGACGCCTTTCATCTTTGCCGa                                                                        >  2:1674435/1‑101 (MQ=17)
                                      aGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat                                                                     <  1:3816435/101‑1 (MQ=21)
                                      aGGCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat                                                                     <  2:647773/101‑1 (MQ=21)
                                       ggCCGAGTAA.GGGGG.GCGGTCGGCCCCGCGCCCGCAAGGGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCCTCTTGCCCGAAata                                                                    <  2:2306986/101‑1 (MQ=17)
                                        gCCGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat                                                                   <  2:1324367/101‑1 (MQ=21)
                                         tcGAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc                                                                  <  1:4054669/100‑1 (MQ=18)
                                         ccGAGTAGGGGGGG.G.GCGCTCCCCCCCGTCCGCAAGCGGACCCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc                                                                  <  2:2063019/101‑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                                                                <  2:3429093/101‑1 (MQ=18)
                                           gAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                <  1:4192533/101‑1 (MQ=18)
                                           gAGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATccc                                                                >  1:3418496/1‑101 (MQ=18)
                                            aGTAAGGGGGG.G.GCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCg                                                               >  1:1000152/1‑101 (MQ=18)
                                             gtaagggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGGTATTTGCAGCAAGATGAAAACGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                              >  1:906249/6‑101 (MQ=17)
                                             gtaagggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                              >  1:1486152/6‑101 (MQ=255)
                                             gtaagggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                              >  1:1504263/6‑101 (MQ=255)
                                             gtaagggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                              <  1:1350477/96‑1 (MQ=255)
                                              taagggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCCGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggg                                                             >  1:3180596/5‑101 (MQ=18)
                                               aagggggg.g.gGCCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                            <  2:1000069/98‑1 (MQ=18)
                                               aagggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                            >  1:1628290/4‑101 (MQ=255)
                                                 gggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggg                                                          >  1:2629620/2‑101 (MQ=255)
                                                  ggggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGt                                                         <  2:960816/101‑1 (MQ=255)
                                                   gggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                        <  1:2172571/101‑1 (MQ=255)
                                                   gggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                        <  1:4511415/101‑1 (MQ=255)
                                                   gggg.g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                        <  1:1668650/101‑1 (MQ=255)
                                                        g.gCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                   <  1:4340028/101‑1 (MQ=255)
                                                            gcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..g                                              >  2:2159991/1‑101 (MQ=255)
                                                             cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGG..G.a                                            >  1:4349045/1‑101 (MQ=18)
                                                             cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGAAGCAAGATGAAACCGGATGCCTTTCATTTTGACCGAAATATCCCGGGGGGTGAA.TGGC..G.a                                            >  2:2703254/1‑101 (MQ=18)
                                                               gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aag                                          <  1:704842/101‑3 (MQ=255)
                                                               gCCCCCCGTCCGCAAGCGGAATCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.GGGC..G.aac                                          >  2:2183415/1‑99 (MQ=18)
                                                                     cgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagag                                    <  1:2418397/99‑9 (MQ=255)
                                                                      gtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagagg                                   >  2:4129599/2‑92 (MQ=255)
                                                                            aaGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGCa                             <  1:814009/101‑1 (MQ=18)
                                                                              gCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGGCTTTCATCTTGCCCGAAATGTCAACGGGGGGGAA.TGGC..GAAGCCAG.AGGGGGGCagc                           >  1:3989609/1‑99 (MQ=17)
                                                                               cGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGCagcg                          >  1:2943247/1‑98 (MQ=18)
                                                                                  aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGG..GAAGCCAG.AGGGGGGC.AGCG.ccc                     >  2:4255159/1‑101 (MQ=17)
                                                                                    tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCG.ccccc                   >  1:2173295/1‑101 (MQ=18)
                                                                                      cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCG.cccccct                 <  2:2243612/101‑3 (MQ=17)
                                                                                       ccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCG.cccccctg                <  2:2984987/101‑4 (MQ=17)
                                                                                          cGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AAGACAGAGGGGGGC.AGCG.cccccctgacc             >  1:1814513/1‑95 (MQ=17)
                                                                                            ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCA.AGAGGGGGCAGCG.CCCCCTGCCCcg            >  1:1282178/1‑101 (MQ=2)
                                                                                            ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCA.AGAGGGGGCAGCG.CCCCCTGCCCcg            >  2:1903843/1‑101 (MQ=0)
                                                                                             gATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgc           >  2:4290240/1‑100 (MQ=2)
                                                                                              atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgcc          >  1:3869648/1‑99 (MQ=0)
                                                                                              atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgcc          <  2:2444885/101‑3 (MQ=2)
                                                                                              atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgcc          <  1:1127028/101‑3 (MQ=2)
                                                                                              atatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGCGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCCCCCTGCCCcggc          >  2:477874/1‑101 (MQ=17)
                                                                                               tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGT..GAAGCCAG.AGGGGGGC.AGCGCCCCCCTGCCCcggcg         >  1:1457259/1‑101 (MQ=17)
                                                                                               tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCCCCCTGACCcggcg         >  2:1611164/1‑101 (MQ=17)
                                                                                                atTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCA.AGAGGGGGCAGCG.CCCCCTGCCCcgccgg        >  1:4479250/1‑101 (MQ=2)
                                                                                                  ttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGCAGCGG.CCCCCGGCCCCGGCGGTc      >  2:714330/1‑101 (MQ=18)
                                                                                                   tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTcc     <  2:1703597/101‑1 (MQ=17)
                                                                                                     cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGt   >  2:4340374/1‑101 (MQ=18)
                                                                                                     cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAG.AGGGGGGC.AGCG.CCCCCCGCCCCCGCGGGcccgg  >  1:4092509/1‑98 (MQ=17)
                                                                                                      |                                                                                                         
GCGTCTGGCCCAGCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA.GGGGG.G.GCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AGCAGA.AGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTG  >  NC_014034/2493904‑2494102

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG