Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,209,675 A→G 100% intergenic (+7/‑264) 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,6750A→G100.0% 103.2 165intergenic (+7/‑264)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (38/91):  ref base (22/13):  total (60/105)
Fisher's exact test strand distribution p-value = 6.01e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGC  >  NC_014034/2209579‑2209765
                                                                                                |                                                                                              
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...g                                                                                       <  2:4520668/101‑1 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...gggg                                                                                    >  1:986916/1‑100 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...gggg                                                                                    <  1:4298825/101‑2 (MQ=255)
    aGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...ggggg                                                                                   >  1:1414623/1‑99 (MQ=255)
     ggcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...GGgggc                                                                                  <  1:3287080/101‑1 (MQ=255)
     ggcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...GGgggc                                                                                  <  1:417646/101‑1 (MQ=255)
     ggcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...GGgggc                                                                                  <  2:3918661/101‑1 (MQ=255)
      gcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...GGgggcg                                                                                 >  2:2536246/1‑100 (MQ=255)
       gggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...GGgggcgc                                                                                <  1:3540196/100‑3 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...gggggggcc                                                                               >  1:1140512/1‑95 (MQ=255)
         gcCGATCTAGTGCCCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...GGgggccctg                                                                              <  2:3951630/101‑5 (MQ=255)
         gcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...gggggggctg                                                                              >  2:4197621/1‑94 (MQ=255)
          ccGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.C...gggggggggg                                                                              >  1:2429738/1‑94 (MQ=255)
          ccGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...GGgggcgctgc                                                                             >  2:314818/1‑96 (MQ=255)
          ccGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.C...GGgggcgctgc                                                                             >  2:4464767/1‑96 (MQ=255)
            gATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.C...gggggggggccc                                                                            >  2:63702/1‑92 (MQ=255)
                tAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGGGCTGCCCCCgg                                                                         >  2:3718537/1‑101 (MQ=255)
                 aGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGGGCTGCCCCCggc                                                                        >  1:2602060/1‑100 (MQ=255)
                 aGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGa                                                                        <  2:4653768/101‑1 (MQ=255)
                  gTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGCAc                                                                       <  1:3445745/101‑1 (MQ=255)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGAccccc                                                                   <  2:582653/101‑1 (MQ=38)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGAccccc                                                                   <  1:1525155/101‑1 (MQ=38)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGAccccc                                                                   <  1:2749778/101‑1 (MQ=38)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGCCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGAccccc                                                                   <  1:3922000/101‑1 (MQ=38)
                       cGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCg                                                                  <  1:2729026/101‑1 (MQ=38)
                       cGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCg                                                                  <  1:1003591/101‑1 (MQ=38)
                        gagaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGGGCTGCCCCCCGCCCCCCgg                                                                 >  2:1952760/1‑101 (MQ=255)
                        gagaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCgg                                                                 <  1:2718737/101‑1 (MQ=38)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCggg                                                                <  1:1398446/101‑1 (MQ=38)
                          gaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGa                                                               <  2:2536986/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGat                                                              <  2:3459915/101‑1 (MQ=38)
                             gggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGatat                                                            <  1:2573253/101‑1 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATAtt                                                           >  2:1083516/1‑101 (MQ=38)
                               gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATAttt                                                          <  2:917687/101‑1 (MQ=38)
                                tCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGC...GGGCGCTGCCCCCGGACCCCCCGGTTGtttt                                                         >  2:800385/1‑100 (MQ=255)
                                 cAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTtgt                                                        <  2:2713770/101‑1 (MQ=37)
                                 cAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTtgt                                                        >  1:4403621/1‑101 (MQ=37)
                                  cggTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTtgtg                                                       <  2:1188953/100‑1 (MQ=37)
                                    gTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCACCAGACCTCCGGGATATTTGTGGc                                                     >  2:262681/1‑101 (MQ=255)
                                      gAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGGGCTGCCCCCGGACCCCCGGGGTATTTGTGGCaa                                                   >  1:2830833/1‑101 (MQ=37)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAg                                                  <  2:3996505/101‑1 (MQ=25)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAg                                                  <  2:3254129/101‑1 (MQ=25)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAg                                                  <  1:3846514/101‑1 (MQ=25)
                                        aGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGa                                                 <  2:1155648/101‑1 (MQ=25)
                                          ctcccGTTCCGGGCAGGCATCAAGGGCCAGTCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATg                                               <  2:3864780/99‑1 (MQ=25)
                                          cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATg                                               <  1:539332/101‑1 (MQ=25)
                                           aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGa                                              <  1:1452081/101‑1 (MQ=25)
                                            cccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCCGCCCCCGGACCCCCGGGATATTTGTGGCAAGACGaa                                             >  1:630074/1‑101 (MQ=25)
                                             ccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGaaa                                            <  1:3029147/101‑1 (MQ=25)
                                              cGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAAc                                           <  1:3815757/101‑1 (MQ=25)
                                               gTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACa                                          <  1:1705176/101‑1 (MQ=25)
                                               gTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACa                                          <  1:2727677/101‑1 (MQ=25)
                                                ttCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACag                                         <  2:3495245/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                        <  2:2417470/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                        <  2:2104184/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                        >  2:1950633/1‑101 (MQ=25)
                                                  ccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGGGGTGCCCCCGGGCCCCCGGGGTATTTTTGGCAAGATGAAACagag                                       >  1:403680/1‑101 (MQ=255)
                                                  ccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCC.CG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCCAGATGAAACagag                                       >  1:2243138/1‑101 (MQ=37)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGc                                      <  1:2279563/101‑1 (MQ=25)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGc                                      <  1:1274824/101‑1 (MQ=25)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGc                                      <  2:1581226/101‑1 (MQ=25)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGCTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGc                                      <  1:3997815/101‑1 (MQ=25)
                                                    gggCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCg                                     <  1:4114484/101‑1 (MQ=25)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGtt                                   <  2:4115067/101‑1 (MQ=25)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGtt                                   <  2:3830781/101‑1 (MQ=25)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGtt                                   <  1:267778/101‑1 (MQ=25)
                                                       caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTa                                  <  1:3409107/101‑1 (MQ=25)
                                                       caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTa                                  <  1:1811178/101‑1 (MQ=25)
                                                       caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCGCAAGGCCGG...GGGCGCTGCCCCCCCACCCCCGGGATATTTGTGGCCAGATGAAACAGAGcccta                                  >  2:3830424/1‑97 (MQ=18)
                                                         ggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAgg                                >  1:2292590/1‑101 (MQ=25)
                                                           caTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGa                              >  2:590998/1‑101 (MQ=25)
                                                           caTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGa                              <  2:1385539/101‑1 (MQ=25)
                                                                 gggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGt                        <  2:646655/101‑1 (MQ=25)
                                                                 gggCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGt                        <  1:2200406/101‑1 (MQ=25)
                                                                   gCCAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGGGGGGGGCGCCCCCCGGCCCCCCGGGGAATTT.TGGCGAGAATAAAAAAAGCGTTACGGAGTTgggga                      >  1:3754447/1‑98 (MQ=255)
                                                                     cAGCCGGGGGTCGAGTTCTGATCTC.AAGGC.CGG.GGGCGCTGCCCCCGGACCCCCCGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTTGTAAcc                    >  1:778155/1‑101 (MQ=25)
                                                                      aGCCGGGGGTCGAGTTCTGATCGCACGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAGCCa                   <  2:1481456/101‑1 (MQ=21)
                                                                       gCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCAt                  <  1:2662954/101‑1 (MQ=25)
                                                                       gCCGGGGGTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCAt                  <  1:3351375/101‑1 (MQ=25)
                                                                         cGGGGGTCGAGTTCTGATCTC.AAGGC.CGG.GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGa                >  2:2834398/1‑101 (MQ=25)
                                                                          gggggTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGAt               <  1:1971111/101‑1 (MQ=25)
                                                                          gggggTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGAt               <  1:4260316/101‑1 (MQ=25)
                                                                          gggggTCGAGTTCTGATCTCAAGGCC.G...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCACGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATc              >  2:1842038/1‑101 (MQ=37)
                                                                           ggtgTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATc              <  2:3095023/98‑1 (MQ=21)
                                                                           ggggTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATc              <  2:2709089/101‑1 (MQ=25)
                                                                           ggggTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATc              <  2:4277392/101‑1 (MQ=25)
                                                                            gggTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCg             <  1:916656/101‑1 (MQ=25)
                                                                            gggTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCg             <  1:1905910/101‑1 (MQ=25)
                                                                              gTCGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCggg           <  1:4211349/101‑1 (MQ=25)
                                                                               tCGAGTTCTGATCTCAACGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGc          <  2:4190896/101‑1 (MQ=25)
                                                                                cGAGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCCGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGGAACCATGATCGGGCa         >  1:3062235/1‑101 (MQ=21)
                                                                                  aGTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAg       <  1:2851799/101‑1 (MQ=25)
                                                                                   gTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGc      <  1:591059/101‑1 (MQ=25)
                                                                                   gTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGc      <  1:1627159/101‑1 (MQ=25)
                                                                                   gTTCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGc      <  2:2068650/101‑1 (MQ=25)
                                                                                   gTTCTGATCTC.AAGGC.CGG.GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGc      >  2:1538136/1‑101 (MQ=25)
                                                                                    ttCTGATCTCAAGGCCGG...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCt     <  1:3308285/101‑1 (MQ=25)
                                                                                     tCTGATCTC.AAGGC.CGG.GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTg    >  2:4639175/1‑101 (MQ=25)
                                                                                     tCTGATCTC.AAGGC.CGG.GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTg    >  2:2252888/1‑101 (MQ=25)
                                                                                       tGATCTCAAGGCC.GG..GGGCGCTGCCCCCGGGCCCCCGGGATATTTGTGGCAAGATTAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAACTgga  >  2:3643447/1‑100 (MQ=21)
                                                                                                |                                                                                              
ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAGGC.C...GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGC  >  NC_014034/2209579‑2209765

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