Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,209,679 1 bp→GG 100% intergenic (+11/‑260) 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,6781.→G100.0% 199.7 218intergenic (+10/‑261)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (42/119):  ref base (28/29):  total (70/148)
Fisher's exact test strand distribution p-value = 2.71e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.55e-01
*NC_0140342,209,6790C→G100.0% 150.0 218intergenic (+11/‑260)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (44/119):  ref base (26/29):  total (70/148)
Fisher's exact test strand distribution p-value = 7.36e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.25e-01

Only 100 of 262 total aligned reads displayed.

ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAG.GC.C..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACAGGGA  >  NC_014034/2209579‑2209780
                                                                                                      ||||                                                                                                    
  cGAGGCGGGCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCTAGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.C..ggg                                                                                                    <  2:1315373/101‑1 (MQ=255)
  cGAGGCGGCCGGTCTAGTGACCGAGATGGGTCAGGTGAAGCCCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.C..ggg                                                                                                    <  2:3744379/101‑1 (MQ=255)
    aGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.C..ggggg                                                                                                  <  2:4603953/101‑3 (MQ=255)
     ggcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.C..GGgggc                                                                                                 <  2:2126741/101‑1 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.C..gcggggggg                                                                                              >  1:4405744/1‑93 (MQ=255)
           cGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.C..GGgggcgctgcc                                                                                           <  1:2266854/101‑7 (MQ=255)
            gATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CC.C..ggggggggggcc                                                                                           >  1:1710917/1‑92 (MQ=255)
            gATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.CGGGGGCGCTGccc                                                                                            <  1:3325693/101‑1 (MQ=255)
            gATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.C..gggggggctggcc                                                                                          >  1:2972048/1‑91 (MQ=255)
             aTCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.CGGGGGCGCTGcccc                                                                                           <  1:1363498/101‑1 (MQ=255)
             aTCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.CGGGGGCGCTGcccc                                                                                           <  1:1242589/101‑1 (MQ=255)
              tCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.CGGGGGCGCTGccccc                                                                                          <  2:3906351/101‑1 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCg                                                                                         <  1:1274500/101‑1 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCACCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCg                                                                                         <  2:4363678/101‑1 (MQ=255)
                tAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.CGGGGGCGCTGCCCCCgg                                                                                        <  2:1488735/101‑1 (MQ=255)
                 aGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CC.CG.GGGGGGTGCCCCCggg                                                                                       >  2:2626368/1‑100 (MQ=255)
                  gTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.CGGGGGCGCTGCCCCCGGAc                                                                                      <  1:1760362/101‑1 (MQ=255)
                   tGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.CGGGGGCGCTGCCCCCGGAcc                                                                                     <  2:4239858/101‑1 (MQ=255)
                    gACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CC.CG.GGGGGCTGCCCCCGGAccc                                                                                    >  1:781236/1‑101 (MQ=255)
                        gagaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCgg                                                                                <  2:4488279/101‑1 (MQ=38)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCTCCGGACCCCCggg                                                                               >  2:4663385/1‑101 (MQ=255)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCggg                                                                               >  2:4071287/1‑101 (MQ=38)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGACATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCggg                                                                               <  1:4349477/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCTCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGat                                                                             <  1:323491/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGat                                                                             <  2:2914232/101‑1 (MQ=38)
                             gggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGGGCTGCCCCCGGACCCCCGGGatat                                                                           >  1:4927205/1‑101 (MQ=38)
                             gggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGatat                                                                           <  2:781134/101‑1 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATAtt                                                                          <  1:3654839/101‑1 (MQ=38)
                              ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATAtt                                                                          >  2:787642/1‑101 (MQ=38)
                               ggcAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATAttt                                                                         <  1:2564645/99‑1 (MQ=37)
                               gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATAttt                                                                         <  2:562148/101‑1 (MQ=38)
                                 cAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGGGCCGCCCCCCGGACCCCGGGATAttttt                                                                       >  2:52516/1‑99 (MQ=255)
                                 cAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTtgt                                                                       <  1:5024313/101‑1 (MQ=37)
                                  aGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGGGCTGCCCCCGGCCCCCCCGGATATTtttg                                                                      >  1:4212080/1‑101 (MQ=255)
                                    gTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGc                                                                    <  1:286542/101‑1 (MQ=37)
                                    gTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGc                                                                    <  2:3626621/101‑1 (MQ=37)
                                     tGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCAG.G..GGGCGCTGCCCCCGGACCCCCGGGCACTTTGTggga                                                                   >  1:3256753/1‑99 (MQ=25)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAg                                                                 <  2:4747456/101‑1 (MQ=25)
                                       aaGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAg                                                                 <  2:3291002/101‑1 (MQ=25)
                                          cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCCGTCCCCGGGCCCCCCGGATATTTGCGGCaacagg                                                              >  2:1259842/1‑97 (MQ=255)
                                            cccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGaa                                                            <  2:2971820/101‑1 (MQ=25)
                                            cccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGaa                                                            <  1:3686560/101‑1 (MQ=25)
                                             ccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGCCCCCCGGGCTATTTGTGGCAAGATGaaa                                                           >  2:2518221/1‑101 (MQ=25)
                                              cGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAAc                                                          <  1:1946190/101‑1 (MQ=25)
                                              cGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAAc                                                          <  2:1138204/101‑1 (MQ=25)
                                              cGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTCGCAAGATGAAAc                                                          <  1:2737580/101‑1 (MQ=25)
                                               gTTCCGGGCAGGCATCAAGGGCCAGCCGGGTGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACa                                                         <  2:1978464/101‑1 (MQ=21)
                                                ttCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACag                                                        <  1:2074408/101‑1 (MQ=25)
                                                ttCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTACTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACag                                                        <  1:4262507/101‑1 (MQ=21)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCCGACCCCCGCGACATTTGGGGCAAGATGAAACaga                                                       >  1:2351952/1‑101 (MQ=21)
                                                 tCCGGGCAGGCATCAAGGGCCAGACGGGGGTCGAGTTCTGATCTCAAGG.AC.CG.GGGTGAAGGAACAGGACGCCCGGGATCTTTGAGGCAAGATGAAACaaa                                                       >  2:2013934/1‑99 (MQ=16)
                                                   ttgcggcgcaACAAGGGCTCGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGGTATTTGTGGCAAGATGCAACAGAGc                                                     >  2:2718362/8‑101 (MQ=12)
                                                   cGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGc                                                     <  2:3174420/101‑1 (MQ=25)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGtt                                                  <  1:3054836/101‑1 (MQ=25)
                                                       caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTa                                                 <  1:1313237/101‑1 (MQ=25)
                                                         ggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGGTATTTGTGGCAAGATGAAACAGAGCGGTAgg                                               >  2:2733639/1‑101 (MQ=21)
                                                         ggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAgg                                               <  1:4146055/101‑1 (MQ=25)
                                                            aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAt                                            >  2:411530/1‑101 (MQ=25)
                                                             tCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAtt                                           <  2:1120863/101‑1 (MQ=25)
                                                             tCAAGGGCCAGCCGGGGGACGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTCGCAAGATGAAACAGAGCGTTAGGGAtt                                           <  2:3369180/101‑1 (MQ=21)
                                                              cAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGGCCCCCGGGCTATTTGTGGGAAGATGACACAGGGGGTTAGGGAttt                                          >  1:3995447/1‑101 (MQ=25)
                                                              cAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAttt                                          <  2:1710787/101‑1 (MQ=25)
                                                               aaGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTg                                         >  2:3099047/1‑101 (MQ=25)
                                                               aaGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTg                                         >  1:1548678/1‑101 (MQ=25)
                                                                  ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTa                                      <  2:4852793/101‑1 (MQ=25)
                                                                   gCCAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTaa                                     <  2:2199251/101‑1 (MQ=25)
                                                                   gCCAGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.CGGGGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTaa                                     >  2:3039645/1‑101 (MQ=25)
                                                                    ccAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAAc                                    <  1:894772/101‑1 (MQ=25)
                                                                    ccAGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAAc                                    <  1:303715/101‑1 (MQ=25)
                                                                      aGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCa                                  <  1:374331/101‑1 (MQ=25)
                                                                      aGCCGGGGGTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCa                                  <  2:2860042/101‑1 (MQ=25)
                                                                      aGCCGGGGGTCGAGTTCTGATCTC.AAG.GC.CGGGGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCa                                  >  2:2459777/1‑101 (MQ=25)
                                                                          gggggTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGAt                              <  2:3473932/101‑1 (MQ=25)
                                                                          gggggTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGAt                              <  1:2113472/101‑1 (MQ=25)
                                                                          gggggTCGAGTTCTGATCTCAAGG.CC.CG.GGGCGCTGCCCCCGGACCCCCCGGATCTTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGGAATCATGAt                              >  1:714191/1‑101 (MQ=25)
                                                                          gggggTCGAGTTCTGATCTC.AAG.GC.CGGGGGCGCTGCCCCCCGGCCCCCGGGGTATTTTTGGCAAGATGAAACAAAGGGGTAGGGGTTTGGTAACCATGAt                              >  1:4830483/1‑101 (MQ=16)
                                                                           ggggTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATc                             <  2:189676/101‑1 (MQ=25)
                                                                           ggggTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATc                             <  1:4037804/101‑1 (MQ=25)
                                                                           ggggTCGAGTTCTGATCTC.AAG.GC.CGGGGGCGCTGCCCCCGGACCCCCCGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAAACATGATc                             >  2:4810892/1‑101 (MQ=25)
                                                                            gggTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCg                            <  2:2773050/101‑1 (MQ=25)
                                                                              gTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCggg                          <  2:870513/101‑1 (MQ=25)
                                                                              gTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCggg                          <  2:2050065/101‑1 (MQ=25)
                                                                              gTCGAGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCggg                          <  1:2424658/101‑1 (MQ=25)
                                                                                  aGTTCTGATCTCAAGG.CCGG..GGGCGCTGCCTCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGAAAg                      <  1:4663708/101‑1 (MQ=21)
                                                                                   gTTCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGc                     <  1:1631357/101‑1 (MQ=25)
                                                                                     tCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTg                   >  2:648945/1‑101 (MQ=25)
                                                                                     tCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTg                   <  2:1600119/101‑1 (MQ=25)
                                                                                     tCTGATCTCAAGG.CCGG..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGAAAGCTg                   >  1:815200/1‑101 (MQ=21)
                                                                                               tcaaggccgggggCGCTGCCCCCGGACCCCCGGGGTATTTGTGGCAAGATGAAACACAGCGTTAGGGATTTGGTACACATGATCGGGCAAGCTGGCAaaac            >  2:1913821/11‑101 (MQ=21)
                                                                                                caaggccgggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACg           <  1:3585859/92‑1 (MQ=255)
                                                                                                  aggccgggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGt         <  1:1135272/94‑1 (MQ=255)
                                                                                                  aggccgggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCCTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGt         <  2:2351788/94‑1 (MQ=255)
                                                                                                    gcctggggCTCTGCCCCCGGTCCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTAc       <  1:3389483/96‑1 (MQ=255)
                                                                                                     ccgggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACa      <  2:3688063/97‑1 (MQ=255)
                                                                                                      cgggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACag     <  1:4806201/98‑1 (MQ=255)
                                                                                                       gggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACagg    >  1:517405/3‑101 (MQ=255)
                                                                                                       gggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACagg    >  2:2396173/3‑101 (MQ=255)
                                                                                                       gggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACagg    >  2:3821961/3‑101 (MQ=255)
                                                                                                        ggggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACaggg   <  1:2135312/100‑1 (MQ=255)
                                                                                                         gggCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACaggga  >  1:2622241/1‑101 (MQ=255)
                                                                                                      ||||                                                                                                    
ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAG.GC.C..GGGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCAACACGGTACAGGGA  >  NC_014034/2209579‑2209780

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG