Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,037,953 T→C 3.6% intergenic (+148/+87) 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,9530T→C 3.6% 9.7 ‑2.9 ‑0.2223intergenic (+148/+87)torR/dorSTorCAD operon transcriptional regulatory protein TorR/DMSO/TMAO‑sensor hybrid histidine kinase
Reads supporting (aligned to +/- strand):  new base (0/8):  ref base (80/135):  total (80/143)
Fisher's exact test strand distribution p-value = 5.31e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.96e-01

Only 100 of 262 total aligned reads displayed.

CCGTGACCGCCGGGGGCCCA.AAG.G.CCCTC.GG.CCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCGCGCCGCCGGGGCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCCGCGATTT  >  NC_014034/3037857‑3038053
                                                                                                       |                                                                                                    
 ccggccccgcgGGGGCCCA.AAGCC.CCCGC.GG.CCGG.CC.CCCCCCCCGCCATGGGCGGGCGCATTCGCGCGGCGGGGGGGGGGCGCAGGCCTCTGTCGTGCtt                                                                                                  <  2:2588040/92‑1 (MQ=9)
  gTGACCGCCGGGGGCCCA.AAG.G.CCCTC.GG.CCTG.CG.CCCGCCCCGCCATCGGCGGGGGCATTCGCGCCCCCCGGGCTGGGGGCAGGGCCCTGTCGTggctgg                                                                                                >  2:2870292/1‑96 (MQ=17)
  ccggacggcgGGGGCCCA.AGG.C.CCCCC.GGCCCCC.CC.CCCCCCCCGCCCGGGGCGGGCGCATGCGCGCCGCGGGGGCGGGGCGCAGGCCTCTGTCACGCTTg                                                                                                 <  2:2422478/93‑1 (MQ=2)
   gggcccggcggggCCCA.AGG.G.CCCCC.GGCCCCC.CC.CCCCCCCCGCCGTGGGCGGGCGCATTCGCGCCGCGGGGGCTGGGCGCAGGCCTCTGTCGTGCTTgg                                                                                                <  1:4141022/92‑1 (MQ=9)
    ggcccgcgGGGGCCCC.AAG.G.CCCCCTGG.CCCG.CGCCCCCCCCCGCCATCGGGGGGCGCATCCGCGCCGCCGGGGCTGGGCGCAGGCCTCTGTCGTGCTTgg                                                                                                <  2:3668158/95‑1 (MQ=17)
    gcccgccgGGGGCCCA.AAGCC.CCCTG.GG.CCCG.CC.CCCGCCCCGCCAGGGGCGGGCGCTTTCGCGCCGCGGGGGCTGGGCGCAGGCCTCTGTCGTGCTTggg                                                                                               <  1:7093808/99‑1 (MQ=14)
    gACCGCCGGGGGCCCA.AAG.G.CCCTC.GG.CCTG.CG.CCCGCCCCGCCATCGGCGGGGGCGATCGCGCCGCCCGGGGGGGGGGGGGGGCGCCGTCGTGCGTgggg                                                                                              >  1:5994700/1‑101 (MQ=2)
    gACCGCCGGGGGCCCA.AAG.G.CCCTC.GG.CCTG.CG.CCCGCCCCGCCATCGGCGGGGGCAATCGCGCCGCCCGGGGGGGGGGGGGGGCCCTGTCGGGCGTgggg                                                                                              >  1:2666273/1‑101 (MQ=11)
      ccggcgGGGGCCCA.GAG.GCCCCGC.GG.CCCG.CC.CCCGCCCCGCCATGGGCGGGCGCATGCGCGCCGCGGGGGCGGGGCGAAGGCCTCTGTCGTGCTTGGGGa                                                                                             <  2:4822410/101‑1 (MQ=255)
       cgcggGGGGCCCA.AGG.C.CCCCC.GGCCCCC.CC.CCCCCCCCGCCGGGGGGGGGCGCTTCCGCGCCGCGGGGGGGGGGCGCAGGCCTCTGTCGTGCTTGGGGaa                                                                                            <  1:2029631/101‑1 (MQ=9)
       ccgggggggCCCA.AGG.G.CCCGCGGG.CGCC.CC.CCCGCCCCGCCGGGGGGGGGCGCATCCGCGCCGGGGGGGCGGGGCGAAGGCCTCTGTCGTGCTTGGGGaa                                                                                            <  2:6974498/97‑1 (MQ=255)
       cccgcggggCCCAGAGG.G.CCCCC.GC.CCCGCCG.CCCCCCCCGCCGTGGGGGGGCGCGTCCGCGCGGCGGGGGGGGGGCGAAGGCCTCTGTCGTGCTTGGGGa                                                                                             <  2:6090839/96‑1 (MQ=255)
          cgGGGGCCCA.AAG.G.CCCTC.GG.CCTG.CG.CCCGCCCCGCCATCGGCGGGCGGCCTCGCGCCCCCCGGGCGGGGGGCGGGGCCCCGCCGGGCTTggggggaagg                                                                                        >  1:3848842/1‑95 (MQ=2)
           gggggCCCA.AAG.G.CCCTC.GG.CCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCGCGCCCCCCGGGCCGGGGGCGGGGCCCTCTCGTGCGTGGGGGAACggg                                                                                       >  2:60272/1‑100 (MQ=14)
            ggggCCCA.AAG.G.CCCTC.GG.CCTG.CG.CCCGCCCCGCCATCGGCGGGGGCCTTCGCGCCGCCCCGGCGGGGGGCAGGCCTCTCTCGTGCGTggggggaacgcc                                                                                      >  2:5334187/1‑93 (MQ=14)
              ggCCCA.AGG.C.CCCCC.GG.CCCGCCC.CCCCCCCCCCCAGGGGCGGGCGCATTCGCGCCGCCGGGGCGGGGCGAAGGCCTCTGTCGTGCTTGGGGAAACGGTCt                                                                                     <  2:2358447/101‑1 (MQ=9)
               cgccA.ACG.G.CCCTC.CG.CCCC.CCCCCCGCCCCGCCGGGGGCGGCCGCTTTCGCGCCGCGGGGGCGGGGCGAAGGCCTCTGTCGTGCTTGGGGAAACGGTCtt                                                                                    <  2:5929443/99‑1 (MQ=9)
                cccA.AAG.G.CCCTC.GG.CCTG.CG.CCCGCCCCGCCATCGGCGGGGGCCTTCGCGCCCCCCGGGCGGGGGGGCGGGCCCTCTCGGGCTTGGGGGGACGGGCTTgg                                                                                  >  1:1752821/1‑101 (MQ=2)
                     aaG.G.CCCTC.GG.CCTG.CG.CCCGCCCCGCCATCGGCGGGGGCACTCGCGCCCCCCGGGCTGGGGGCCGGCCCCCGTCGGGCCTGGGGGAACgggcgttgggca                                                                              >  2:4999602/1‑91 (MQ=9)
                       g.g.CCCTC.GG.CCTG.CG.CCCGCCCCGCCATCGGCGGGCGGATTCGCGCCGCCGGGGGGGGGGGGCGGGCCCTGTCGTGCTTTGGGGAACGGGGTTGGCGCAgg                                                                            >  1:5230537/1‑101 (MQ=9)
                               c.GG.CCTG.CG.CCCGCCCCGCCATCGGCGGGGGCCTTCGCGCCGCCCGGGGGGGGGGCCGGCCCCTGTCGTGCCTGGGGGGAAGGGGTTGGCACAgggggtgg                                                                      >  2:455717/1‑95 (MQ=9)
                                  g.CCTG.CG.CCCGCCCCGCCATCGGCGGGGGCATTCGCGCCGCCCGGGGTGGGGGCAGGGCCCCGTCGTGGTTGGGGGAAAGGGCTTGGGGCAGGGGATgggc                                                                    >  2:6195167/1‑101 (MQ=12)
                                   cgcccccccccccccccccCCGCGGCGGGCCCTTCGGCGCCGCGGGGGGGGGGCACAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATgc                                                                      <  1:1328688/92‑1 (MQ=12)
                                     gcccgg.cccccccGCCCATCGGGGGGCGCTTTCGCGCCGCGGGGGGGGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGcc                                                                   <  2:166538/96‑1 (MQ=18)
                                      ccccccccccccccccCGGCGGGGGGCGCTTCCGCCCCGCGGGGGGGGGGCAAAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGcc                                                                   <  1:6556879/95‑1 (MQ=16)
                                       cccccccccccccgcCGTGGGCGGGCGCTTCCGCGCCGGGGGGGGGGGGCGAAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCa                                                                  <  1:3002562/96‑1 (MQ=18)
                                         cccccccccccgcCAGGGGGGGGCGCATTCGCGCGGCGGGGGCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAgc                                                                <  2:910083/98‑1 (MQ=21)
                                          g.cccgccaCGCCATCGGCGGGGGCATTCGCGCCGCCCGGGGTGGGGGCAGGGCCCTGTCGTGCTTGGGGGAAAGGGCTTGGGCCCGGGGATGGGGCAgcgg                                                              >  2:4214593/1‑101 (MQ=255)
                                               cccccgccGTCGGGGGCCGCTGCCGGGCGGGGGGGGCGGGCCCAAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAgcgggcgg                                                          <  1:2870401/100‑1 (MQ=16)
                                               cccccgccGTCGGCGGGCGCTCCCGCGCCGCGGGGGCGGGGCAAAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAgcgggcgg                                                          <  1:3637955/100‑1 (MQ=18)
                                                ccccgccGGCGGCGGCCGCATCCGCGCGGCGGGGGGGGGGCAAAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGcgggcggg                                                         <  2:3607811/101‑1 (MQ=21)
                                                 cccgccTTGGTCGGGCGCGTTCGCGGGGCGGGGGCGGGGCCCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCgggg                                                        <  1:3455725/101‑1 (MQ=21)
                                                   cgcggggggcgggcgGTTCCGCGCGGGGGGGGCGGGGCGAAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGaa                                                      <  1:4700951/94‑1 (MQ=21)
                                                   cgccAGGGGGGGGCGCATCCCGGCGGCCGGGGTGGGGCACAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGaa                                                      <  2:1082927/101‑1 (MQ=25)
                                                    gccgcgggggggcgCTTCCGCGCGGGGGGGGCGGGGCAAAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGaaa                                                     <  1:1455213/92‑1 (MQ=21)
                                                     ccgggggcggccgCATCCGCGCCGCGGGGGCGGGGCGAAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAg                                                    <  2:4857564/96‑1 (MQ=25)
                                                       gtgggggggcgCGTTCGCGCCGCGGGGTGGGGGCGAAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGc                                                  <  2:2725239/95‑1 (MQ=25)
                                                        ggggcgggcgCGTCCGCGCGGCGGGGGGGGGGCAAAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCa                                                 <  2:4043997/99‑1 (MQ=25)
                                                         cggcggGCGCATTCGCGCCGCCGGGGCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGGAAGGCAg                                                >  2:1708429/1‑101 (MQ=255)
                                                         cggcggGCGCATTCGCGCCGCCGGGGCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAg                                                >  1:4506481/1‑101 (MQ=255)
                                                            cgggcgCATCGGCGCCGCGGGGGGGGGGCACAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGc                                             <  1:835378/101‑1 (MQ=25)
                                                             gggcgCATTCGCGCCGCCGGGGCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGcc                                            >  2:6875855/1‑101 (MQ=255)
                                                               gcgcATTCGCGCCGCCGGGGCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCtt                                          >  2:4365275/1‑101 (MQ=255)
                                                               gcccggccgccCCGCGGGGGGGGGGCAAAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCtt                                          <  1:4285612/94‑1 (MQ=25)
                                                                   ctccgcgGCGCGGGGGCGGGGCGAAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTccc                                      <  1:3352824/98‑1 (MQ=37)
                                                                   aTTCGCGCCGCCGGGGCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTccc                                      >  2:1142596/1‑101 (MQ=255)
                                                                   aTTCGCGCCGCCGGGGCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTccc                                      >  1:1803954/1‑101 (MQ=255)
                                                                    tgcgggccgcgGGGGCGGGGCAAAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCg                                     <  2:1942687/96‑1 (MQ=25)
                                                                     tcccccccccgGGGGGGGGCGCCGGCCTCTGCTGCGCCTGGGGAAACGGTCTTGGCACCGCCCGGGCTCCGGCGGGCGGGGAAAGGCAGTGCCTTTCCTGa                                    <  1:7002272/93‑1 (MQ=25)
                                                                      ggggccgggggggCGGGGCAAAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGAt                                   <  2:3908434/92‑1 (MQ=25)
                                                                         gccgccGGGGCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCg                                >  1:2325222/1‑101 (MQ=255)
                                                                           cgcggGGGCGGGGCACAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGcc                              <  1:5411448/101‑1 (MQ=255)
                                                                           atgtggggCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGcc                              >  2:2819751/5‑101 (MQ=255)
                                                                            gccgGGGCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGCGCCTTTCCCGATCCGccc                             >  2:6452168/1‑101 (MQ=255)
                                                                              cgGGGCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGa                           >  2:3409434/1‑101 (MQ=255)
                                                                               ggggCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGaa                          >  2:2687340/1‑101 (MQ=255)
                                                                               ggggCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGaa                          >  2:1101180/1‑101 (MQ=255)
                                                                                gggCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGaaa                         >  1:6102738/1‑101 (MQ=255)
                                                                                gggCGGGGCGAAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGaaa                         <  2:3543005/101‑1 (MQ=255)
                                                                                 ggCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGTCCGAAAc                        >  1:3981220/1‑101 (MQ=255)
                                                                                 ggCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAAc                        >  1:3239912/1‑101 (MQ=255)
                                                                                 ggCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAAc                        >  1:3687282/1‑101 (MQ=255)
                                                                                    tGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTg                     >  2:410066/1‑101 (MQ=255)
                                                                                    tGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTg                     >  1:1866612/1‑101 (MQ=255)
                                                                                    tGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCCCCCGAAACCTg                     >  1:5832249/1‑101 (MQ=255)
                                                                                    ggggcaaaGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTg                     <  2:6568422/94‑1 (MQ=255)
                                                                                    gggacgcAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTg                     <  2:564212/97‑1 (MQ=255)
                                                                                     gggCGCAGGCCTCTGTCGTGCTTGGGCAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTgg                    <  1:1142644/101‑1 (MQ=255)
                                                                                       gcgcAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTgggg                  >  1:5922047/1‑101 (MQ=255)
                                                                                        gaaaGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGa                 <  1:267791/98‑1 (MQ=255)
                                                                                        cgcAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGa                 >  2:1493591/1‑101 (MQ=255)
                                                                                         gcAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGAc                >  2:6711220/1‑101 (MQ=255)
                                                                                         aaaGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGAc                <  2:5114358/99‑1 (MQ=255)
                                                                                          cAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACt               <  1:1084318/101‑1 (MQ=255)
                                                                                           aGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTc              >  1:5659262/1‑101 (MQ=255)
                                                                                           aGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTc              <  1:682049/101‑1 (MQ=255)
                                                                                           aGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTc              <  2:1843203/101‑1 (MQ=255)
                                                                                           aGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTc              >  2:1339270/1‑101 (MQ=255)
                                                                                           aGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTc              <  1:981185/101‑1 (MQ=255)
                                                                                            ggCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCa             <  1:4857760/101‑1 (MQ=255)
                                                                                            ggCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCa             <  2:4633510/101‑1 (MQ=255)
                                                                                             gCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAg            >  1:4812435/1‑101 (MQ=255)
                                                                                               ctctGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAgag          >  1:436411/1‑99 (MQ=255)
                                                                                               ctctGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGcc          <  1:5553523/101‑1 (MQ=255)
                                                                                               ctctGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGcc          >  2:6358766/1‑101 (MQ=255)
                                                                                               ctctGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGcc          >  1:662418/1‑101 (MQ=255)
                                                                                               ctctGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGcc          <  1:1843870/101‑1 (MQ=255)
                                                                                                tctGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGccc         >  2:4703718/1‑101 (MQ=255)
                                                                                                 ctGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCcg        <  2:5168069/101‑1 (MQ=255)
                                                                                                 ctGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCcg        <  2:5103450/101‑1 (MQ=255)
                                                                                                  tGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCcgc       >  1:6311061/1‑101 (MQ=255)
                                                                                                   gTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCcgcg      >  1:2162588/1‑101 (MQ=255)
                                                                                                     cGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCCGCGAt    <  1:1622675/101‑1 (MQ=255)
                                                                                                     cGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCCGCGAt    <  2:3809722/101‑1 (MQ=255)
                                                                                                     cGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGAGCCAGCGGGCGGGGGAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCACCACCCGAt    >  2:2222697/1‑101 (MQ=25)
                                                                                                      gTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCCGCGAtt   >  1:3538157/1‑101 (MQ=255)
                                                                                                      gTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCCGCGAtt   <  2:3239796/101‑1 (MQ=255)
                                                                                                       tGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCCGCGAttt  >  1:1185437/1‑101 (MQ=255)
                                                                                                       tGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCCGCGAttt  <  1:388505/101‑1 (MQ=255)
                                                                                                       tGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCCGCGAttt  <  1:6195486/101‑1 (MQ=255)
                                                                                                       |                                                                                                    
CCGTGACCGCCGGGGGCCCA.AAG.G.CCCTC.GG.CCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCGCGCCGCCGGGGCTGGGCGCAGGCCTCTGTCGTGCTTGGGGAAACGGTCTTGGCACAGGTGATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTGGGGACTCAGCCCGCGATTT  >  NC_014034/3037857‑3038053

Base quality scores:  ATCG  < 3 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG