Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,327 +C 6.5% intergenic (‑123/‑182) 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,3273.→C 6.5% 25.3 ‑2.3 ‑0.4154intergenic (‑123/‑182)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (2/8):  ref base (74/70):  total (76/78)
Fisher's exact test strand distribution p-value = 9.82e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.06e-01

Only 100 of 177 total aligned reads displayed.

TTGAGATGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.G.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAC  >  NC_014034/2249228‑2249426
                                                                                                       |                                                                                                       
ttGAGATGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGaa                                                                                                          <  2:272470/101‑1 (MQ=255)
   agaTGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...ggg                                                                                                    >  2:83736/1‑101 (MQ=255)
   agaTGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...ggg                                                                                                    >  2:4263275/1‑101 (MQ=255)
    gaTGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGGAAGATGAA...gggg                                                                                                   >  2:2608384/1‑101 (MQ=255)
     aTGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...ggggg                                                                                                  <  1:906413/101‑1 (MQ=255)
        cacaGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCt                                                                                               <  2:2499306/101‑1 (MQ=255)
          caGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCttt                                                                                             <  2:3655832/101‑1 (MQ=255)
           aGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTc                                                                                            <  1:1594615/101‑1 (MQ=255)
           aGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTc                                                                                            <  1:4216614/101‑1 (MQ=255)
            gATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCa                                                                                           >  2:2338390/1‑101 (MQ=255)
            gATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCa                                                                                           <  2:1854648/101‑1 (MQ=255)
             aTACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCAt                                                                                          >  2:1026100/1‑101 (MQ=255)
             aTACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCAt                                                                                          <  1:546313/101‑1 (MQ=255)
                ccGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCtt                                                                                       <  1:1806711/101‑1 (MQ=255)
                ccGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCtt                                                                                       <  1:3744580/101‑1 (MQ=255)
                 cGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTg                                                                                      <  2:1451728/101‑1 (MQ=255)
                  gCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGc                                                                                     >  1:276229/1‑101 (MQ=255)
                    aGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCg                                                                                   <  2:4228665/101‑1 (MQ=255)
                    aGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCg                                                                                   <  2:905584/101‑1 (MQ=255)
                    aGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCg                                                                                   >  1:3377750/1‑101 (MQ=255)
                     gCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGa                                                                                  <  2:2066088/101‑1 (MQ=255)
                      cTTGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGaa                                                                                 <  1:2927752/101‑1 (MQ=255)
                       ttGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGaaa                                                                                >  2:1538005/1‑101 (MQ=255)
                       ttGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGaaa                                                                                <  2:4225319/101‑1 (MQ=255)
                        tGGGCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGaaaa                                                                               <  1:986743/101‑1 (MQ=255)
                         gggCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAat                                                                              <  2:4468/101‑1 (MQ=255)
                          ggCACGCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCTGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAata                                                                             <  2:1455475/101‑1 (MQ=38)
                          ggCACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAata                                                                             <  1:3156191/101‑1 (MQ=255)
                            cACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATc                                                                           <  2:3982962/101‑1 (MQ=255)
                            cACCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATc                                                                           <  2:2180251/101‑1 (MQ=255)
                             aCCCGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATcc                                                                          <  1:2418784/101‑1 (MQ=255)
                               ccGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCg                                                                        >  1:2346695/1‑101 (MQ=255)
                               ccGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCg                                                                        >  1:4321817/1‑101 (MQ=255)
                                cGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCgg                                                                       >  2:3468556/1‑101 (MQ=255)
                                cGGGGCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCgg                                                                       >  1:1909028/1‑101 (MQ=255)
                                  gggCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCgggg                                                                     <  1:4246974/101‑1 (MQ=255)
                                   ggCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCggggg                                                                    >  2:4366761/1‑101 (MQ=255)
                                    gCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCgggggg                                                                   <  1:2587984/101‑1 (MQ=255)
                                    gCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCgggggg                                                                   <  2:3245345/101‑1 (MQ=255)
                                    gCAAGAAGGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCgggggg                                                                   >  2:1158809/1‑101 (MQ=255)
                                         aagGGGGCTCTCC.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAAt                                                              <  2:951015/101‑1 (MQ=255)
                                         aagGGGGCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAAt                                                              >  2:1626727/1‑101 (MQ=255)
                                            ggggCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGgc                                                           <  1:69665/101‑1 (MQ=255)
                                               gCTCT.GCCCCCGTCCGGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGTCTTTTCATGTTGCCGAAACTATCCCGGTCGCTGAATGgcggc                                                        >  2:295579/1‑99 (MQ=21)
                                               gCTCT.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCa                                                        <  1:1006474/101‑1 (MQ=255)
                                                ctct.GCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAg                                                       >  2:247635/1‑101 (MQ=255)
                                                  tctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATTAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGc                                                      >  1:2175671/6‑101 (MQ=255)
                                                      cccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCaga                                                  <  1:4090986/100‑1 (MQ=255)
                                                       ccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGAAGCAagag                                                 >  2:1764339/1‑101 (MQ=255)
                                                            tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGGTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAgggggg                                            >  1:3163885/1‑101 (MQ=255)
                                                            tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAgggggg                                            >  1:320558/1‑101 (MQ=255)
                                                            tcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAgggggg                                            >  2:2576331/1‑101 (MQ=255)
                                                             cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGc                                           >  1:2453951/1‑101 (MQ=255)
                                                              ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCa                                          <  1:3924492/101‑1 (MQ=255)
                                                               cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGAAGCCAGAggggg                                             <  1:1337573/94‑1 (MQ=18)
                                                                gcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAATGGCGAAGCCAGAgggggg                                            >  2:744283/7‑101 (MQ=18)
                                                                 caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGAAGCCAGAGGGGGGc                                           <  1:2327243/96‑1 (MQ=18)
                                                                 caagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGAAGCCAGAGGGGGGc                                           <  1:2900124/96‑1 (MQ=18)
                                                                 caagcggACTCCCCCCGGGATATTGGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGAAGCCAGAGGGGGGc                                           <  1:3882924/96‑1 (MQ=18)
                                                                   ccGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cc                                   <  1:499750/101‑1 (MQ=255)
                                                                    cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.ccc                                  >  1:515202/1‑101 (MQ=255)
                                                                     gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cccc                                 >  2:630121/1‑101 (MQ=255)
                                                                     gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cccc                                 >  2:508026/1‑101 (MQ=255)
                                                                     gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cccc                                 >  2:110774/1‑101 (MQ=255)
                                                                      ggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.ccccc                                >  1:4337870/1‑101 (MQ=255)
                                                                        aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cccccct                              >  1:9319/1‑99 (MQ=255)
                                                                        aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cccccct                              <  2:3414293/101‑3 (MQ=255)
                                                                         cT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cccccctt                             >  2:1921007/1‑98 (MQ=255)
                                                                          t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.ccccccttc                            >  1:194212/1‑97 (MQ=255)
                                                                          t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.ccccccttc                            >  2:594067/1‑97 (MQ=255)
                                                                            cccccGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.ccccccttcg                           <  1:591924/101‑6 (MQ=255)
                                                                            cccccGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.ccccccttcg                           <  2:3163756/101‑6 (MQ=255)
                                                                             ccccGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTcgc                           <  2:1700602/101‑1 (MQ=255)
                                                                               ccGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCgcgg                         >  1:1603494/1‑101 (MQ=255)
                                                                                 gggATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCgcgggc                       >  1:954626/1‑101 (MQ=255)
                                                                                 gggATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCgcgggc                       >  1:3009782/1‑101 (MQ=255)
                                                                                 gggATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCgcgggc                       >  2:2718678/1‑101 (MQ=255)
                                                                                  ggATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCgcgggcg                      >  1:3601854/1‑101 (MQ=255)
                                                                                   gATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCACGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCgccggcgg                     >  1:3619563/1‑101 (MQ=25)
                                                                                      atTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCcg                  <  1:1026188/101‑1 (MQ=255)
                                                                                       tttGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCcgc                 >  1:814224/1‑101 (MQ=255)
                                                                                            ggCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgc            <  2:144875/101‑1 (MQ=255)
                                                                                              cAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgc          >  2:2892422/1‑101 (MQ=255)
                                                                                               aaGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGCGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgcg         >  2:671280/1‑101 (MQ=255)
                                                                                               aaGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgcg         >  2:434029/1‑101 (MQ=255)
                                                                                               aaGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgcg         >  2:1480789/1‑101 (MQ=255)
                                                                                               aaGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgcg         >  2:1578778/1‑101 (MQ=255)
                                                                                                  aTGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcggcg      >  2:539962/1‑101 (MQ=255)
                                                                                                  aTGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcgccg      <  1:2627658/101‑1 (MQ=255)
                                                                                                  aTGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGAGcgcct      >  2:4391130/1‑100 (MQ=255)
                                                                                                  aTGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcgccg      >  1:4210006/1‑101 (MQ=255)
                                                                                                   tGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcgccgc     >  2:4191414/1‑101 (MQ=255)
                                                                                                   tGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcgccgc     >  2:2393415/1‑101 (MQ=255)
                                                                                                    gAA...GGGGGGCTTTCATATTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCa    >  2:485481/1‑101 (MQ=255)
                                                                                                     aa...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCaa   <  2:1254047/101‑1 (MQ=255)
                                                                                                      a...GGGGGGCTTTCATCTTGCGGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc  <  2:1976457/101‑1 (MQ=255)
                                                                                                      a...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAGCCC.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc  >  1:1524429/1‑101 (MQ=255)
                                                                                                      a...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc  >  2:1113979/1‑101 (MQ=255)
                                                                                                      a...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc  <  1:1887955/101‑1 (MQ=255)
                                                                                                      a...GGGGGCCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc  <  2:4321548/101‑1 (MQ=255)
                                                                                                       |                                                                                                       
TTGAGATGCACAGATACCGCAGCTTGGGCACCCGGGGCAAGAAGGGGGCTCT.G.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAC  >  NC_014034/2249228‑2249426

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