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% 123.7 171intergenic (+9/‑262)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (35/93):  ref base (23/20):  total (58/113)
Fisher's exact test strand distribution p-value = 2.69e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTG  >  NC_014034/2209579‑2209763
                                                                                                  |                                                                                          
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g                                                                                      <  1:1109642/101‑1 (MQ=255)
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..g                                                                                      <  1:3196282/101‑1 (MQ=255)
 tCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gg                                                                                     <  2:2040983/101‑1 (MQ=255)
  cGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.C..gg                                                                                     >  1:179292/1‑101 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggg                                                                                   >  2:2825330/1‑100 (MQ=255)
    aGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.C..gggg                                                                                   >  2:3043306/1‑100 (MQ=255)
    aGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..ggggg                                                                                  >  2:2598228/1‑99 (MQ=255)
       nggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gnnnngnn                                                                               >  1:617181/2‑94 (MQ=255)
       cggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGgggcgc                                                                               >  2:3430590/1‑99 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggggggct                                                                              >  2:1859007/1‑95 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGgggcgct                                                                              >  1:2700949/1‑98 (MQ=255)
         gcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..gggggggggg                                                                             >  2:2994868/1‑94 (MQ=255)
           cGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGgggcgctgcc                                                                           <  2:2242895/101‑7 (MQ=255)
           cGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGgggcgccgcc                                                                           >  1:3728886/1‑95 (MQ=255)
            gATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGc.cc                                                                           >  2:3423477/1‑101 (MQ=255)
            gATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGgggcgcgcgcc                                                                          >  2:39902/1‑94 (MQ=255)
            gATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C..GGGGGCGGc.gccc                                                                         >  2:1485945/1‑101 (MQ=255)
             aTCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGc.cca                                                                          <  2:1236533/101‑2 (MQ=255)
              tCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGc.cccc                                                                         >  1:1405025/1‑101 (MQ=255)
              tCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGc.cccc                                                                         <  1:1058664/101‑1 (MQ=255)
              tCTAGGGGCGGAGATGGGTCAGGTGAAGCCCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGc.cccc                                                                         <  2:744056/101‑1 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCg                                                                        >  2:2533287/1‑101 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGC.CCCCg                                                                        <  2:2465617/101‑1 (MQ=255)
                tAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGC.CCCCgg                                                                       <  1:1350096/101‑1 (MQ=255)
                   tGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGT.CCCCGGAcc                                                                    >  1:2641713/1‑101 (MQ=255)
                   tGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGAcc                                                                    >  1:1297425/1‑101 (MQ=255)
                   tGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGC.CCCCGGAcc                                                                    <  2:4208403/101‑1 (MQ=255)
                     aCCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGC.CCCCGGAcccc                                                                  <  1:206580/101‑1 (MQ=255)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGAccccc                                                                 <  1:3552908/101‑1 (MQ=38)
                       cGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCg                                                                >  1:1308885/1‑101 (MQ=38)
                        gagaTGGGTCAGGTGAAGCGCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCgg                                                               <  1:4086549/101‑1 (MQ=38)
                        gagaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCgg                                                               >  2:3362928/1‑101 (MQ=38)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCggg                                                              >  2:606764/1‑101 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGat                                                            <  2:1701531/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGat                                                            <  1:2600575/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGat                                                            <  1:3795052/101‑1 (MQ=38)
                            tGGGGCAAGTGAAGCCCCCGTTCCGGGCAGGCATGAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGata                                                           <  2:2690145/101‑1 (MQ=25)
                             gggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGatat                                                          <  1:712475/101‑1 (MQ=38)
                             gggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGatat                                                          <  2:3714189/101‑1 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATAtt                                                         <  2:2409563/101‑1 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATAtt                                                         >  2:412672/1‑101 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATAtt                                                         >  1:1922993/1‑101 (MQ=38)
                                 cAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTtgt                                                      >  2:450888/1‑101 (MQ=37)
                                 cAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTtgt                                                      <  2:2859474/101‑1 (MQ=37)
                                   ggTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTgg                                                    <  1:2825549/101‑1 (MQ=37)
                                   ggTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTgg                                                    <  2:1081468/101‑1 (MQ=37)
                                    gTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGc                                                   <  2:830647/101‑1 (MQ=37)
                                    gTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGc                                                   <  1:1874437/101‑1 (MQ=37)
                                     tGAAGCCCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCa                                                  <  2:4104028/101‑1 (MQ=25)
                                     tGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCa                                                  <  1:3803313/101‑1 (MQ=37)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAg                                                <  2:571198/101‑1 (MQ=25)
                                        aGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGa                                               <  1:852365/101‑1 (MQ=25)
                                        aGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGa                                               <  2:373735/101‑1 (MQ=25)
                                         gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGAt                                              <  1:3454090/101‑1 (MQ=25)
                                           aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCACGATGa                                            >  2:531327/1‑101 (MQ=25)
                                           aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGa                                            >  2:1918427/1‑101 (MQ=25)
                                             ccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGaaa                                          <  2:958562/101‑1 (MQ=25)
                                             ccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGaaa                                          <  1:3268364/101‑1 (MQ=25)
                                              cGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAAc                                         >  2:2777115/1‑101 (MQ=25)
                                               gTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACa                                        >  1:3709920/1‑101 (MQ=25)
                                               gTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGACGG..GGGGGGTGCACCCCGGCCCCCCGGATAATTTGTGGCaaaagaaaaa                                         >  1:4092268/1‑93 (MQ=255)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                      <  2:1149150/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                      <  2:1744390/101‑1 (MQ=25)
                                                    gggCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGGGCTGC.CCCCGGGCCCCCCGGATATTTGTGGGAAGATGAAACagaggg                                   >  1:3469470/1‑99 (MQ=25)
                                                    gggCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCg                                   <  2:1803450/101‑1 (MQ=25)
                                                     ggcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGGGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGt                                  >  2:1938589/1‑101 (MQ=25)
                                                       caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTa                                <  2:3919565/101‑1 (MQ=25)
                                                       caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTa                                <  2:2627715/101‑1 (MQ=25)
                                                       caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTa                                <  2:114596/101‑1 (MQ=25)
                                                        aggcaTTAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAg                               <  1:2236342/101‑1 (MQ=21)
                                                        aggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAg                               <  2:179254/101‑1 (MQ=25)
                                                        aggcaCCAACGGCCAGACAGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCCGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAg                               <  2:3806787/101‑1 (MQ=18)
                                                         ggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAgg                              <  1:1164002/101‑1 (MQ=25)
                                                          gcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCATGGCCGG..GGGCGCTGC.CCCCTGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAggg                             <  1:3613371/101‑1 (MQ=21)
                                                          gcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAggg                             <  1:2379385/101‑1 (MQ=25)
                                                            aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGGGCTGC.CCCCGGACCCCCCGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAt                           >  2:1001901/1‑101 (MQ=25)
                                                             tCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAtt                          <  1:2615286/101‑1 (MQ=25)
                                                               aaGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTg                        <  2:3008964/101‑1 (MQ=25)
                                                                 nnnnnnnnnnggggTCGNGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGt                      <  2:617117/91‑1 (MQ=18)
                                                                 gggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGt                      <  2:1297329/101‑1 (MQ=25)
                                                                 gggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGt                      <  1:450945/101‑1 (MQ=25)
                                                                  ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTa                     <  1:169210/101‑1 (MQ=25)
                                                                  ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTa                     <  1:2926399/101‑1 (MQ=25)
                                                                  ggCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGGGCCGC.CCCCCGGCCCCCCGGGGATTTGTGGCCAGAAGAAACAGAGCGGTAGGGGttttgga                     >  1:413487/1‑97 (MQ=255)
                                                                   gCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTaa                    <  2:1959642/101‑1 (MQ=25)
                                                                    ccAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAAc                   <  2:3785549/101‑1 (MQ=25)
                                                                     cAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAAcc                  <  2:3660598/101‑1 (MQ=25)
                                                                     cAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAAcc                  >  2:3168301/1‑101 (MQ=25)
                                                                       ggcGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCAt                <  1:1635365/99‑1 (MQ=25)
                                                                       gCCGGGGGTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCAt                <  1:2598419/101‑1 (MQ=25)
                                                                          gggggTCGAGTTCTGATCTC.AAGGC.CGGGGGGGCTGC.CCCCGGGCCCCCCGGGTATTTGTGGGAAGATGAAACAGAGCGTTAGGGGTTTGGGAACCATGAt             >  1:474814/1‑101 (MQ=21)
                                                                          atctgTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGAt             <  2:823121/97‑1 (MQ=21)
                                                                             ggTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCgg          <  1:2097841/101‑1 (MQ=25)
                                                                              gTCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCggg         <  1:1039616/101‑1 (MQ=25)
                                                                               tCGAGTTCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGc        <  1:3043531/101‑1 (MQ=25)
                                                                                cGAGTTCTGATCTC.AAGGC.CGGGGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCa       >  2:4170324/1‑101 (MQ=25)
                                                                                  aGTTCTGATCTC.AAGGC.CGGGGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAg     >  2:2844692/1‑101 (MQ=25)
                                                                                    ttCTGATCTCGATGCC.GG.GGGCGCTGC.CCCAGGACACCAGGGCTATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATTGGGCAAGCt   <  2:1669894/101‑1 (MQ=21)
                                                                                     tCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTg  <  1:2777334/101‑1 (MQ=25)
                                                                                     tCTGATCTCAAGGCCGG..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTg  <  1:883988/101‑1 (MQ=25)
                                                                                                  |                                                                                          
ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C..GGGCGCTGC.CCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTG  >  NC_014034/2209579‑2209763

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