Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,209,677 G→C 74.9% 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→C74.9% 164.7 ‑4.6 0.0227intergenic (+9/‑262)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (51/119):  ref base (28/29):  total (79/148)
Fisher's exact test strand distribution p-value = 1.04e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGC  >  NC_014034/2209579‑2209765
                                                                                                  |                                                                                            
  cGAGGCGGGCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCTAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..ggg                                                                                      <  2:1315373/101‑1 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggg                                                                                     >  1:647464/1‑100 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggg                                                                                     >  2:1946053/1‑100 (MQ=255)
     ggcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggggg                                                                                   >  2:1900462/1‑98 (MQ=255)
     ggcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGgggc                                                                                   <  2:1210774/101‑1 (MQ=255)
      gcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..ggggggg                                                                                  >  2:4598732/1‑97 (MQ=255)
      gcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGgggcg                                                                                  <  1:4936555/101‑2 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGgggcgct                                                                                >  1:2352084/1‑98 (MQ=255)
             aTCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGccc.c                                                                            <  2:3765481/101‑1 (MQ=255)
             aTCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGccc.c                                                                            <  1:1363498/101‑1 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCg                                                                          <  1:1274500/101‑1 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCGGCCC.CCg                                                                          >  2:2042629/1‑101 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCC.CCg                                                                          <  1:4660440/101‑1 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCC.CCg                                                                          <  2:4459816/101‑1 (MQ=255)
                 aGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGGTGCCC.CCGGa                                                                        >  1:4818611/1‑101 (MQ=255)
                 aGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.CG.GGGGGGTGCCC.CCggg                                                                        >  2:2626368/1‑100 (MQ=255)
                  gTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCC.CCGGAc                                                                       <  1:1760362/101‑1 (MQ=255)
                  gTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCC.CCGGAc                                                                       <  2:3272398/101‑1 (MQ=255)
                   tGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCC.CCGGAcc                                                                      <  2:4239858/101‑1 (MQ=255)
                    gACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGAccc                                                                     >  2:1631225/1‑101 (MQ=255)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGAccccc                                                                   <  1:2972618/101‑1 (MQ=38)
                      ccGAGATGGGTCAAGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGAccccc                                                                   <  1:2487251/101‑1 (MQ=38)
                       cGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCg                                                                  <  1:3271101/101‑1 (MQ=38)
                        gagaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCgg                                                                 <  1:190325/101‑1 (MQ=38)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCT.CCGGACCCCCggg                                                                >  2:4663385/1‑101 (MQ=255)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCggg                                                                >  2:4071287/1‑101 (MQ=38)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCCGCCCGCCGGCCCCCCgg                                                                 >  1:1978608/1‑101 (MQ=255)
                           aTGGGTCAGGTGAAGCTCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGat                                                              <  1:323491/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGat                                                              <  2:2496510/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGat                                                              <  1:4219853/101‑1 (MQ=38)
                             gggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGatat                                                            <  1:226010/101‑1 (MQ=38)
                             gggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGatat                                                            <  2:2837495/101‑1 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATAtt                                                           <  1:3654839/101‑1 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATAtt                                                           <  1:4192438/101‑1 (MQ=38)
                               gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCATCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGGCCCCCGGGATCttt                                                          >  1:1137592/1‑101 (MQ=37)
                               gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGGGCTGCCC.CCGGACCCCCGGGAGAttt                                                          >  1:3688341/1‑101 (MQ=38)
                               gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATAttt                                                          <  1:858822/101‑1 (MQ=38)
                               gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATAttt                                                          <  2:562148/101‑1 (MQ=38)
                               cttccgatctgCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATAttt                                                          <  2:2084869/91‑1 (MQ=25)
                                tCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGGCCCCCGGGATATTtg                                                         >  1:1318567/1‑101 (MQ=38)
                                tCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTtg                                                         >  1:2305365/1‑101 (MQ=37)
                                tCAGCTGAACCCCCCTTTCCCGTCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.CG.GGGCGCTGCCC.CCGGACCCCCGGGATATTtg                                                         <  2:2892618/101‑1 (MQ=25)
                                 cAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTtgt                                                        <  2:2737931/101‑1 (MQ=37)
                                 cAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTtgt                                                        <  1:5024313/101‑1 (MQ=37)
                                     tGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTTTGGCa                                                    >  2:3836889/1‑101 (MQ=37)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAg                                                  <  2:4747456/101‑1 (MQ=25)
                                         gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATAtttagatcggaag                                                >  1:2085017/1‑91 (MQ=25)
                                         gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGAt                                                <  1:3622650/101‑1 (MQ=25)
                                           aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTGCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGa                                              <  2:347478/101‑1 (MQ=25)
                                              cGTTCCGGGCAGGCATCATGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGACGAAAc                                           <  2:3380673/101‑1 (MQ=21)
                                              cGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAAc                                           <  2:1138204/101‑1 (MQ=25)
                                               gTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGaaaaa                                          >  1:758939/1‑99 (MQ=25)
                                               gTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACa                                          <  1:4301082/101‑1 (MQ=25)
                                                ttCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACag                                         <  2:1717542/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGACGGGGGTCGAGTTCTGATCTCAAGGAC.CG.GGGTGAAGGAA.CAGGACGCCCGGGATCTTTGAGGCAAGATGAAACaaa                                        >  2:2013934/1‑99 (MQ=16)
                                                  ccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGGTATTTGTGGCAAGATGAAACagag                                       >  2:3299155/1‑101 (MQ=25)
                                                  ccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACagag                                       <  2:3003179/101‑1 (MQ=25)
                                                  ccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACagag                                       <  2:2921559/101‑1 (MQ=25)
                                                     ggcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.CG.GGGGGGCGGCC.CCGGGCCCCCGGGGGATTTGTGGGAAGAAGaaaacaagcgg                                    >  1:2116798/1‑93 (MQ=255)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGtt                                   <  2:3731651/101‑1 (MQ=25)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGtt                                   <  1:3054836/101‑1 (MQ=25)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCGGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGtt                                   <  1:4370326/101‑1 (MQ=21)
                                                       caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTa                                  <  1:4872781/101‑1 (MQ=25)
                                                        aggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAg                                 >  2:793818/1‑101 (MQ=25)
                                                        aggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAg                                 <  1:33621/101‑1 (MQ=25)
                                                         ggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.CG.GGGCGCTCCCC.CCCGACCCCCGCGCCATTTTTGGGGAGACGAAAAagagactgcgg                                >  1:1249915/1‑94 (MQ=255)
                                                           caTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGCTATTTGTGGCAAAATGAAACAGAGCGTTAGGGa                              >  2:610044/1‑101 (MQ=25)
                                                            aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAt                             >  2:411530/1‑101 (MQ=25)
                                                              cAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGGCCCCCGGGCTATTTGTGGGAAGATGACACAGGGGGTTAGGGAttt                           >  1:3995447/1‑101 (MQ=25)
                                                              cAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAANCAGAGCGTTAGGGGttt                           >  2:2223730/1‑101 (MQ=21)
                                                              cAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAttt                           <  2:715724/101‑1 (MQ=25)
                                                               aaGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTg                          <  1:2669324/101‑1 (MQ=25)
                                                               aaGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTg                          >  2:3099047/1‑101 (MQ=25)
                                                                aGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTgg                         <  2:2351920/101‑1 (MQ=25)
                                                                 gggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGt                        <  2:2662591/101‑1 (MQ=25)
                                                                  ggCGAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTa                       <  1:4775853/101‑1 (MQ=21)
                                                                  ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTa                       <  2:4852793/101‑1 (MQ=25)
                                                                  ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTa                       <  2:3627428/101‑1 (MQ=25)
                                                                    ccAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAAc                     <  1:894772/101‑1 (MQ=25)
                                                                    ccAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAAc                     <  1:2598085/101‑1 (MQ=25)
                                                                    ccAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAAc                     <  1:303715/101‑1 (MQ=25)
                                                                    ccAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAAc                     >  2:1459037/1‑101 (MQ=25)
                                                                     cAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTACcc                    <  1:2883830/101‑1 (MQ=21)
                                                                     cAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAAcc                    <  1:793920/101‑1 (MQ=25)
                                                                      aGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCa                   <  1:374331/101‑1 (MQ=25)
                                                                      aGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCCGCCC.CCCGGCCCCCCGGGTATTTTTGGCAAGACGAAACAAAGCGGTAGGGGttttggaaacc                   >  1:2854547/1‑93 (MQ=255)
                                                                      aGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGGCC.CCGGACCCCCGGGATATTTGTGGCAAGATGCAACAGAGCGCTAGGGATTTGGTAACCa                   >  2:512554/1‑101 (MQ=21)
                                                                      aGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCa                   >  2:2459777/1‑101 (MQ=25)
                                                                      aGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCC.CCGGACCCCCCCGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCa                   >  1:1521041/1‑101 (MQ=25)
                                                                       gCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCAt                  >  2:4805875/1‑101 (MQ=25)
                                                                          gggggTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGAt               <  2:2719608/101‑1 (MQ=25)
                                                                          gggggTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGAt               >  2:2135164/1‑101 (MQ=25)
                                                                           ggggTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATc              <  1:4037804/101‑1 (MQ=25)
                                                                                cGAGTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCa         <  1:167578/101‑1 (MQ=25)
                                                                                cGAGTTCTGATCTC.AAGGC.CGGGGGGGCGGCCC.CCCGGCCCCCCGGGTATTTGTGGGGAGATGAAAAAGAGGGGTAGGGGTTTGGGGACCACGATCgggga         >  1:4098077/1‑99 (MQ=255)
                                                                                 gAGTTCTGATCTC.AAGGC.CGGGGGCGCTGCCC.CCCGGCCCCCCGGGTATTTGTGGGCAGAAGAAAAAGAGCGGTAGGGGTTTTGGAACCATGATCGGGGaa        >  2:1135158/1‑101 (MQ=255)
                                                                                   gTTCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGc      <  1:1631357/101‑1 (MQ=25)
                                                                                     tCTGATCTCAAGGCCGG..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTg    <  2:1600119/101‑1 (MQ=25)
                                                                                       tGATCTCAAGGCC.GG.GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGc  >  2:783774/1‑101 (MQ=25)
                                                                                       tGATCTCAAGGCC.GG.GGGCGCTGCCC.CCCGGCCCCCCGGGTATTTTTGGGAAGATGAAACAGAGCGGTAGGGGTTTTGTAAACAAGATCGGGCAAACTggg  >  2:313868/1‑100 (MQ=12)
                                                                                                  |                                                                                            
ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C..GGGCGCTGCCC.CCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGC  >  NC_014034/2209579‑2209765

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG