Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,249,367 +T 2.4% intergenic (‑163/‑142) 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,3671.→T 2.4% 12.0 0.0 ‑0.5209intergenic (‑163/‑142)parE/RCAP_rcc02080DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand):  new base (3/2):  ref base (109/95):  total (112/97)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.13e-01

Only 100 of 233 total aligned reads displayed.

GAAGGGGGCTCTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATA  >  NC_014034/2249268‑2249466
                                                                                                                                                                                                                
gaagGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGaa                                                                                                        <  1:3706845/101‑1 (MQ=255)
 aagGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGaag                                                                                                       >  1:1350309/1‑100 (MQ=255)
 aagGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAAt                                                                                                       <  1:2366384/101‑1 (MQ=255)
 aagGGGGCTCGGCCCCCGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAAt                                                                                                       <  1:3553316/101‑1 (MQ=255)
  agGGGGCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATg                                                                                                      <  1:2222555/101‑1 (MQ=255)
  agGCGGCTCTCCCCCCGCCCCGTCCCGGGCCC.CCCCCGGT.ATTTTGGAGGCAAGATGAG...GGGGGCCTTTCTCCTTGCCAAAAATATCCCGGGGGGTGAATg                                                                                                      <  2:883553/101‑1 (MQ=16)
    ggggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGgc                                                                                                    <  1:1297624/101‑1 (MQ=255)
    ggggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGgc                                                                                                    >  2:3586307/1‑101 (MQ=255)
    gggcctcgccccccGTCCCGTGCGGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGgc                                                                                                    <  2:3995811/91‑1 (MQ=37)
      ggCTCTGCCCCCGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGgcgc                                                                                                  <  1:110420/101‑1 (MQ=255)
         tctgcccccGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGc                                                                                               >  1:999976/6‑101 (MQ=255)
         tctgcccccGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGc                                                                                               >  1:264106/6‑101 (MQ=255)
          cttcccccGTCCCGTGCCGGGTCT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGcc                                                                                              <  1:4425168/97‑1 (MQ=255)
            gcccccGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCag                                                                                            >  2:60642/3‑101 (MQ=255)
             cccccGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCaga                                                                                           <  2:1794228/100‑1 (MQ=255)
              ccccGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCagag                                                                                          <  1:3116734/101‑1 (MQ=255)
               cccgtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAgg                                                                                         >  1:2206632/1‑101 (MQ=255)
                 cgtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAgggg                                                                                       >  2:81113/1‑101 (MQ=255)
                 cgtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAgggg                                                                                       >  2:45967/1‑101 (MQ=255)
                 cgtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGGGAATGGCGCAGCCAGAgggg                                                                                       >  2:2556298/1‑101 (MQ=255)
                  gtcccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAggggg                                                                                      <  2:1738376/101‑1 (MQ=255)
                    cccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGc                                                                                    <  2:1634417/101‑1 (MQ=255)
                     ccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCa                                                                                   >  2:2215383/1‑101 (MQ=255)
                     ccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCa                                                                                   >  2:2747181/1‑101 (MQ=255)
                     ccgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGAGGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCa                                                                                   <  2:3214748/101‑1 (MQ=255)
                      cgtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGTGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAg                                                                                  <  2:3509829/101‑1 (MQ=255)
                      cgcaagcggACTCCCCCCGGG.ATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGAAGCCAGAggggg                                                                                      <  1:3922950/94‑1 (MQ=18)
                       gtGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGTAAATATCCCGGGGGGGGAATGGCGCAGCCAGAGGGGGGCAg.c                                                                                >  2:4617883/1‑101 (MQ=255)
                       gcaagcggACTCCCCCCGGG.ATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGAAGCCAGAgggggg                                                                                     <  1:3576625/95‑1 (MQ=18)
                        tGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAg.cg                                                                               <  2:4149620/101‑1 (MQ=255)
                           gcggACTCCCCCCGGG.ATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGAAGCCAGAGGGGGGCAg.c                                                                                <  1:991578/99‑1 (MQ=18)
                           ccattcaccccccGGG.ATATTTTCGGCAAGATGAA...AGCCCCCCTTCATCTTGCCGAAAATATCCCGGGGGGCGAATGGCGCAGCCATAGGGGGGCAG.CG.cc                                                                            >  2:3809676/9‑101 (MQ=11)
                             ggACTCCCCCCGGG.ATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAATGGCGAAGCCAGAGGGGGGCAg.cg.c                                                                             <  1:527047/101‑1 (MQ=17)
                              gACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cc.cccc                                                                       <  1:737818/101‑2 (MQ=255)
                              gACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cc.cccc                                                                       >  1:4029250/1‑100 (MQ=255)
                               aCT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cc.cccct                                                                      >  2:4380080/1‑99 (MQ=255)
                                cT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cc.cccctt                                                                     >  1:358914/1‑98 (MQ=255)
                                 t.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cc.ccccttc                                                                    <  2:3059914/101‑5 (MQ=255)
                                 t.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cc.ccccttc                                                                    <  1:1557905/101‑5 (MQ=255)
                                 cccccccGGGAATTTTGGAGGAAAGAGGAG...GGGGGGTTTTCTTTTTGCAAAAATTACCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cc.cccct                                                                      <  1:3981906/99‑3 (MQ=12)
                                   cccccGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cc.ccccttcg                                                                   <  1:4528552/101‑6 (MQ=255)
                                   cccccGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.cc.ccccttcg                                                                   >  1:812051/1‑96 (MQ=255)
                                     cccGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTcgcg                                                                  >  2:1272730/1‑101 (MQ=255)
                                        ggg.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCgcgggc                                                               >  2:487361/1‑101 (MQ=255)
                                             tatTTGAGGCAAGATGAG...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGcc                                                           <  2:3019066/101‑1 (MQ=255)
                                             tatTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGcc                                                           >  2:2176177/1‑101 (MQ=255)
                                              atTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGccc                                                          >  1:106802/1‑100 (MQ=255)
                                              atTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCcg                                                          >  2:2656498/1‑101 (MQ=255)
                                               tttGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCcgc                                                         >  1:277983/1‑101 (MQ=255)
                                                  gAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGAc                                                      <  1:1953156/101‑1 (MQ=255)
                                                   gggaaaGATGAG...GGGGGGTTTTCATCTTCCCGAAAATACCGGGGGGGGGGAAGGGCGCAGCCAGGGGGGGGCAGCCC.CC.CCCTTCGCGGGCGGCCGCGAcg                                                     <  1:3329894/97‑1 (MQ=255)
                                                      aaaGATGAG...GGGGGGTTTTCATCTTGCCGAAAATATCCGGGGGGGTGAATGGCGCAGCCAGGGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGAcgcgc                                                  <  1:4168614/100‑1 (MQ=25)
                                                       aaGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGAcgcgcg                                                 >  2:3307546/1‑101 (MQ=255)
                                                       aaGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGAcgcgcg                                                 >  2:1824702/1‑101 (MQ=255)
                                                       aaGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGAcgcgcg                                                 >  1:2517704/1‑101 (MQ=255)
                                                        aGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACgcgcgc                                                >  2:3871958/1‑101 (MQ=255)
                                                        aGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACgcgcgc                                                >  1:3305551/1‑101 (MQ=255)
                                                         gATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGcgcc                                               <  1:4351254/101‑1 (MQ=255)
                                                          aTGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGGGGCCGCGACGCGggccg                                              >  1:2914682/1‑101 (MQ=25)
                                                          aTGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGcgccg                                              >  2:3922272/1‑101 (MQ=255)
                                                              a...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                          >  1:2704383/1‑101 (MQ=255)
                                                              a...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                          <  2:244620/101‑1 (MQ=255)
                                                              a...GGGGGGCTTTCATCTCGCAGAAAATAGCCGGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                          <  1:230433/101‑1 (MQ=37)
                                                                  ggggggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                         <  1:2747345/101‑1 (MQ=255)
                                                                  ggggggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                         >  1:2461999/1‑101 (MQ=255)
                                                                   gggggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                        >  2:658265/1‑101 (MQ=255)
                                                                    ggggCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAgg                                       >  2:4531718/1‑101 (MQ=255)
                                                                      ggCTTTGATCTTGCCGAAACTAGCCCGGGGGGTGAATGGCGCAGCCAGGGGGGGGCGG.CGCCC.CCCTTCGCGGGCGGCCGCGAGGCGCGCCGCAACAGGTc                                     <  2:4687873/101‑1 (MQ=25)
                                                                            cATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                               <  2:358880/101‑1 (MQ=255)
                                                                            cATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGccc                               >  2:2856412/1‑101 (MQ=255)
                                                                               cTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                            >  2:3402257/1‑101 (MQ=255)
                                                                               cTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                            >  1:1769054/1‑101 (MQ=255)
                                                                                  gCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGt                         <  1:674676/101‑1 (MQ=255)
                                                                                     gAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAGCCC.CC.CCCTTCGCGGGGGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                      <  1:1468801/101‑1 (MQ=38)
                                                                                     gAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                      <  1:3002827/101‑1 (MQ=255)
                                                                                     gAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAt                      <  1:3448326/101‑1 (MQ=255)
                                                                                      aaaaTATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                     <  1:3643399/101‑1 (MQ=255)
                                                                                      aaaaTATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.CCACCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCAtt                     >  2:1904205/1‑101 (MQ=255)
                                                                                        aaTATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                   <  1:2082068/101‑1 (MQ=255)
                                                                                        aaTATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                   >  1:4276132/1‑101 (MQ=255)
                                                                                         atatCCCGGGGGGGGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTGGTGCCCCTTTGTCCTTCTg                  >  1:422611/1‑101 (MQ=37)
                                                                                             cccGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCCACCGGTCGTGCCCCGTTGTCATTCTGTCGc              >  2:4243977/1‑101 (MQ=38)
                                                                                             cccGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGc              >  1:4316619/1‑101 (MQ=255)
                                                                                                ggggggTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCGGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc           >  2:2810997/1‑101 (MQ=255)
                                                                                                ggggggTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc           >  2:1074043/1‑101 (MQ=255)
                                                                                                 gggggTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAccc          >  2:4075529/1‑101 (MQ=255)
                                                                                                  ggggTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa         <  2:2285914/101‑1 (MQ=255)
                                                                                                   gggTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt        >  2:4208839/1‑101 (MQ=255)
                                                                                                   gggTGAAGGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt        <  1:382215/101‑1 (MQ=255)
                                                                                                    ggTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg       >  2:239263/1‑101 (MQ=255)
                                                                                                    ggTGAATGGCGCAGACGGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCCACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg       <  1:4618136/101‑1 (MQ=255)
                                                                                                     gggAATGGCGCCGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt      <  2:1970990/99‑1 (MQ=255)
                                                                                                     gTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt      >  1:1820113/1‑101 (MQ=255)
                                                                                                     gTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt      >  2:1314970/1‑101 (MQ=255)
                                                                                                     gTGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt      >  1:2161967/1‑101 (MQ=255)
                                                                                                      tGAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTc     >  2:1388287/1‑101 (MQ=255)
                                                                                                       gAATGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa    >  2:4012377/1‑101 (MQ=255)
                                                                                                        aaTGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCAt   >  2:68389/1‑101 (MQ=255)
                                                                                                         aTGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa  <  2:4355028/101‑1 (MQ=255)
                                                                                                         aTGGCGCAGCCAGAGGGGGGCAG.CGCCC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAAAGGTCGTGCCCCTTTGTCTTTCTGTCGCACCCATGTCATa  >  1:1706688/1‑101 (MQ=255)
                                                                                                                                                                                                                
GAAGGGGGCTCTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGG.ATATTTGAGGCAAGATGAA...GGGGGGCTTTCATCTTGCCGAAAATATCCCGGGGGGTGAATGGCGCAGCCAGAGGGGGGCAG.CG.CC.CCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATA  >  NC_014034/2249268‑2249466

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