Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,048 G→A 6.7% intergenic (+96/‑146) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0480G→A 6.7% 20.2 ‑0.6 ‑0.0164intergenic (+96/‑146)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (7/4):  ref base (78/75):  total (85/79)
Fisher's exact test strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.78e-01

Only 100 of 188 total aligned reads displayed.

GTAAGGGGGGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.A.G.CAGAAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCA  >  NC_014034/2493949‑2494147
                                                                                                   |                                                                                                          
taaggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggg                                                                                                           >  2:1381320/5‑101 (MQ=255)
taaggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggg                                                                                                           >  2:4753295/5‑101 (MQ=255)
 aaggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                                                          >  2:3737533/4‑101 (MQ=255)
 aaggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                                                          >  1:347592/4‑101 (MQ=255)
  aggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggg                                                                                                         <  1:4354082/99‑1 (MQ=255)
  aggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggg                                                                                                         >  1:4201577/3‑101 (MQ=255)
   ggggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGACGCCTTTCATCTTGCCCGAAATATCCCgggggg                                                                                                        >  1:4450895/2‑101 (MQ=255)
   ggggggcgCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggg                                                                                                        <  1:783717/100‑1 (MQ=255)
    gggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGCGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGt                                                                                                       >  2:4551779/1‑101 (MQ=255)
     ggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                                                                      >  1:1625525/1‑101 (MQ=255)
     ggggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTg                                                                                                      >  2:3922708/1‑101 (MQ=255)
      gggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAACGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                                                                     >  1:191082/1‑101 (MQ=255)
       ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGaa                                                                                                    <  1:4475775/101‑1 (MQ=255)
       ggggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGaa                                                                                                    >  1:4682235/1‑101 (MQ=255)
         ggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                                                                 <  2:3759132/101‑1 (MQ=255)
         ggCGCTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                                                                 <  2:1976251/101‑1 (MQ=255)
              tGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aa                                                                                         >  2:98397/1‑100 (MQ=255)
                cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcc                                                                                      <  1:1381414/101‑5 (MQ=255)
                cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcc                                                                                      >  2:3871363/1‑97 (MQ=255)
                cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.A.T.CAg                                                                                    <  1:2790042/101‑1 (MQ=255)
                cccccCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagc                                                                                       <  1:4551889/101‑4 (MQ=255)
                cccccCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagc                                                                                       <  1:3577036/101‑4 (MQ=255)
                 ccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcca                                                                                     <  2:2546496/101‑6 (MQ=255)
                    ccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccaga                                                                                   >  2:617701/4‑94 (MQ=255)
                     cgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagag                                                                                  <  2:2067283/99‑9 (MQ=255)
                       tCCGCAAGCGGACTCCCCCCGGGAGATTTGCAGCAAGATGAAACCGGATACCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..G.aagccagaggg                                                                                >  1:43063/1‑91 (MQ=18)
                            aaGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGCa                                                                          <  1:1431276/101‑1 (MQ=18)
                             aGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGCag                                                                         <  1:4753551/101‑2 (MQ=18)
                                ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAC.TGGC..GAA.G.CCAGAGGGGGGCagcgc                                                                      >  1:4638215/1‑97 (MQ=17)
                                  aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGCagcggcc                                                                    >  2:2348947/1‑95 (MQ=17)
                                  aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCG.ccc                                                                  >  1:283443/1‑101 (MQ=18)
                                  aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCG.ccc                                                                  >  2:3025943/1‑101 (MQ=18)
                                   cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.GCCAG.AGGGGGGC.AGCG.cccc                                                                 <  1:402027/101‑1 (MQ=17)
                                   cTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCG.cccc                                                                 <  2:196547/101‑1 (MQ=17)
                                   cTCCCCCCGGGATAGTTGCAGCAAGATGAATCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCG.cccc                                                                 <  2:43056/101‑1 (MQ=17)
                                    tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCG.ccccc                                                                >  2:128706/1‑101 (MQ=17)
                                    tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCG.ccccc                                                                >  1:806307/1‑101 (MQ=17)
                                    tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCG.ccccc                                                                >  1:1565968/1‑101 (MQ=18)
                                      cccccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCG.cccccct                                                              <  2:848858/101‑3 (MQ=17)
                                         ccGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCG.cccccctgcc                                                           <  2:4018385/101‑6 (MQ=17)
                                          cGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCATAGGGGGGC.AGCGCCCCCCTGccc                                                           >  1:2925682/1‑101 (MQ=18)
                                          cGGAATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGccc                                                           <  1:2203745/101‑1 (MQ=17)
                                           gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGcccc                                                          <  2:191057/101‑1 (MQ=17)
                                           gggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGcccc                                                          >  1:3846964/1‑101 (MQ=17)
                                            ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcg                                                         >  1:659612/1‑101 (MQ=2)
                                             gATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgc                                                        >  1:4404457/1‑100 (MQ=2)
                                              atatTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgcc                                                       >  1:2418385/1‑99 (MQ=2)
                                              atatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCcggg                                                       >  1:3337824/1‑100 (MQ=17)
                                                 tttGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgccgga                                                    >  2:2878051/1‑100 (MQ=2)
                                                 tttGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.A.G.GCCAAGAGGGGGCAGCG.CCCCCTGCCCcgccgga                                                    >  2:2605858/1‑100 (MQ=2)
                                                   tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTcc                                                  >  2:263673/1‑101 (MQ=17)
                                                    gcagcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCg                                                 >  1:4463377/1‑101 (MQ=17)
                                                          aGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTggcgg                                           >  2:1161814/1‑101 (MQ=18)
                                                              gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGcc                                       <  2:4421536/101‑1 (MQ=17)
                                                              gAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGcc                                       >  1:2626861/1‑101 (MQ=17)
                                                                aaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCtt                                     <  1:4650818/101‑1 (MQ=17)
                                                                 aCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg                                    >  1:2685816/1‑101 (MQ=18)
                                                                 aCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg                                    <  1:3335693/101‑1 (MQ=17)
                                                                   ccgATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc                                  <  2:1774723/99‑1 (MQ=17)
                                                                   cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc                                  <  2:4640867/101‑1 (MQ=17)
                                                                   cGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAc                                  <  2:1041092/101‑1 (MQ=17)
                                                                    ggATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGCcg                                 >  1:1223598/1‑101 (MQ=17)
                                                                     gATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgg                                <  1:2311185/101‑1 (MQ=17)
                                                                      aTGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcggg                               <  1:1172984/101‑1 (MQ=17)
                                                                      aTGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcggg                               <  2:1430771/101‑1 (MQ=17)
                                                                         ccTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcgg                            <  1:1254947/101‑1 (MQ=17)
                                                                          cTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGAcgggcggg                           >  1:3228256/1‑101 (MQ=18)
                                                                           tttCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTTCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCgggg                          >  1:2587620/1‑101 (MQ=17)
                                                                           tttCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCgggg                          <  1:4605986/101‑1 (MQ=17)
                                                                           tttCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCgggg                          <  2:944117/101‑1 (MQ=17)
                                                                           tttCATCTTCCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCgggg                          <  2:4386003/101‑1 (MQ=17)
                                                                            ttCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCggggg                         >  2:1718966/1‑101 (MQ=18)
                                                                             tCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGgc                        <  2:4272800/101‑1 (MQ=18)
                                                                             tCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGAGGCGGCGCCTTGACGGGCGGGGgc                        >  2:2012970/1‑101 (MQ=17)
                                                                              cATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGgcg                       >  1:3806231/1‑101 (MQ=18)
                                                                               aTCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCGCGGCGGTCCGTGGCGGCGCCTTGACGGGCgggggggg                      >  1:4009401/1‑98 (MQ=17)
                                                                                tCTTTCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGcgcg                     <  2:2176871/101‑1 (MQ=18)
                                                                                 cTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCgg                    >  2:1633663/1‑101 (MQ=18)
                                                                                   tGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGGGCCACGAAGGGCCGGGGCGCgggg                  >  1:835590/1‑100 (MQ=9)
                                                                                   tGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGc                  <  1:617765/101‑1 (MQ=18)
                                                                                   tGCCCGAAATACCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGc                  <  2:4637959/101‑1 (MQ=18)
                                                                                     cccGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCAt                <  2:3538976/101‑1 (MQ=18)
                                                                                     cccGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCAt                <  1:172579/101‑1 (MQ=18)
                                                                                      ccGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGGGCCTTGACGGGCGGGGGCGCGGGCAtt               >  1:2628535/1‑101 (MQ=18)
                                                                                      ccGAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCAtt               <  2:3593736/101‑1 (MQ=18)
                                                                                      ccGAAATACCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCAtt               <  1:2122/101‑1 (MQ=18)
                                                                                        gAAATATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTaa             <  1:98413/101‑1 (MQ=18)
                                                                                         aaaTATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAAc            <  1:459589/101‑1 (MQ=18)
                                                                                          aaTATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAAcc           >  1:1127494/1‑101 (MQ=21)
                                                                                          aaTATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAAcc           >  2:90304/1‑101 (MQ=21)
                                                                                          aaTATCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAAcc           <  1:2740830/101‑1 (MQ=21)
                                                                                            taccccGGGGGGTGAA.TGGC..G.A.A.GCCAGAGGGGGCCAGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGa         <  2:1329142/98‑1 (MQ=21)
                                                                                             atCCGGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGaa        <  2:81842/101‑1 (MQ=21)
                                                                                             atCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGaa        <  1:2278845/101‑1 (MQ=21)
                                                                                             atCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGG.GAGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGGCGCTTGACGGGCGGGGGCGCGGGCATTAACCGaa        >  1:538066/1‑101 (MQ=16)
                                                                                             atCCCGGGGGGGGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGGGCGGGCATTAACCGaa        >  2:2821257/1‑101 (MQ=21)
                                                                                              tCCCGGGGGGTGAA.TGGC..GAA.G.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGaaa       <  1:4150127/101‑1 (MQ=21)
                                                                                                  ggggggTGAA.TGGC..GAA.G.CCAGAGGGCGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGc   >  1:4254827/1‑101 (MQ=25)
                                                                                                   gggggTGAA.TGGC..GAA.GCCAG.AGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCa  >  1:316332/1‑101 (MQ=25)
                                                                                                   gggggTGAA.TGGC..G.AGG.CCAGAGGGGGGC.AGCGCCCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCa  <  2:417052/101‑1 (MQ=25)
                                                                                                   |                                                                                                          
GTAAGGGGGGGCGCTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.A.G.CAGAAGGGGGGCAAGCG.CCCCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGACGGGCGGGGGCGCGGGCATTAACCGAAATTGCA  >  NC_014034/2493949‑2494147

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