Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,037,934 T→G 13.9% intergenic (+129/+106) 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,9340T→G13.9% 4.2 ‑4.2 ‑0.054intergenic (+129/+106)torR/dorSTorCAD operon transcriptional regulatory protein TorR/DMSO/TMAO‑sensor hybrid histidine kinase
Reads supporting (aligned to +/- strand):  new base (0/8):  ref base (23/23):  total (23/31)
Fisher's exact test strand distribution p-value = 1.54e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.75e-01

Only 100 of 125 total aligned reads displayed.

GAGGTCATGCGGACCAGAAGTCCCGTGACCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTG  >  NC_014034/3037835‑3038034
                                                                                                               |                                                                                                        
gAGGTCATGCGGACCAGAAGTCCCGTGACCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTg                                                                                                         <  2:2489417/101‑1 (MQ=255)
 gggTCATGCGGACCAGAAGTCCCGTGCCCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTgg                                                                                                        <  1:4089996/100‑1 (MQ=255)
 gggTCATGCGGACCAGAAGCCCCGTGACCGCCG.GGGGCCCAAGG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTgg                                                                                                        <  2:1348898/100‑1 (MQ=255)
  ggccATGCGGCCCAGAAGCCCCGTGCCCGCCG.GGGGCCCAAAG.G.CCCT.C.GCCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTgg.g                                                                                                      <  2:3544476/98‑1 (MQ=17)
   ggcATGCGGACCAGAAGCCCCGTGCCCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.gc                                                                                                     <  1:4523695/99‑1 (MQ=255)
   gTCATGCGGCCCAGAAGTCCCGTGACCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.gc                                                                                                     <  2:1289/101‑1 (MQ=255)
   gTCATGCGGACCAGAAGTCCCGTGACCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGGGGGTTCC.C.G.CGCC.CCGGGGCGgg.g                                                                                                      >  1:3873427/1‑101 (MQ=18)
   gTCAGGCGGACCAGAAGTCCGGTGACCGCCG.GGGGCCCAAGG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.gc                                                                                                     <  1:2460468/101‑1 (MQ=255)
      aggCGGCCCAGAAGTCCCGTGCCCGCGG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CC.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCa                                                                                                  <  1:3428374/99‑1 (MQ=17)
      aggCGGACCAGAAGTCCCGTGACCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCa                                                                                                  <  1:5399572/99‑1 (MQ=255)
      aggCGGACCAGAAGCCCCGTGCCCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCa                                                                                                  <  1:2289967/99‑1 (MQ=255)
      aTGCGGACCAGAAGTCCCGTGACCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGGGGGGTCG.C.G.C.CGCCCGGGGCGgg.gggg                                                                                                   >  1:2978873/1‑98 (MQ=17)
       tgcgtcccAGAGCTCCCGTGACCGCCG.GGGCCCCAACG.C.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAg                                                                                                 <  2:819108/95‑1 (MQ=17)
       tGCGGACCAGAAGCCCCGTGACCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAg                                                                                                 <  1:738136/101‑1 (MQ=255)
         cGGACCAGGAGTCCGGGGCCCGCGG.GGGCCCCAAAG.G.CCCC.C.GGCCTG.CC.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGc                                                                                               <  1:4599634/101‑1 (MQ=14)
          ggACCAGAAGTCCCGTGACCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGGGCGCTCG.C.G.C.CC.CCCGGGCGGG.GCGCGgggc                                                                                              >  2:1785286/1‑99 (MQ=17)
           gcccAGAAGCCCCGTGCCCGCGG.GGGCCCCAAGG.C.CCCT.C.GGCCGG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCct                                                                                             <  2:300534/99‑1 (MQ=14)
           gcccAGAAGCCCCGCGCCCGCCG.GGGCCCCAAGG.G.CCCT.C.GCCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCct                                                                                             <  2:5015901/99‑1 (MQ=14)
           gACCAGAAGCCCCGGGACCGCGG.GGGGCCCAAGG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCct                                                                                             <  2:795850/101‑1 (MQ=17)
            cccAGAAGCCCGGTGCCCGCGG.GGGCCCCAAGG.C.CCCT.C.GGCCTG.CC.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCctc                                                                                            <  1:373715/100‑1 (MQ=14)
              cAGAGGTCCTGTGACCGCCG.AGGGCCCAAAG.G.CACT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGGGCCGGCG.C.G.C.CC.CCCCGGCGGG.GGGCGggggcccc                                                                                          >  2:2151884/1‑95 (MQ=17)
              cAGAAGTCCCGTGACCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGGGCCGTCG.C.G.C.CC.CCCGGGGTCG.GCGCAGGGctccg                                                                                          >  2:4997362/1‑99 (MQ=17)
               aGAAGTCCCGTGACCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGGGGGTTCGCC.G.C.CC.CCGGGGGGGG.GGGCGGGcgccct                                                                                          >  1:5254120/1‑96 (MQ=14)
               aGAAGTCCCGTGACCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGGGCCAACG.C.C.C.CC.CCCGGGGCGG.GGGCAggggccctc                                                                                         >  2:1551023/1‑94 (MQ=14)
               aGAAGTCCCGTGACCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCACCCG.C.G.CGCC.CCGGGGCGGG.GGGCCgggcccag                                                                                          >  1:734583/1‑95 (MQ=14)
                tttgTCCCATGACCGCCGAGGGCCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCCATGGCGGGGGCATTCC.G.G.C.CC.CCCGGGCGGG.GGGCCgggccctct                                                                                         >  1:2896946/4‑94 (MQ=14)
                  aGTCCCGTGACCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCC.C.G.C.CGCCCGGGGCGGG.GCCCGGGccccctcgg                                                                                       >  2:3566663/1‑94 (MQ=14)
                    tCCCGTGACCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGGGCATTTC.C.G.C.CG.CCGGGGGCGG.GGGGCGGGCTCTGTCGT.G.c                                                                                  >  2:2548801/1‑101 (MQ=17)
                             ccgcg.GGGCCCCAAGG.G.CCCC.C.GGCCCGCCC.CCCCCCCCGCCATCGGCGGGCGCATCCG.C.G.C.CG.CGGGGGCGGG.GCGAAGGCCTCTGTCGT.G.CTTGGGGaa                                                                          <  2:5699079/98‑1 (MQ=9)
                              gccg.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGGGGGTTTG.C.G.C.CG.CCGGGGCGGG.GGGCGGGGCCCTGTCCT.G.CTTtgggggaa                                                                        >  2:2744082/1‑97 (MQ=11)
                              ccgg.ggggCCCAAAGCC.CCCT.G.GCCCGG.CC.CCCCCCCCCCCTGGGGCGGGCTCATTCG.C.G.C.CG.CGGGGGCGGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGaaa                                                                         <  2:5461778/98‑1 (MQ=9)
                                gg.ggggCCACAGG.G.CCCTGC.GGCCGC.CC.CCCGCCCCGCGAGCGGGGGGGGCATTCG.C.G.C.GG.CGGGGGGTGG.GCGAAGGCCTCTGTCGT.G.CTTGGGGAAACg                                                                       <  2:5136178/100‑1 (MQ=255)
                                gg.ggggCAC.AAG.G.CCCC.C.GGCCCC.CCCCCCCCCCCGCCGGCGGCGGGCGCGTCCC.C.C.C.CG.GGGGGGCTGG.GCGAAGGCCTCTGTCGT.G.CTTGGGGAAACgg                                                                      <  1:4596901/100‑1 (MQ=255)
                                cg.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGGGGACTCC.C.G.C.CC.CCCGGGCGGG.GGGGGGGGCCCCGTCGT.G.CTTtgggggacag                                                                      >  1:915015/1‑95 (MQ=2)
                                    gggCCCAAGG.C.CCCC.CGGGCCGG.CC.CCCCCCCCGCCAGGGGGGGGCGCATTCG.C.G.C.CG.CGGGGGCGGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTc                                                                    <  1:2862500/101‑1 (MQ=9)
                                      gCCAAAAG.GCCCCC.C.GCCCGG.CG.CCCGCCCCGCCTGGGGGGGGCGCTTACG.C.G.A.CG.CCGGCGCGGG.GCACAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCtt                                                                  <  2:935857/101‑1 (MQ=2)
                                       cccAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGGGGGCTCG.C.G.C.CG.CCCGGGGCGG.GGGCGGGGCCCTCTCGT.G.GCTGGGGGGACGGGCTTgg                                                                >  2:4402817/1‑101 (MQ=2)
                                          caagccccccC.C.GCCCCG.CC.CCCCCCCCCCCTTCGGCGGGCGCATGCG.C.G.C.CG.CGGGGGCGGG.GCCAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGc                                                               <  2:4202039/94‑1 (MQ=9)
                                          aaaG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGGGGGATCC.G.G.G.CG.CCGGGGCGGG.GCGGGGGGCGCCCTCGT.G.GTTGGGGGAAAggcttttgccc                                                             >  1:797564/1‑92 (MQ=2)
                                             g.g.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGGGGGTTCG.C.G.C.CG.CCCGGGCTGG.GGGCCGGGCCCCTTCGG.G.GTTGGGGGGACAGGGTTGGGACAgg                                                          >  2:899291/1‑101 (MQ=2)
                                                   cT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGGGCCGGCG.G.G.G.CG.CCGGGGCGGG.GGGAAGGGCTCTGTTGT.G.CTTTTGGTAACGGGCGTgggggcggcgcg                                                      >  1:1006312/1‑91 (MQ=2)
                                                   cT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGGGCCCTCG.C.GCC.CC.CCGGGGGGGG.GGGCGGGCC.CTCTCCT.GCCTTGGGGAAA.GGGCTTgggcccgggggg                                                      >  2:4133150/1‑91 (MQ=0)
                                                        ggCCTG.CG.CCCGCCCCGCCATCGGCGGGGGCATTCG.C.G.C.CGCCCGGGGCGGG.GGGCGGGCC.CCGTCGG.G.CTTTGGGGAAAGGGCCTGGCACCGGGG.ATggg                                                  >  2:4213055/1‑99 (MQ=2)
                                                        cgggcgcgcccccgccccgcCAGGGGCGGGCGCATTCG.C.G.C.CG.CGGGGGCGGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATg                                                    <  1:1317878/91‑1 (MQ=18)
                                                        cgcccggcccccccccccgcCGGGGGGGGGCGCATGCG.C.G.C.GG.GGGGGGGGGG.GCGAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATg                                                    <  2:2978706/91‑1 (MQ=16)
                                                           cccccccccccccccgcCAGGGGGGGGCGCTTCCG.C.G.C.CG.GGGGGGGGGG.CCCAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATgcgc                                                 <  1:3609614/94‑1 (MQ=16)
                                                           cTG.CG.CCCGCCCCGCCATCGGCGGGGGCATTCG.C.G.C.CG.CCGGGGGTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGGAAAGGGCTTGGCACAGGGG.ATgagcga                                               >  2:901943/1‑99 (MQ=18)
                                                            ccccccccccccccgcCGGCGGGGGGCGCATCCG.C.G.C.CG.CGGGGGCGGG.GCAAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGcc                                                <  1:1904367/95‑1 (MQ=18)
                                                            ccccccccccccccgcCAGCGGCGGGCGCTTGCG.C.G.G.CG.GGGGGGGGGG.GCGAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGcc                                                <  1:4271028/95‑1 (MQ=17)
                                                             ccgccccccccccgcCTGCGGGGGGCGCGTTCG.C.G.C.CG.GGGGGGCGGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCa                                               <  2:733418/96‑1 (MQ=18)
                                                                g.cccgcccCGCCATCGGCGGGCGCATTCG.C.G.C.CC.CCCGGGCGGG.GGGCAGGCCTCTGTCGT.G.CTTGGGGGAACGGGCTTGGCCCAGGGG.ATGCGCCAgccg                                           >  2:2407027/1‑99 (MQ=17)
                                                                cccccgccccgcCGGGGGGGGGCGCTTTCG.G.G.C.GG.CGGGGGCGGG.GCGAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAgcg                                            <  2:4106114/99‑1 (MQ=18)
                                                                 cccccccccgcCGGGGGGGGGCCCGTCCG.C.G.C.GG.GGGGGGGGGG.GCAAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAgcgg                                           <  1:4117262/100‑1 (MQ=255)
                                                                 cccccccccgcCCTGGGGGGGCGCATCCG.C.G.C.CG.CGGGGGGGGG.GCAAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAgcgg                                           <  1:2457791/100‑1 (MQ=21)
                                                                  cccgccccgcCATCGGCGGGGGCATTTG.CGG.C.CC.CCGGGGGGGG.GGGCGGGCC.CTGTCGTGG.TTTGGGGAAAGGGGC.TGGGACAGGGG.GTGCGCCAgcggg                                          >  2:111988/1‑101 (MQ=9)
                                                                      ccccgccATCGACGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGcgggcggg                                      >  1:5679891/1‑101 (MQ=255)
                                                                       cccgccATCGGGGGGCGCATCCG.C.G.C.GG.CGGGGGCGGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCgggg                                     <  1:5705646/101‑1 (MQ=21)
                                                                       cccgccATCGGCGGGCGCATTCG.C.G.C.CG.CCCGGGCGGG.GGGCCGGCCCCTGTCGT.G.CTTGGGGGAACGGGCTTGGCACAGGGG.ATGCGCCCGCGGGGgggg                                     >  2:3285068/1‑101 (MQ=18)
                                                                        ccgccGTGGGGGGGCGCTTCCG.C.G.C.CG.CGGGGGGGGG.GCGAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGa                                    <  1:410471/101‑1 (MQ=21)
                                                                        ccgccGTCGGCGGCCGCGTCGG.C.G.C.CG.GGGGGGCTGG.GCGAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGa                                    <  2:3031009/101‑1 (MQ=25)
                                                                        ccgccATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGGTGG.GGGGCGGGCCCCGTCGT.G.CTTAGGGGAACAGTCTTTGGACAGGGG.GTGGGCCCGCGGGGggggg                                    >  2:2275457/1‑100 (MQ=12)
                                                                        ccgccATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGGGGG.GGGCCGGGCCCTGTCGT.G.CTTTGGGGAAAGGGGTTGGCACAGGGG.ATGGGCCCGGGGGGggggg                                    >  2:5062460/1‑100 (MQ=12)
                                                                         cgcccgcggcggGCGCTTTCG.C.G.C.CG.CGGGGGGGGG.GCGAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGaa                                   <  2:914939/95‑1 (MQ=21)
                                                                         cgccATCGGCGGGCGCATTCC.G.G.C.CG.CCCGGGGTGG.GGGCCGGGCTCTTTCGT.G.CTTCGGGGAACGGGCTTGGGCCAGGGGCATG.GCCCGCGGGCggggga                                   >  2:467062/1‑99 (MQ=255)
                                                                         cgccAGGGGGGGGCGCTTGCG.C.G.C.GG.CGGGGGCGGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGaa                                   <  1:4162165/101‑1 (MQ=21)
                                                                          gccGTCGGCGGTCGCGTTCG.C.G.C.CG.CGGGGGTTGG.CCACAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGaaa                                  <  2:4929597/101‑1 (MQ=25)
                                                                          cccgtgggcgggcgCATTCG.C.G.C.CG.CGGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGaaa                                  <  2:5301357/95‑1 (MQ=255)
                                                                           ccATCGGCGGGCGCATTCG.C.G.C.CG.CCCGGGCCGG.GGGCCGGCCTCTGTCGT.G.CTTGGGGAAACGGGCTTTGCACAGGGG.ATGCGCCCGCGGGCGGGGGAAg                                 >  2:1743828/1‑101 (MQ=21)
                                                                            ctggggcgggcgCCTGGG.C.G.G.CG.CGGGGGCGGG.GCGAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAgg                                <  1:3854159/97‑1 (MQ=25)
                                                                              tCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAggga                              >  2:5102227/1‑99 (MQ=255)
                                                                               cggcggGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GGGCCGGCCTCTGTCGT.G.CTTGGGGAAACGGGCTTGGGACAGGGG.GGGCGCCAGGGGGCGGGGGAAGGCAg                             >  2:2526442/1‑101 (MQ=21)
                                                                               cggcggGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGGAAGGCAg                             >  2:486122/1‑101 (MQ=255)
                                                                                ggggggcgCAGGCG.C.G.C.GG.CCGGCGGGGG.GCAAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGt                            <  2:4618119/98‑1 (MQ=25)
                                                                                ggcgggcgCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGt                            >  1:2040926/1‑101 (MQ=255)
                                                                                 gggggcgCGTCCG.C.G.C.GG.GGGGGGGGGG.GCGAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTg                           <  2:5104480/99‑1 (MQ=25)
                                                                                 gcgggcgCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTg                           >  1:374523/1‑101 (MQ=255)
                                                                                  cgggcgCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGc                          >  1:4027338/1‑101 (MQ=255)
                                                                                  cgggcgCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGAGGGGAAAGGCAGTGc                          >  1:4798794/1‑101 (MQ=255)
                                                                                     gcgcATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCtt                       >  1:318675/1‑101 (MQ=255)
                                                                                      cgcATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAgatcggaa                      >  2:3178452/1‑94 (MQ=255)
                                                                                       gcctgcg.c.g.C.CG.CGGGGGGGGG.GCAAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTc                     <  1:1975823/96‑1 (MQ=37)
                                                                                       gcATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCCGGGAAAGGCAGTGCCTTTc                     >  1:531558/1‑101 (MQ=255)
                                                                                         gtccg.c.g.C.CG.GGGGGGCGGG.GCAAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTccc                   <  1:890997/98‑1 (MQ=37)
                                                                                         aTTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTccc                   >  1:2609584/1‑101 (MQ=255)
                                                                                         aTTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTccc                   >  2:638096/1‑101 (MQ=255)
                                                                                          ttCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCg                  >  2:531405/1‑101 (MQ=255)
                                                                                           gcg.c.g.C.CG.GGGGGGGTGG.GAAAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGa                 <  2:4403320/100‑1 (MQ=37)
                                                                                            cg.c.g.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGAt                >  2:2237664/1‑101 (MQ=255)
                                                                                             cgcgcccggg.ccgGGGCGGGCGCAAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTcc                    <  1:4158828/92‑1 (MQ=21)
                                                                                               c.g.c.gg.cGGGGGCGGG.GGGAAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATcc              <  1:2366470/101‑1 (MQ=25)
                                                                                                 g.c.cg.ccGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAACGGCAGTGCCTTTCCCGATCCg             <  1:3285268/101‑1 (MQ=255)
                                                                                                 g.c.cg.ccGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCg             <  1:945825/101‑1 (MQ=255)
                                                                                                   c.cg.ccgGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGc            <  1:1743927/101‑1 (MQ=255)
                                                                                                     cg.ccgGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGcc           >  2:1907086/1‑101 (MQ=255)
                                                                                                         cgGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGa        >  1:4321054/1‑101 (MQ=255)
                                                                                                          ggggCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAATGGTCTTGGCACAGGTG.ATGCGCCAGCGCGCGGGGAAAGGCAGTGCCTTCCCCGATCCGCCCGaa       >  2:3814066/1‑101 (MQ=255)
                                                                                                          ggggCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGaa       >  1:1486207/1‑101 (MQ=255)
                                                                                                              cTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCt   >  2:1829109/1‑101 (MQ=255)
                                                                                                               tGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTg  >  1:1118923/1‑101 (MQ=255)
                                                                                                               tGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTg  >  1:4519505/1‑101 (MQ=255)
                                                                                                               |                                                                                                        
GAGGTCATGCGGACCAGAAGTCCCGTGACCGCCG.GGGGCCCAAAG.G.CCCT.C.GGCCTG.CG.CCCGCCCCGCCATCGGCGGGCGCATTCG.C.G.C.CG.CCGGGGCTGG.GCGCAGGCCTCTGTCGT.G.CTTGGGGAAACGGTCTTGGCACAGGTG.ATGCGCCAGCGGGCGGGGAAAGGCAGTGCCTTTCCCGATCCGCCCGAAACCTG  >  NC_014034/3037835‑3038034

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