Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,379 G→A 2.4% intergenic (‑175/‑130) parE ← / → RCAP_rcc02080 DNA topoisomerase IV subunit B/HAD superfamily hydrolase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,249,3790G→A 2.4% 17.3 ‑1.5 ‑0.7454intergenic (‑175/‑130)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (8/3):  ref base (227/216):  total (235/219)
Fisher's exact test strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.15e-01

Only 100 of 512 total aligned reads displayed.

TGCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.CC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAT  >  NC_014034/2249279‑2249478
                                                                                                             |                                                                                                       
gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                                         >  2:3917124/3‑101 (MQ=255)
gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCag                                                                                                         >  1:5875644/3‑101 (MQ=255)
    ccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggg                                                                                                     <  1:3352527/101‑1 (MQ=255)
      gtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAggggg                                                                                                   >  2:9158144/1‑101 (MQ=255)
       tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                                                  <  1:4053185/101‑1 (MQ=255)
       tcccgtGCCGGGACT.CCCCCCGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAgggggg                                                                                                  <  1:7246252/101‑1 (MQ=255)
          cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAg                                                                                              <  2:8502363/101‑1 (MQ=255)
           gtGCCGGGACT.CCCCCGGGAAATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAgc                                                                                             >  1:2645482/1‑101 (MQ=255)
           gggCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAgc                                                                                             <  1:7629877/99‑1 (MQ=255)
            caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATACCCCGGGGGGTGAA.TGGCGAA..GCCAGAGGGGGG.c                                                                                                <  1:6454527/96‑1 (MQ=18)
               cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.cc.c                                                                                       <  1:7664704/101‑1 (MQ=255)
                    cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.cc.c.ccctt                                                                                 <  1:3038057/101‑4 (MQ=255)
                     t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.cc.c.cccttc                                                                                >  1:7059912/1‑97 (MQ=255)
                     t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.cc.c.cccttc                                                                                <  2:108225/101‑5 (MQ=255)
                     t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.cc.c.cccttc                                                                                >  1:94288/1‑97 (MQ=255)
                        ccccGGGATATTTGAGGCAAGATGAA...GGTG..TCCTTTCATCTTGCCGAAAATATCCCAGGCGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTcgc                                                                               <  2:7325280/101‑1 (MQ=21)
                         cccGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAAAGGCCAAGAGGGGG.CAGCG.CC.C.CCtgc                                                                                 >  2:8121891/1‑99 (MQ=2)
                          ccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCgcgg                                                                             >  1:7989119/1‑101 (MQ=255)
                            gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCgcgggc                                                                           >  2:3383730/1‑101 (MQ=255)
                             ggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCgcgggcg                                                                          >  1:4991200/1‑101 (MQ=255)
                               atatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGggcggc                                                                        >  2:3121539/1‑101 (MQ=255)
                               atatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGggcggc                                                                        >  2:3720733/1‑101 (MQ=255)
                                tatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGcc                                                                       <  1:3496917/101‑1 (MQ=255)
                                tatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGcc                                                                       >  2:5563506/1‑101 (MQ=255)
                                tatTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGcc                                                                       >  1:201674/1‑101 (MQ=255)
                                 atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCcg                                                                      >  2:4167831/1‑101 (MQ=255)
                                    tGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGa                                                                   <  1:7902126/101‑1 (MQ=255)
                                    tGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGa                                                                   >  2:6930111/1‑101 (MQ=255)
                                     gAGGCAAGATGAA...GGGG..GTCTTTCATCTTGCCGAAAATATCGCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGAc                                                                  <  2:9023061/101‑1 (MQ=37)
                                      aGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGAcg                                                                 >  2:6290077/1‑101 (MQ=255)
                                          aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGCcgcgcg                                                             >  1:2233317/1‑101 (MQ=255)
                                              tGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGcgccgc                                                         >  2:1892084/1‑101 (MQ=255)
                                               gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                        >  2:5877176/1‑101 (MQ=255)
                                               gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                        <  2:8915125/101‑1 (MQ=255)
                                                aa...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                       >  2:8296706/1‑101 (MQ=255)
                                                aa...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCaa                                                       >  1:8437312/1‑101 (MQ=255)
                                                 a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGGGACGCGCGCCGCacc                                                      >  1:8866058/1‑99 (MQ=255)
                                                 a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                      <  2:143857/101‑1 (MQ=255)
                                                 a...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                      <  2:4716052/101‑1 (MQ=255)
                                                      ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                    <  2:2704485/101‑1 (MQ=255)
                                                      ggg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                                    <  2:7194864/101‑1 (MQ=255)
                                                       gg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAgg                                                   <  2:4959999/101‑1 (MQ=255)
                                                       gg..ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAgg                                                   <  1:6664534/101‑1 (MQ=255)
                                                           ggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCCCGACGCCCGCCGCAACAGGTc                                                 >  1:104517/1‑101 (MQ=25)
                                                            gCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCg                                                >  2:2052331/1‑101 (MQ=255)
                                                              tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                              >  2:46374/1‑101 (MQ=255)
                                                              tttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                              <  2:2085355/101‑1 (MQ=255)
                                                               ttCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGc                                             >  1:470451/1‑101 (MQ=255)
                                                                    cTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                        >  2:1207273/1‑101 (MQ=255)
                                                                     ttGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                       <  1:9108057/101‑1 (MQ=255)
                                                                      tGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                      >  2:5994234/1‑101 (MQ=255)
                                                                       gCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGt                                     <  2:2202575/101‑1 (MQ=255)
                                                                        ccGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                                    >  1:745313/1‑101 (MQ=255)
                                                                         cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                                   >  2:6984678/1‑101 (MQ=255)
                                                                         cGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                                   <  2:7407471/101‑1 (MQ=255)
                                                                          gAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                  <  1:3216960/101‑1 (MQ=255)
                                                                          gAAAATATCCCGGGGGGGGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                                  >  2:6147952/1‑101 (MQ=255)
                                                                           aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                                 >  2:1321964/1‑101 (MQ=255)
                                                                           aaaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                                 >  2:6523545/1‑101 (MQ=255)
                                                                            aaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                                <  2:6730736/101‑1 (MQ=255)
                                                                            aaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                                >  2:5333992/1‑101 (MQ=255)
                                                                            aaaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                                >  1:8511560/1‑101 (MQ=255)
                                                                             aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGAGGGCG.CC.C.CTTCCGCGGG.GGGCGCGACGCGCGCCGCAACAGGTCCGGCCCCTTTGTCAttttt                              >  2:4449441/1‑98 (MQ=25)
                                                                             aaTATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                               <  1:168561/101‑1 (MQ=255)
                                                                              atatCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.CC.CACCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGGGCCCCTTTGTCATTGTg                              >  1:4225264/1‑101 (MQ=38)
                                                                              atagCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTg                              <  1:7729471/101‑1 (MQ=255)
                                                                                gtgccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                            <  2:17513/98‑1 (MQ=255)
                                                                                atCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTc                            >  2:4768171/1‑101 (MQ=255)
                                                                                 tCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCg                           <  2:5769311/101‑1 (MQ=255)
                                                                                  cccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGGCGc                          >  1:2953672/1‑101 (MQ=255)
                                                                                   ccGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCa                         >  1:1193650/1‑101 (MQ=255)
                                                                                    cGGGGGCTGAA.TGGGGAA..GCTAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                        <  2:5507599/101‑1 (MQ=255)
                                                                                     ggggggTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                       <  1:3820811/101‑1 (MQ=255)
                                                                                      gggggTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                      >  1:5520198/1‑101 (MQ=255)
                                                                                      gggggTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc                      >  2:3603947/1‑101 (MQ=255)
                                                                                       ggggTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa                     <  2:7714973/101‑1 (MQ=255)
                                                                                       ggggTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa                     <  1:6292409/101‑1 (MQ=255)
                                                                                        gggGGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt                    <  2:2892870/101‑1 (MQ=255)
                                                                                          gTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt                  >  1:3640337/1‑101 (MQ=255)
                                                                                           ttggactgcgcA..TCGAGCGGGAGG.GAGCG.CCGC.TCGTCGGTGGGGGCCGGGAGGCGCGCCGCAACAGGTCGTTCCCCTTTTTCAATCTGTCGCGCCCATGt                  <  2:7961388/94‑1 (MQ=255)
                                                                                            gAA.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa                >  1:4823683/1‑101 (MQ=255)
                                                                                             aa.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt               <  1:5397756/101‑1 (MQ=255)
                                                                                              a.TGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa              >  1:5697934/1‑101 (MQ=255)
                                                                                               agggCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa              <  1:6240183/99‑1 (MQ=255)
                                                                                                tGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa             >  1:4361087/1‑101 (MQ=255)
                                                                                                tGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa             >  2:5111777/1‑101 (MQ=255)
                                                                                                tGGCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa             <  2:7571567/101‑1 (MQ=255)
                                                                                                 tctgcA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg            <  2:2346446/98‑1 (MQ=255)
                                                                                                 ggCGCA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg            <  1:1532694/101‑1 (MQ=255)
                                                                                                 ggCGCA..GCCAGAGGGGGG.CAGCG.CCTC.CCTGCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg            <  2:1532035/101‑1 (MQ=255)
                                                                                                   cgcA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCt          <  2:2604963/101‑1 (MQ=255)
                                                                                                    gcA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAgaga         >  1:2346513/1‑98 (MQ=255)
                                                                                                     cA..GCCAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTgg        <  2:6621588/101‑1 (MQ=255)
                                                                                                          ccAGAGGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCg     >  1:8743714/1‑101 (MQ=255)
                                                                                                            agagGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa   >  1:5232983/1‑101 (MQ=255)
                                                                                                            agagGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa   >  1:1882767/1‑101 (MQ=255)
                                                                                                            agagGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCa   >  2:8654391/1‑101 (MQ=255)
                                                                                                             gagGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt  >  1:2849322/1‑101 (MQ=255)
                                                                                                             gagGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt  >  1:6753273/1‑101 (MQ=255)
                                                                                                             gagGGGGG.CAGCGCCC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAt  >  1:1410492/1‑101 (MQ=255)
                                                                                                             |                                                                                                       
TGCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGG.CAGCG.CC.C.CCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGGCCGCAT  >  NC_014034/2249279‑2249478

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