Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,037,952 G→A 2.7% intergenic (+147/+88) torR → / ← dorS TorCAD operon transcriptional regulatory protein TorR/DMSO/TMAO‑sensor hybrid histidine kinase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,037,9520G→A 2.7% 7.7 ‑2.4 ‑0.2219intergenic (+147/+88)torR/dorSTorCAD operon transcriptional regulatory protein TorR/DMSO/TMAO‑sensor hybrid histidine kinase
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (79/134):  total (79/140)
Fisher's exact test strand distribution p-value = 8.96e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.21e-01

Only 100 of 256 total aligned reads displayed.

CCGTGACCGCCGGGGGCCCA.AAG.G..CCCTC.GG.CC.TG.CG.CCCGCCCCGCCAT.CGGCGGGCGCATTCGCGCCG.CCGGGGCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCCGCGATT  >  NC_014034/3037857‑3038052
                                                                                                           |                                                                                                     
 cggggcccgggggggCCCA.AGC.G..CCCCC.CG.CC.CGCCC.CCCCCCCCCCCCG.GGGGGGGCCCATTCCCGCCG.CGGGGGGGGGGCGCAGGCCTCTGT.CGTGC.tt                                                                                                 <  1:3167732/91‑1 (MQ=2)
 ccggccccgcgGGGGCCCA.AAGCC..CCCGC.GG.CC.GG.CC.CCCCCCCCGCCAT.GGGCGGGCGCATTCGCGCGG.CGGGGGGGGGGCGCAGGCCTCTGT.CGTGC.tt                                                                                                 <  2:2588040/92‑1 (MQ=9)
 ccccccccgcgGGGGCCAA.AGG.GC.CCCGC.GC.CC.GC.CC.CCCGCCCCGCCCA.GGGCGGGCGCATTCGCGCCG.CCGGGGGTGGGCGCAGGCCTCTGT.CACGC.tt                                                                                                 <  2:4100333/92‑1 (MQ=9)
  gccgcgcgggggggCCCG.CGG.C..CCCTG.GC.CC.GG.CG.CCCGCCCCCCCCG.GGGGGGGCGCATTCGCGCCG.CGGGGGGTGGGCGCAGGCCTCTGT.CACGC.TTgc                                                                                               <  2:2773928/92‑2 (MQ=9)
  gTGACCGCCGGGGGCCCA.AAG.G..CCCTC.GG.CC.TG.CG.CCCGCCCCGCCAT.CGGCGGGGGCATTCGCGCCC.CCCGGGCTGGGGGCAGGGCCCTGT.CGTggctgg                                                                                                >  2:2870292/1‑96 (MQ=17)
  ccggacggcgGGGGCCCA.AGG.C..CCCCC.GGCCC.CC.CC.CCCCCCCCGCCCG.GGGCGGGCGCATGCGCGCCG.CGGGGGCGGGGCGCAGGCCTCTGT.CACGC.TTg                                                                                                <  2:2422478/93‑1 (MQ=2)
   gggcccggcggggCCCA.AGG.G..CCCCC.GGCCC.CC.CC.CCCCCCCCGCCGT.GGGCGGGCGCATTCGCGCCG.CGGGGGCTGGGCGCAGGCCTCTGT.CGTGC.TTgg                                                                                               <  1:4141022/92‑1 (MQ=9)
    gACCGCCGGGGGCCCA.AAG.G..CCCTC.GG.CC.TG.CG.CCCGCCCCGCCAT.CGGCGGGGGGCTTCGCCCCG.CCCCGGGCGGGGGCGGGGCCCCCT.CGTGG.TTgggg                                                                                             >  2:4867699/1‑101 (MQ=2)
     gggcccgGGGGCCCACAAG.GCCCCCGC.GG.CC.CC.CC.CCCCCCCCCCCGT.GGGGGGGCGCTTCCGCGCCG.CGGGGGGGGGGCGCAGGCCTCTGT.CGTGC.TTgg                                                                                               <  2:5565444/97‑1 (MQ=9)
       ccgggggggCCCA.AGG.G..CCCGCGGG.CG.CC.CC.CCCGCCCCGCCGG.GGGGGGGCGCATCCGCGCCG.GGGGGGCGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGaa                                                                                           <  2:6974498/97‑1 (MQ=255)
       cccgcggggCCCAGAGG.G..CCCCC.GC.CC.CGCCG.CCCCCCCCGCCGT.GGGGGGGCGCGTCCGCGCGG.CGGGGGGGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGa                                                                                            <  2:6090839/96‑1 (MQ=255)
          ggggggCCCA.ACC.G..CCCCC.CG.CC.CC.CCCCCCCCCCCCCCGG.CGGGGGGCGCTTTCGCGCGG.CGGGGGGGGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACg                                                                                        <  1:4794598/100‑1 (MQ=9)
          cgGGGGCCCA.AAG.G..CCCTC.GG.CC.TG.CG.CCCGCCCCGCCAT.CGGCGGGGGCATTCGCGCCG.CCCGGGGGGGGGGCGGGGCCCTGT.CCTGG.TTgggggggcgg                                                                                       >  1:3809883/1‑95 (MQ=11)
             gggCCCA.AAG.G..CCCTC.GG.CC.TG.CG.CCCGCCCCGCCAT.CGGCGGGGGCCTTTCCGCCGCCCGGGGCGGGGGGC.GGCCCCCGT.CGGGG.TTGGGGGAAAggggt                                                                                    >  1:2620585/1‑98 (MQ=2)
              ggCCCA.AGG.C..CCCCC.GG.CC.CGCCC.CCCCCCCCCCCAG.GGGCGGGCGCATTCGCGCCG.CCGGGGCGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCt                                                                                    <  2:2358447/101‑1 (MQ=9)
               gCCCC.ACG.G..CCCCC.GG.CCGCG.CC.CCCCCCCCGCCGG.CGGGGGGCGCATTCGCGCCG.CGGGGGGGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCtt                                                                                   <  1:5217043/101‑1 (MQ=9)
               gCCCA.AAG.G..CCCTC.GG.CC.TG.CG.CCCGCCCCGCCAT.CGGCGGGGGGGCTTGCGCCG.CCGGGGGCGGGGGCGGGGCTCCGTCCGTGC.TTtggggaggggcgtg                                                                                   >  2:267787/1‑91 (MQ=2)
               cgccA.ACG.G..CCCTC.CG.CC.CC.CCCCCCGCCCCGCCGG.GGGCGGCCGCTTTCGCGCCG.CGGGGGCGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCtt                                                                                   <  2:5929443/99‑1 (MQ=9)
                ccacgcgg.g..CCCT..GG.CC.CC.CG.CCCCCCCCCCCCTGCGGCGGGCGCATGCGCGCCG.CGGGGGCGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTg                                                                                  <  2:365556/94‑1 (MQ=9)
                    ccagaggccccGC.GC.CC.GG.CC.CCCCCCCCGCCAT.GGGCGGGCGCTTTCGCGCCG.CGGGGGCTGGGCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTgg                                                                                 <  2:5024391/93‑1 (MQ=14)
                       g.g..CCCTC.GG.CC.TG.CG.CCCGCCCCGCCAT.CGGCGGGGGCCTTCGCCCCC.CCGGGGGTGGGGGCCGGGCCCCGG.CCTGC.TTGGGGGGACGGGCTTGGGACCgg                                                                           >  1:1776787/1‑101 (MQ=9)
                                   g.CC.TG.CG.CCCGCCCCGCCAT.CGGCGGGGGCATTCGCGCCG.CCCGGGGTGGGGGCAGGGCCCCGT.CGTGG.TTGGGGGAAAGGGCTTGGGGCAGGGGATgggc                                                                   >  2:6195167/1‑101 (MQ=12)
                                    cggcccccccccccccccgcCGG.GGGCGGGCGCTTGCGCGCGG.GGGGGGGGGGCCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATg                                                                      <  1:2010279/91‑1 (MQ=12)
                                       gcccgg.cccccccGCCCAT.CGGGGGGCGCTTTCGCGCCG.CGGGGGGGGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGcc                                                                  <  2:166538/96‑1 (MQ=18)
                                        tG.CG.CCCGCCCCGCCAT.CGGCGGGCGCATTCGCGCCG.CCGGGGGGGGGCGCCGGCCTCTGT.CGTGC.TTGGGGGAACGGGCTTGGGACAGGGGATGGGCCac                                                                >  1:6015289/1‑100 (MQ=17)
                                        ccgcccccccccccgcCGG.GGGCGGGCGCGGCGGCGCGG.CGGGGGGGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGcc                                                                  <  1:5332485/95‑1 (MQ=12)
                                        ccccccccccccccccCGG.CGGGGGGCGCTTCCGCCCCG.CGGGGGGGGGGCAAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGcc                                                                  <  1:6556879/95‑1 (MQ=16)
                                         cccccccccccccgcCGT.GGGCGGGCGCTTCCGCGCCG.GGGGGGGGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCa                                                                 <  1:3002562/96‑1 (MQ=18)
                                          gccccccgccccgcCAT.CGGCGGGCGCGTGCGCGCGG.GGGGGGGGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAg                                                                <  2:4804106/97‑1 (MQ=18)
                                          gcccccccccccgcCCT.CGGCGGGCGCATGCGCGCCG.CGGGGGCGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAg                                                                <  1:3531685/97‑1 (MQ=18)
                                           gccccccccccccCCG.CGGGGGGCGCATGCGCGCCG.CGGGGGGGGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAgc                                                               <  1:883615/98‑1 (MQ=21)
                                           cg.cccgccCCGCCAT.CGGCGGGCGCAATCGCGCCC.CCCGGGGGGGAGGCCGGGCCCCGG.CGTGC.TTGGGGGAAAGGGCTTGGCACAGGGGGTGCGCCCgcg                                                              >  2:4857516/1‑101 (MQ=9)
                                           cccccccccccgcCAG.GGGGGGGCGCATTCGCGCGG.CGGGGGCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAgc                                                               <  2:910083/98‑1 (MQ=21)
                                            g.cccgcccCGCCAT.CGGCGGGGGCATTCGCGCCG.CCCGGGGGGGGGGCAGGGCCCCGT.CGTGCTTTGGGGAAAGGGCTTTGGC.CACGGGGAGCGCCAgcgg                                                             >  2:1408973/1‑101 (MQ=9)
                                            g.cccgcccCGCCAT.CGGCGGGGGCATTCGCGCCG.CCCGGGCCGGGGGGCGGCCCCTGT.CGTGC.TTTGGGGAAAGGGCGTGGCACAGGGGATGGGCCAgcgg                                                             >  2:5105594/1‑101 (MQ=16)
                                            ccccccccccgcCGG.GGGCGGGCGCCTCCGCGCGG.GGGGGGCGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAgcg                                                              <  1:2768838/99‑1 (MQ=21)
                                               cccccccgccGC.CGGCGGCCCCGCGCGCGCCG.CGGGGGGGGGGAAAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAgcgggc                                                           <  1:1481104/98‑1 (MQ=19)
                                               cccccccgccCG.CGGCGGGCGCTCCCGCGCGC.CGGGGGCGGGGGGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAgcgggc                                                           <  1:5642046/98‑1 (MQ=18)
                                                   cccgccTT.GGTCGGGCGCGTTCGCGGGG.CGGGGGCGGGGCCCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCgggg                                                       <  1:3455725/101‑1 (MQ=21)
                                                      cgcggggggcgggcgGTTCCGCGCGG.GGGGGGCGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGaa                                                     <  1:4700951/94‑1 (MQ=21)
                                                       gccccgggggggcgCGTTCGCGCGG.CGGGGGGGGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGaaa                                                    <  1:1223747/92‑1 (MQ=25)
                                                          gtgggggggcgCGTTCGCGCCG.CGGGGTGGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGc                                                 <  2:2725239/95‑1 (MQ=25)
                                                            gggcgggcgCGCGCGCGCCG.GGGGGGGGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAg                                               <  2:3569870/100‑1 (MQ=25)
                                                            cggcggGCGCATTCGCGCCG.CCGGGGCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGGGGGGGGGGAAAGGGAg                                               >  2:3960337/1‑101 (MQ=255)
                                                            cggcggGCGCATTCGCGCCG.CCGGGGCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGGAAGGCAg                                               >  2:1708429/1‑101 (MQ=255)
                                                            cggcggGCGCATTCGCGCCG.CCGGGGCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAg                                               >  1:4506481/1‑101 (MQ=255)
                                                             ggcgggcgCATTCGCGCCG.CCGGGGCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGt                                              >  1:3850930/1‑101 (MQ=255)
                                                               cgggcgCATCGGCGCCG.CGGGGGGGGGGCACAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGc                                            <  1:835378/101‑1 (MQ=25)
                                                               cgggcgCACGCGCGCCG.GGGGGGCTGGCAGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGc                                            <  2:3288973/101‑1 (MQ=25)
                                                                gggcgCATTCGCGCCG.CCGGGGCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGcc                                           >  2:6875855/1‑101 (MQ=255)
                                                                gggcgCATTCGCGCCG.CCGGGGCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGcc                                           >  1:6112177/1‑101 (MQ=255)
                                                                 gtggcATGCGGGCCG.CGGGGGTTGGGCAAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCt                                          <  2:4195880/98‑1 (MQ=25)
                                                                  gcccggccgccCCG.CGGGGGGGGGGCAAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCtt                                         <  1:4285612/94‑1 (MQ=25)
                                                                     cATTCGCGCCG.CCGGGGCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTcc                                      >  1:5984203/1‑101 (MQ=255)
                                                                      aTTCGCGCCG.CCGGGGCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTccc                                     >  1:1803954/1‑101 (MQ=255)
                                                                       tgcgggccg.cgGGGGCGGGGCAAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCg                                    <  2:1942687/96‑1 (MQ=25)
                                                                       gccccgccg.cGGGGGCGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCg                                    <  2:6440010/97‑1 (MQ=25)
                                                                         tcccccccccgGGGGGGGGCGCCGGCCTCTGC.TGCGC.CTGGGGAAACGGTCTTGGCACCGCCCGGGCTCCGGCGGGCGGGGAAAGGCAGTGCCTTTCCTGa                                   <  1:7002272/93‑1 (MQ=25)
                                                                         cgcgCCG.CGGGGGGGGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGAt                                  <  2:1035904/101‑1 (MQ=255)
                                                                          ccgccg.cGGGGGGGGGGCAAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATc                                 <  2:3848679/100‑1 (MQ=25)
                                                                            gccggggggggCGGGGCAAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATcc                                <  2:5319322/94‑1 (MQ=25)
                                                                               g.ccgGGGCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGCTCCGccc                            >  1:1389549/1‑101 (MQ=255)
                                                                               g.ccgGGGCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGCGCCTTTCCCGATCCGccc                            >  2:6452168/1‑101 (MQ=255)
                                                                                   ggggCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGaa                         >  2:5718504/1‑101 (MQ=255)
                                                                                   ggggCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGaa                         >  2:5133921/1‑101 (MQ=255)
                                                                                   ggggCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGaa                         >  2:1101180/1‑101 (MQ=255)
                                                                                    gggCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGaaa                        >  1:3655715/1‑101 (MQ=255)
                                                                                    gggCGGGGCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGaaa                        <  2:3543005/101‑1 (MQ=255)
                                                                                     ggCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGTCCGAAAc                       >  1:3981220/1‑101 (MQ=255)
                                                                                      ggggggCGAAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAAcc                      <  1:1853194/98‑1 (MQ=255)
                                                                                       gggggaaaaGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCt                     <  2:6633755/93‑1 (MQ=255)
                                                                                        tGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCCCCCGAAACCTg                    >  1:5832249/1‑101 (MQ=255)
                                                                                        gggacgcAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTg                    <  1:5906432/97‑1 (MQ=255)
                                                                                         gggCGCAGGCCTCTGT.CGTGC.TTGGGCAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTgg                   <  1:1142644/101‑1 (MQ=255)
                                                                                           gcaaaGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTgggg                 <  1:3624363/97‑1 (MQ=255)
                                                                                           gcaaaGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTgggg                 <  1:7065874/97‑1 (MQ=255)
                                                                                            gaaaGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGa                <  1:267791/98‑1 (MQ=255)
                                                                                             gcAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGAc               <  2:451751/101‑1 (MQ=255)
                                                                                             gcAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGAc               >  2:6711220/1‑101 (MQ=255)
                                                                                             gcAGGCCTCTGT.CGTGC.GTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGAc               <  2:1139848/101‑1 (MQ=255)
                                                                                              cAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACt              >  1:2138153/1‑101 (MQ=255)
                                                                                              cAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACt              <  1:1084318/101‑1 (MQ=255)
                                                                                              aaGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACt              <  1:3271507/100‑1 (MQ=255)
                                                                                               aGGCCTCTGT.CGTGC.TTGGGGCAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTc             <  2:4575192/101‑1 (MQ=255)
                                                                                               aGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTc             <  2:1843203/101‑1 (MQ=255)
                                                                                               aGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTc             >  2:1339270/1‑101 (MQ=255)
                                                                                               aGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTc             <  1:682049/101‑1 (MQ=255)
                                                                                                ggCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCa            <  1:4857760/101‑1 (MQ=255)
                                                                                                 gCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAg           <  1:6283556/101‑1 (MQ=255)
                                                                                                 ctctctGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAg           <  2:436399/99‑1 (MQ=255)
                                                                                                   ctctGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAgag         >  1:436411/1‑99 (MQ=255)
                                                                                                   ctctGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGcc         >  2:3099594/1‑101 (MQ=255)
                                                                                                   ctctGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGcc         >  1:662418/1‑101 (MQ=255)
                                                                                                   ctctGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGcc         <  1:5553523/101‑1 (MQ=255)
                                                                                                     ctGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCcg       <  2:3688591/101‑1 (MQ=255)
                                                                                                     ctGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCACGATCTGCCCGACACCTAGTGACTCAGCCcg       >  2:4652478/1‑101 (MQ=37)
                                                                                                        t.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCCGCGa    <  2:6031007/101‑1 (MQ=255)
                                                                                                          cGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCCGCGAt   <  1:1622675/101‑1 (MQ=255)
                                                                                                          cGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGAGCCAGCGGGCGGGGGAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCACCACCCGAt   >  2:2222697/1‑101 (MQ=25)
                                                                                                           gTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCCGCGAtt  <  1:2228891/101‑1 (MQ=255)
                                                                                                           |                                                                                                     
CCGTGACCGCCGGGGGCCCA.AAG.G..CCCTC.GG.CC.TG.CG.CCCGCCCCGCCAT.CGGCGGGCGCATTCGCGCCG.CCGGGGCTGGGCGCAGGCCTCTGT.CGTGC.TTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCCGCGATT  >  NC_014034/3037857‑3038052

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG