Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,061 G→A 9.4% intergenic (+109/‑133) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0610G→A 9.4% 26.4 ‑2.0 ‑0.3128intergenic (+109/‑133)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/8):  ref base (66/50):  total (70/58)
Fisher's exact test strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.47e-01

Only 100 of 135 total aligned reads displayed.

GCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.A.G.CAGA.AGGGGGGCAAGCG.CC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGA  >  NC_014034/2493961‑2494152
                                                                                                       |                                                                                                 
gcTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..g                                                                                                   <  1:3097965/101‑1 (MQ=255)
 cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.a                                                                                                 <  2:1195484/101‑1 (MQ=255)
  tGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aa                                                                                                <  1:3568460/101‑2 (MQ=255)
  tGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aa                                                                                                <  2:2495794/101‑2 (MQ=255)
   gCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aag                                                                                               <  1:1081080/101‑3 (MQ=255)
    cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCCGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGAGAA.TGGC..G.aagcc                                                                                             >  1:2946572/1‑97 (MQ=18)
    cccccCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagc                                                                                              <  1:2808866/101‑4 (MQ=255)
    cccccCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagc                                                                                              <  2:3332561/101‑4 (MQ=255)
    cccccCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..G.aagc                                                                                              >  2:3866505/1‑98 (MQ=18)
     ccccG.TCCTCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATTCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcca                                                                                            <  1:1331605/101‑6 (MQ=18)
     ccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcca                                                                                            <  1:1897539/101‑6 (MQ=255)
       cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccag                                                                                           <  1:3438156/97‑7 (MQ=255)
       cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccag                                                                                           >  1:3570867/5‑95 (MQ=255)
       cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccag                                                                                           <  1:769318/97‑7 (MQ=255)
       cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccag                                                                                           <  2:1226321/97‑7 (MQ=255)
         cgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagag                                                                                         >  1:3401873/3‑93 (MQ=255)
          gtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..G.aagccagagg                                                                                        >  1:1827773/2‑92 (MQ=21)
           tCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagaggg                                                                                       >  2:2511595/1‑91 (MQ=255)
           tCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagaggg                                                                                       <  1:1724067/101‑11 (MQ=255)
                aaGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGCa                                                                                >  2:2302328/1‑101 (MQ=18)
                aaGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAA.G.CCAG.AGGGGGGCa                                                                                >  1:4364755/1‑101 (MQ=18)
                 aGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGCag                                                                               <  1:2028967/101‑2 (MQ=18)
                 aGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGCag                                                                               <  2:413502/101‑2 (MQ=18)
                  gCGGACTCCCCCCGGGATTTTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGCagc                                                                              <  2:4145191/101‑3 (MQ=18)
                  gCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGCagc                                                                              >  2:450441/1‑99 (MQ=18)
                    ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAA.G.CCAG.AGGGGGGCagcgc                                                                            >  1:815323/1‑97 (MQ=17)
                       cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCG.cc.cc                                                                      >  2:688702/1‑101 (MQ=17)
                        tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGCAATATCCCGGGGGGGGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCG.cc.ccc                                                                     >  1:2395499/1‑101 (MQ=17)
                        tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.GCCAG..AGGGGGGC.AGCG.cc.ccc                                                                     <  2:4142574/101‑1 (MQ=17)
                         ccccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCG.cc.cccc                                                                    <  2:1952332/101‑2 (MQ=17)
                          cccccGGGATATTTGCAGCAAGATGAAACCGGATCCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCG.cc.cccct                                                                   <  1:19700/101‑3 (MQ=17)
                           ccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAT.TGGC..GAA.G.CCGG.GGGGGGGCGGGCC.cc.ccccc                                                                   >  2:3136668/1‑99 (MQ=17)
                           ccccGGGATATGTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCG.cc.cccctg                                                                  <  2:1322038/101‑4 (MQ=17)
                             ccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCG.cc.cccctgcc                                                                <  2:2951645/101‑6 (MQ=17)
                              cGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGACC.CCCTGccc                                                                >  2:1214419/1‑101 (MQ=17)
                              cGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATACCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGccc                                                                >  1:2563268/1‑101 (MQ=17)
                               gggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCA.AGAGGGGGCAGCG.CC.CCCTGcccc                                                               <  1:2186554/101‑1 (MQ=0)
                               gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGcccc                                                               <  1:3208022/101‑1 (MQ=17)
                               gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGcccc                                                               >  2:1325244/1‑101 (MQ=18)
                               gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGcccc                                                               >  1:576951/1‑101 (MQ=17)
                                ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcg                                                              >  2:1447132/1‑101 (MQ=0)
                                ggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGAGAA.GGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCcg                                                              >  1:4009471/1‑101 (MQ=17)
                                ggATANNTGCAGAAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCcg                                                              <  1:1594892/101‑1 (MQ=17)
                                 gATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcgc                                                             <  1:2793530/101‑2 (MQ=0)
                                 gATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCcgg                                                             >  1:2613171/1‑101 (MQ=17)
                                  gtatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcgcc                                                            <  2:1429237/100‑3 (MQ=2)
                                  atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcgcc                                                            <  1:1155047/101‑3 (MQ=2)
                                   tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCcggcg                                                           >  2:3118786/1‑101 (MQ=18)
                                   tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGCGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCACcggcg                                                           >  2:341801/1‑101 (MQ=17)
                                    atTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.A.A.GACAGAGGGGGGC.AGCGCCC.CCCTGCCCcggcgg                                                          >  1:477110/1‑101 (MQ=18)
                                    atTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAA.G.CCAG.AGGGGGGCAGCAC.CC.CCCTGCCCcggcgg                                                          >  1:501559/1‑101 (MQ=17)
                                    anNTGCNGCAAGATGAANGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcgccgg                                                          <  2:1707762/101‑1 (MQ=2)
                                      ttGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCA.AGAGGGGGCAGCG.CC.CCCTGCCCcgccggat                                                        >  2:91130/1‑99 (MQ=2)
                                        gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCg                                                      >  2:3410585/1‑101 (MQ=18)
                                         cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGTGC.AGCGCCC.CCCTGCCCCGGCGGTCCGt                                                     >  1:1490453/1‑101 (MQ=18)
                                         cagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCCGCGGGCCgg                                                     >  2:3847046/1‑100 (MQ=14)
                                          agcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTg                                                    >  2:1075039/1‑101 (MQ=18)
                                            caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTggc                                                  >  2:503747/1‑101 (MQ=17)
                                            caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTggc                                                  >  1:3860326/1‑101 (MQ=18)
                                             aaGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTggcg                                                 >  2:2922315/1‑101 (MQ=18)
                                             aaGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTggcg                                                 <  2:1110206/101‑1 (MQ=17)
                                               gATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAA.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTggcggc                                               <  1:2171140/101‑1 (MQ=18)
                                                 tGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGgcgc                                             <  2:225445/101‑1 (MQ=17)
                                                 tGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCG.CCTCCCTGCCCCGGCGGTCCGTGGCGgggc                                             >  1:414463/1‑101 (MQ=17)
                                                    aaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCtt                                          <  1:1249351/101‑1 (MQ=17)
                                                     aCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg                                         <  1:2335769/101‑1 (MQ=17)
                                                      ccGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTgc                                        >  1:3157711/1‑100 (MQ=17)
                                                       cGGATTCCTTTCATCTTGCCGGAAATACCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc                                       <  2:3625771/101‑1 (MQ=14)
                                                       cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGAGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc                                       >  2:3217971/1‑101 (MQ=18)
                                                       cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.ACG.CCAG.AGGGGGGC.CGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc                                       >  2:2470170/1‑101 (MQ=18)
                                                       cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGGGGCGGCGCCTTGAc                                       >  1:2467268/1‑101 (MQ=17)
                                                       cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGCGAA.TGGC..GAA.G.CCAG.AGGGGGGC.CGCGCCC.CCCTGCCCCGGCGGCCCGGGGCGGGGCCTTgcc                                       >  1:3069832/1‑99 (MQ=12)
                                                        ggATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGaag                                      >  2:590309/1‑99 (MQ=17)
                                                         gATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGGGCCTTGAcgg                                     >  1:1750029/1‑101 (MQ=17)
                                                             ccTTTCATCTTGCCCGAAATAGCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcgg                                 <  2:2365113/101‑1 (MQ=17)
                                                             ccTTTCATCTGGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.GGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcgg                                 <  1:4447110/101‑1 (MQ=17)
                                                                ttCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCggggg                              >  2:109218/1‑101 (MQ=18)
                                                                 tCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGgc                             >  2:940965/1‑101 (MQ=18)
                                                                  cATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGgcg                            >  2:1887942/1‑101 (MQ=18)
                                                                   aTCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCgggggggg                           >  1:3265848/1‑98 (MQ=18)
                                                                   aTCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGgcgc                           >  1:4316494/1‑101 (MQ=18)
                                                                   aTCTTGCCCGAAATATCCCGGGGGGTGAA.TAGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGgcgc                           <  1:621085/101‑1 (MQ=18)
                                                                     cTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAGGGGCGGGGGGGCgg                         >  1:2783674/1‑101 (MQ=17)
                                                                     cTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCgg                         <  2:1427906/101‑1 (MQ=18)
                                                                        gCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCggggggatcggaa                      >  2:3740502/1‑93 (MQ=17)
                                                                          ccGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCAtt                    <  1:3502492/101‑1 (MQ=18)
                                                                          ccGAAATATCCCGGGGGGGGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCAtt                    >  2:2103695/1‑101 (MQ=18)
                                                                           ggAAGTATCCGGGGGGGTGAA.TGGC..GAA.G.CCGG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTa                   <  1:1523195/100‑1 (MQ=21)
                                                                                  tcggggggggTGAA.TGGC..G.AGG.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGaaa            <  1:3249841/97‑1 (MQ=25)
                                                                                    ccGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAAtt          >  2:1775408/1‑101 (MQ=21)
                                                                                    ccGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAAtt          <  2:3162691/101‑1 (MQ=21)
                                                                                     cGGGGGGTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAAACGAAAttc         >  2:683660/1‑100 (MQ=21)
                                                                                      ggggggTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGGGGGGGCGCGGGGATTAACCGAAAATGc        >  1:810231/1‑101 (MQ=18)
                                                                                       gggggTGAA.TGGC..GAA.GCCAG..AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCa       >  2:1844584/1‑101 (MQ=25)
                                                                                       gggggTGAA.TGGC..GAA.GCCAG..AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCa       >  2:762742/1‑101 (MQ=25)
                                                                                        ggggTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCaa      >  1:2041783/1‑101 (MQ=25)
                                                                                          ggTGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACa    <  1:688782/101‑1 (MQ=25)
                                                                                           gTGAA.TGGC..GAA.GCCAG..AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAg   >  1:1934737/1‑101 (MQ=25)
                                                                                            tGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGa  >  1:4331794/1‑101 (MQ=25)
                                                                                            tGAA.TGGC..GAA.G.CCAG.AGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGa  <  2:2936647/101‑1 (MQ=25)
                                                                                                       |                                                                                                 
GCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.A.G.CAGA.AGGGGGGCAAGCG.CC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCAACAGA  >  NC_014034/2493961‑2494152

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