Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,209,677 G→C 100% intergenic (+9/‑262) cobO → / → cobQ1 cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,209,6770G→C100.0% 272.1 435intergenic (+9/‑262)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (125/214):  ref base (64/32):  total (189/246)
Fisher's exact test strand distribution p-value = 3.62e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 442 total aligned reads displayed.

ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGC  >  NC_014034/2209579‑2209765
                                                                                                  |                                                                                             
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g                                                                                         >  2:2618127/1‑101 (MQ=255)
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g                                                                                         >  2:682604/1‑101 (MQ=255)
 tCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gg                                                                                        <  1:3207829/101‑1 (MQ=255)
  cGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..ggg                                                                                       >  2:2435032/1‑101 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggg                                                                                      >  1:3421884/1‑100 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggg                                                                                      <  1:1006365/101‑2 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggggggcg                                                                                 >  1:6131230/1‑95 (MQ=255)
           cGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGgggcgctgcc                                                                              <  2:5677108/101‑7 (MQ=255)
           cGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGgggcgctgcc                                                                              >  2:1887153/1‑95 (MQ=255)
             aTCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGcccc                                                                              <  1:8311704/101‑1 (MQ=255)
             aTCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGcccc                                                                              <  1:8483489/101‑1 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.CG.GGGGGCTGCCCCCg                                                                            >  2:1830312/1‑101 (MQ=255)
                tAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGATTTCTGATCTC.AAGGC.CGGGGGCGCTGCCCCCgg                                                                           <  1:6459552/101‑1 (MQ=255)
                 aGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGGCCCCggc                                                                          >  2:803958/1‑100 (MQ=255)
                 aGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGa                                                                          >  1:4645725/1‑101 (MQ=255)
                 aGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGa                                                                          >  1:4080933/1‑101 (MQ=255)
                  gTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGAc                                                                         >  1:2717552/1‑101 (MQ=255)
                  gTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCCCCGGAc                                                                         <  1:3702086/101‑1 (MQ=255)
                   tGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCCCCGGAcc                                                                        <  1:7687268/101‑1 (MQ=255)
                    gACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGAccc                                                                       >  2:8408257/1‑101 (MQ=255)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGAccccc                                                                     <  1:1927837/101‑1 (MQ=38)
                       cGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGCCCCCCg                                                                    >  2:7811948/1‑101 (MQ=255)
                        gagaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCgg                                                                   <  2:2334774/101‑1 (MQ=38)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCggg                                                                  >  2:6291696/1‑101 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGat                                                                <  2:5832811/101‑1 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.tt                                                            <  2:9086628/101‑1 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.tt                                                            <  1:9052330/101‑1 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.tt                                                            <  1:5543922/101‑1 (MQ=38)
                               gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.ttt                                                           <  2:1218448/101‑1 (MQ=38)
                               gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.ttt                                                           <  2:160670/101‑1 (MQ=38)
                                    gTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGGGCTGCCCCCGGACCCCCGGGATA.TTTGTGGc                                                      >  1:7994096/1‑101 (MQ=37)
                                     tGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCa                                                     >  1:6638944/1‑101 (MQ=37)
                                     tGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCa                                                     >  1:7309460/1‑101 (MQ=37)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAg                                                   >  2:4895063/1‑101 (MQ=25)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAg                                                   <  2:8624692/101‑1 (MQ=25)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAg                                                   >  1:3602680/1‑101 (MQ=25)
                                        aGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGa                                                  <  2:7942026/101‑1 (MQ=25)
                                        aGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGa                                                  <  2:753088/101‑1 (MQ=25)
                                         gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGAt                                                 >  1:5209696/1‑101 (MQ=25)
                                         gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGAt                                                 >  1:7383867/1‑101 (MQ=25)
                                         gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGAt                                                 <  1:7185252/101‑1 (MQ=25)
                                          cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGATTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATg                                                <  1:7538072/101‑1 (MQ=25)
                                          cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATg                                                <  2:2742826/101‑1 (MQ=25)
                                           cccctgtgcGGCCAGGCATCAAGGGCCAGGCGGGTGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGa                                               <  2:4235763/93‑1 (MQ=18)
                                            cccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGaa                                              <  1:5651540/101‑1 (MQ=25)
                                              cGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAAc                                            >  2:3802271/1‑101 (MQ=25)
                                              cGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAAc                                            <  2:3279016/101‑1 (MQ=25)
                                                ttCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACag                                          <  2:3200921/101‑1 (MQ=25)
                                                ttCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACag                                          <  1:5834550/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACaga                                         <  2:4498503/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACaga                                         <  1:542916/101‑1 (MQ=25)
                                                  ccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACagag                                        <  1:5941628/101‑1 (MQ=25)
                                                  ccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACagag                                        <  1:488252/101‑1 (MQ=25)
                                                  ccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGAGGCAAGATGAAACagag                                        <  2:5369898/101‑1 (MQ=21)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGc                                       >  1:4492696/1‑101 (MQ=25)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGc                                       <  2:2082851/101‑1 (MQ=25)
                                                    gggCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCg                                      >  2:605788/1‑101 (MQ=25)
                                                    gggCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCg                                      <  1:927978/101‑1 (MQ=25)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGtt                                    <  2:5474910/101‑1 (MQ=25)
                                                       caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTa                                   <  2:6561535/101‑1 (MQ=25)
                                                        aggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.g                                 >  2:1538523/1‑101 (MQ=25)
                                                        aggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.g                                 <  1:5583615/101‑1 (MQ=25)
                                                          gcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.ggg                               <  1:1596642/101‑1 (MQ=25)
                                                          gcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.ggg                               >  2:359940/1‑101 (MQ=25)
                                                           caTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGa                              >  1:4273855/1‑101 (MQ=25)
                                                           caTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCCGCCCCCGGACCCCCCGGGGA.TTTGTGGCAAGATGAAAAAAAGGGGTA.gggg                              >  1:6265971/1‑100 (MQ=16)
                                                            aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGCTC.TTTGTGGCAAGATGAAACAGAGCGTGA.ggggt                             >  2:2974633/1‑99 (MQ=21)
                                                            aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGAt                             <  2:7016400/101‑1 (MQ=25)
                                                            aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGAt                             <  2:504090/101‑1 (MQ=25)
                                                             tCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGAtt                            >  2:3460637/1‑101 (MQ=25)
                                                             tCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGAtt                            <  2:6707487/101‑1 (MQ=25)
                                                               aaGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTg                          <  2:1066185/101‑1 (MQ=25)
                                                               aaGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.CG.GGGGGGTGGCCCCGGACCCCCCGGGTATTTTGTGGGAGAATAAAACA.AGCGGTA.gggggttg                          >  2:3047173/1‑96 (MQ=255)
                                                                aGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTgg                         <  2:6425115/101‑1 (MQ=25)
                                                                 gggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGt                        <  1:6963045/101‑1 (MQ=25)
                                                                 gggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGt                        <  2:6170330/101‑1 (MQ=25)
                                                                 gggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGt                        <  2:2895256/101‑1 (MQ=25)
                                                                  ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTa                       <  1:7066742/101‑1 (MQ=25)
                                                                  ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTa                       <  1:8759839/101‑1 (MQ=25)
                                                                    ccAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAAc                     <  2:7323115/101‑1 (MQ=25)
                                                                    ccAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAAc                     >  1:7001560/1‑101 (MQ=25)
                                                                     cAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAAcc                    <  2:3996009/101‑1 (MQ=25)
                                                                         cGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGCATTTGGTAACCATGa                <  1:8838355/101‑1 (MQ=21)
                                                                         cGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGa                >  1:1395711/1‑101 (MQ=25)
                                                                         cGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGa                >  2:7166087/1‑101 (MQ=25)
                                                                           ggggTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATc              <  1:1830365/101‑1 (MQ=25)
                                                                             ggTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCgg            <  2:8720393/101‑1 (MQ=25)
                                                                              gTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCCCCGGACCCCCCGGATA.TTTGTGGCAAGATGAAAAAGAGCGTTA.GGGGTTTGGTAACCATGATCggg           >  1:280647/1‑101 (MQ=21)
                                                                               tCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGc          >  1:3882394/1‑101 (MQ=25)
                                                                               tCGAGTTCTGATCTC.AAGGC.CGGGGGGGCGGCCCCCCGGCCCCCCGGGGA.TTTGTGGGAAGATGGAAAAGAGCGGTAGGGGGTTTGGAAACCAtatccggg           >  1:262243/1‑94 (MQ=255)
                                                                               gcGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGc          <  2:8176325/100‑1 (MQ=25)
                                                                                cGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCa         <  1:1744982/101‑1 (MQ=25)
                                                                                cGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCa         >  1:5022696/1‑101 (MQ=25)
                                                                                 gAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCaa        <  1:1172235/101‑1 (MQ=25)
                                                                                  aGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAg       >  1:6758283/1‑101 (MQ=25)
                                                                                  aGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAg       >  2:1101127/1‑101 (MQ=25)
                                                                                   gTTCTGATCTC.AAGGC.CGGGGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGc      >  1:777656/1‑101 (MQ=25)
                                                                                    ttCTGATCTC.AAGGC.CGGGGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCt     >  1:2253194/1‑101 (MQ=25)
                                                                                       tGATCTCAAGGCC.GG.GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGc  <  2:7594750/101‑1 (MQ=25)
                                                                                       tGATCTC.AAGGC.CGGGGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGc  >  2:8815834/1‑101 (MQ=25)
                                                                                                  |                                                                                             
ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C..GGGCGCTGCCCCCGGACCCCCGGGATA.TTTGTGGCAAGATGAAACAGAGCGTTA.GGGATTTGGTAACCATGATCGGGCAAGCTGGC  >  NC_014034/2209579‑2209765

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