Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,014 A→G 4.2% intergenic (+62/‑180) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0140A→G 4.2% 18.6 ‑0.6 ‑1.0288intergenic (+62/‑180)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (6/6):  ref base (114/162):  total (120/168)
Fisher's exact test strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.07e-01

Only 100 of 318 total aligned reads displayed.

AGCGCCGCGGCGGGCGCGCGGCACGACAGGCC.GAGTAA..GGGG.GGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AG.CAGAAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTG  >  NC_014034/2493915‑2494113
                                                                                                       |                                                                                                         
aGCGCCGCGGCGGGCGCGCGGCACGACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCCAGATGaaa                                                                                                           >  1:4687584/1‑101 (MQ=21)
aGCGCCGCGGCGGGCGCGCGGCACGACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGaaa                                                                                                           <  1:4345290/101‑1 (MQ=21)
 gcgcCGCGGCGGGCGCGCGGCACGACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAc                                                                                                          <  2:5806286/101‑1 (MQ=21)
    ccgcGGCGGGCGCGCGGCACGACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCgg                                                                                                       <  1:1632771/101‑1 (MQ=21)
     cgggcggggcgCGCGGC.CGGCCGGCC.GGGTAAGGGGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGa                                                                                                      <  2:4891201/95‑1 (MQ=14)
      gcggcgGGCGCGCGGCACGACAGGCCGGGGAAA..GGGG.GGGGGCTGCCCCCCGTCCGGCAGCGGGCTCCCCCCCGGCAATTTGCCGCACAAAGAAACCGGaa                                                                                                     >  1:1170684/1‑100 (MQ=255)
         gcgggcgCGCGGCACGACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGcc                                                                                                  <  2:1668214/101‑1 (MQ=21)
         gccggcgCGCGGCACGACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGcc                                                                                                  <  1:6940405/98‑1 (MQ=18)
             gcgcgcGGCACGACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                                                                              <  2:1128414/101‑1 (MQ=21)
                 gcgGCACGACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCt                                                                                          <  1:458653/101‑1 (MQ=21)
                  cgGCAGGACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                         <  1:1844552/101‑1 (MQ=18)
                  cgGCACGACAGGCC.GAGTAAG.GGGG.GGGCGGTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAACCGATGCATTTCAttct                                                                                         >  2:3292121/1‑98 (MQ=17)
                  cgGCACGACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                         <  1:4932642/101‑1 (MQ=21)
                  cgGCACGACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                         <  1:576379/101‑1 (MQ=21)
                   ggCACGACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTg                                                                                        <  1:5481342/101‑1 (MQ=21)
                     cACGACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGcc                                                                                      <  1:884545/101‑1 (MQ=21)
                        gACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGGCCGa                                                                                   >  1:3212653/1‑101 (MQ=18)
                        gACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGa                                                                                   >  1:967734/1‑101 (MQ=21)
                        gACAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGa                                                                                   >  2:6135284/1‑101 (MQ=21)
                          gaGGTC.GAGTAAG.GGGG.GGGCGCCGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                 <  2:4485748/100‑1 (MQ=18)
                          cAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                 <  1:3644490/101‑1 (MQ=21)
                          cAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                 <  1:5408148/101‑1 (MQ=21)
                          cAGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                 <  2:6934405/101‑1 (MQ=21)
                           aGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat                                                                                >  1:530803/1‑101 (MQ=21)
                           aGGCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat                                                                                <  1:2055620/101‑1 (MQ=21)
                            ggCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                               >  2:6480146/1‑101 (MQ=21)
                            ggCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                               <  2:2710869/101‑1 (MQ=21)
                            ggCC.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                               >  1:6581394/1‑101 (MQ=21)
                              cc.GATTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc                                                                             <  2:40144/101‑1 (MQ=18)
                              cc.GAGTAA..GGGGTGAGCGCTGCCCCCCGTCCGCAAGCGTACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc                                                                             <  2:2836554/101‑1 (MQ=21)
                               c.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                            <  1:7103377/101‑1 (MQ=21)
                               c.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                            <  1:4164654/101‑1 (MQ=21)
                               c.GAGTAAG.GGGG.GGGCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                            <  2:1032815/101‑1 (MQ=18)
                                 cgactaaggggg.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatatac                                                                            >  1:120244/9‑99 (MQ=18)
                                      aaggggg.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                       >  2:707962/4‑101 (MQ=255)
                                      aaggggg.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                       <  1:2951550/98‑1 (MQ=255)
                                      aagcggg.gggCGCTCCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                       <  2:5713439/97‑1 (MQ=18)
                                       aggggg.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggg                                                                      >  2:202817/3‑101 (MQ=255)
                                       aggggg.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggg                                                                      >  1:1594922/3‑101 (MQ=255)
                                       aggggg.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATACCggggg                                                                      >  1:1715591/3‑101 (MQ=255)
                                        ggggg.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggg                                                                     >  2:1828193/2‑101 (MQ=255)
                                        ggggg.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggg                                                                     >  1:3122284/2‑101 (MQ=255)
                                        ggggg.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATACCCCgggggg                                                                     <  1:6828304/100‑1 (MQ=255)
                                         gggg.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggggg                                                                    >  2:4799099/1‑100 (MQ=255)
                                          ggg.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggggg                                                                   >  1:3456816/1‑99 (MQ=255)
                                          ggg.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                                   <  2:2973957/101‑1 (MQ=255)
                                          ggg.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                                   <  1:1978697/101‑1 (MQ=255)
                                            g.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGaa                                                                 <  1:6962873/101‑1 (MQ=255)
                                            g.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGaa                                                                 >  2:4137507/1‑101 (MQ=255)
                                            g.gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGAAGCAAGATGACACCGGATGCCTTTCATCTTGCCCGAAATATACCGGGGGGTGaa                                                                 >  2:4653910/1‑101 (MQ=17)
                                              gggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.t                                                               <  1:6616493/101‑1 (MQ=255)
                                                  gcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..g                                                         >  2:5334729/1‑101 (MQ=255)
                                                  gcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..g                                                         <  2:2157931/101‑1 (MQ=255)
                                                   cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.a                                                       <  2:4569864/101‑1 (MQ=255)
                                                    tccccccCGTCCGAAAGCGGCCT.CCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aag                                                     <  1:4904669/99‑3 (MQ=18)
                                                     gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aag                                                     >  1:4097084/1‑99 (MQ=255)
                                                       ccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcca                                                  <  1:6196766/101‑6 (MQ=255)
                                                       ccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGACCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcca                                                  >  2:6712796/1‑96 (MQ=21)
                                                       ccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGAAAGATGAAACCGGATGCCTTTCATCTTGCCCGAGATATCCCGGGGGGTGAA.TGGC..ggacgca                                                   >  2:1215755/1‑95 (MQ=18)
                                                         cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccag                                                 <  1:2613364/97‑7 (MQ=255)
                                                                    gCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCagc                                      <  1:3524300/101‑3 (MQ=18)
                                                                    gCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCagc                                      >  2:5420299/1‑99 (MQ=18)
                                                                    gCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCagc                                      >  1:5933131/1‑99 (MQ=18)
                                                                     cGGCCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCagcg                                     <  2:4391520/101‑4 (MQ=18)
                                                                     cGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGCagcg                                     >  1:381771/1‑98 (MQ=18)
                                                                        aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.ccc                                >  2:6585192/1‑101 (MQ=18)
                                                                         cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccc                               <  2:718542/101‑1 (MQ=17)
                                                                           ccccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccccc                             >  1:2407036/1‑100 (MQ=17)
                                                                           ccccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccccc                             <  2:1570929/101‑2 (MQ=17)
                                                                           ccccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGG.AGGCG.cccccc                             >  2:2902251/1‑100 (MQ=17)
                                                                            cccccGTGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccccct                            <  1:3292237/101‑3 (MQ=17)
                                                                            cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccccct                            <  2:3550655/101‑3 (MQ=17)
                                                                            cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccccct                            <  2:3617111/101‑3 (MQ=17)
                                                                             ccccGGGATATTGGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGCGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCG.cccccctg                           >  1:6308386/1‑98 (MQ=17)
                                                                                 gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGcccc                        <  2:4096890/101‑1 (MQ=17)
                                                                                 gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGcccc                        <  2:1733916/101‑1 (MQ=17)
                                                                                  ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcg                       <  2:6340765/101‑1 (MQ=0)
                                                                                  ggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcg                       >  2:3475203/1‑101 (MQ=18)
                                                                                  ggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcg                       <  2:5653150/101‑1 (MQ=17)
                                                                                  ggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcg                       >  1:5969404/1‑101 (MQ=18)
                                                                                   gATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgc                      >  1:6081322/1‑100 (MQ=2)
                                                                                   gATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcgg                      >  2:31145/1‑101 (MQ=18)
                                                                                     tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcggcg                    >  2:2512571/1‑101 (MQ=17)
                                                                                     tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcgccg                    >  2:3164905/1‑101 (MQ=17)
                                                                                      atTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgccgg                   >  1:3954242/1‑101 (MQ=2)
                                                                                       tttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcggcggg                  >  1:3783275/1‑100 (MQ=17)
                                                                                       tttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGt                  <  1:6135596/101‑1 (MQ=17)
                                                                                       tttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGt                  <  2:6449675/101‑1 (MQ=17)
                                                                                        ttGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AG.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgccggat                 <  2:5078699/101‑3 (MQ=2)
                                                                                         tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATACCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGACGGGC..CCCCCGGCCCCGGCGGTcc                >  2:4252375/1‑101 (MQ=17)
                                                                                           cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AA.GCCAGGGGGGGGCAGCGCCCCCCTGCCCCGGCGGTCCGt              <  1:1021714/101‑1 (MQ=18)
                                                                                               aaGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGCCCGTggcg          >  1:3383889/1‑101 (MQ=17)
                                                                                                aGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGGGGTCCGGgccgg         >  1:2858290/1‑101 (MQ=9)
                                                                                                 tggaaaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggcggc        >  2:1944999/5‑101 (MQ=17)
                                                                                                   tGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCAGTGGCGgcgc      >  2:605370/1‑101 (MQ=17)
                                                                                                    gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGcc     >  1:334723/1‑101 (MQ=18)
                                                                                                      aaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCtt   >  1:2830502/1‑101 (MQ=18)
                                                                                                       aCCGGATGCCTTTCATCTTGCCCGAAATATGCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg  <  1:229743/101‑1 (MQ=17)
                                                                                                       aCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg  <  2:375861/101‑1 (MQ=17)
                                                                                                       aCCGGATGCCTTTCATCTTGCCCGAAATAGCCCGGGGGGTGAA.TGGC..GAAG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg  <  2:1348192/101‑1 (MQ=17)
                                                                                                       |                                                                                                         
AGCGCCGCGGCGGGCGCGCGGCACGACAGGCC.GAGTAA..GGGG.GGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AG.CAGAAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTG  >  NC_014034/2493915‑2494113

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG