Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,209,679 C→G 73.4% 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,6790C→G73.4% 105.1 ‑2.8 ‑0.1146intergenic (+11/‑260)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (26/81):  ref base (16/23):  total (42/104)
Fisher's exact test strand distribution p-value = 6.29e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.60e-01

Only 100 of 155 total aligned reads displayed.

ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAG..GC.C..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGC  >  NC_014034/2209579‑2209765
                                                                                                       |                                                                                         
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g                                                                                        <  2:1156663/101‑1 (MQ=255)
 ncGNNNCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..gg                                                                                       <  1:1452416/100‑1 (MQ=255)
  cGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..ggg                                                                                      >  2:2814362/1‑101 (MQ=255)
  cGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..ggg                                                                                      <  2:3634568/101‑1 (MQ=255)
  cGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..ggg                                                                                      >  2:1689248/1‑101 (MQ=255)
  cGAGGCGGCCGATATAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..ggg                                                                                      >  2:3542076/1‑101 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCAGCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..gggg                                                                                     >  2:3519761/1‑100 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CC.C..ggg                                                                                      >  1:3212029/1‑101 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..gggg                                                                                     >  2:4211116/1‑100 (MQ=255)
    aGGCGGCCTATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..ggggg                                                                                    <  2:4012140/101‑3 (MQ=255)
     gtcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..GGgggc                                                                                   <  1:2168742/99‑1 (MQ=255)
     ggcggtCCATCTGGTGGCGGACATGGCTCGGGTGAAGCGCCCGATCCGGGAAGGCGTCGAGGGCCATCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..GGgggc                                                                                   <  2:3978187/101‑1 (MQ=255)
     ggcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..GC.G..ggggg                                                                                    >  1:2307235/1‑99 (MQ=255)
       cggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..gggggggg                                                                                 >  2:1014754/1‑96 (MQ=255)
       cggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..GGgggcgc                                                                                 <  1:2231869/101‑3 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..ggggggggt                                                                                >  2:3831272/1‑95 (MQ=255)
         gcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..gggggggctg                                                                               >  1:1727672/1‑94 (MQ=255)
         gcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..GGgggcgctg                                                                               <  2:1220697/101‑5 (MQ=255)
           cGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..GGgggcgctgcc                                                                             <  2:3476819/101‑7 (MQ=255)
             aTCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.Gcccc                                                                            >  1:2729496/1‑101 (MQ=255)
              tCTAGTGGCCGAGATGGGTCAGGTGAAGCCCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGGGGCGCT.Gccccc                                                                           <  2:3209662/101‑1 (MQ=255)
              tCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGGGGCGCT.Gccccc                                                                           <  2:515466/101‑1 (MQ=255)
              tCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGGGGCGCT.Gccccc                                                                           <  2:1244895/101‑1 (MQ=255)
               cTAGTGACCGAGATTGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGGGGCGCT.GCCCCCg                                                                          <  2:1299032/101‑1 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCg                                                                          >  1:1443117/1‑101 (MQ=255)
                tAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGGGGCGCT.GCCCCCgg                                                                         <  1:1173166/101‑1 (MQ=255)
                   tGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGAcc                                                                      >  2:1031567/1‑101 (MQ=255)
                    gACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGGGCT.GCCCCCCGAccc                                                                     >  1:1024402/1‑101 (MQ=255)
                    gACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CC.CG.GGGCGCT.GCCCCCAGAccc                                                                     >  2:576334/1‑101 (MQ=255)
                    gACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGGGGCGCT.GCCCCCGGAccc                                                                     <  1:3060226/101‑1 (MQ=255)
                    gACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGGGGCGCT.GCCCCCGGAccc                                                                     <  1:1107137/101‑1 (MQ=255)
                      ccGAGATGGGTCAGGTGNAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGAccccc                                                                   <  2:83088/101‑1 (MQ=38)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGAccccc                                                                   <  1:1946486/101‑1 (MQ=38)
                      ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGAccccc                                                                   <  1:2857568/101‑1 (MQ=38)
                        gagaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCgg                                                                 >  2:2848499/1‑101 (MQ=38)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCGGGCATCAACGGCCAGCCGGGGGTCGAGTTCTGATCTCGCGG..CCGG..GGGCGCT.GCCCCCGGACCCCCggg                                                                <  1:3972133/101‑1 (MQ=25)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCggg                                                                >  2:693281/1‑101 (MQ=38)
                         agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCggg                                                                <  1:812689/101‑1 (MQ=38)
                           aggggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGat                                                              <  1:3153411/99‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGat                                                              <  2:2686132/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGat                                                              <  2:1544796/101‑1 (MQ=38)
                           aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGat                                                              <  2:1071005/101‑1 (MQ=38)
                            ggggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGata                                                             <  2:1505415/100‑1 (MQ=38)
                               gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATAttt                                                          >  1:673469/1‑101 (MQ=38)
                                   ggTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTgg                                                      >  2:1458145/1‑101 (MQ=37)
                                    gTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGc                                                     <  2:1019952/101‑1 (MQ=37)
                                        aGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGa                                                 <  1:426219/101‑1 (MQ=25)
                                         gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGAt                                                >  1:1325037/1‑101 (MQ=25)
                                          cccccGTGCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATg                                               <  2:1559291/99‑1 (MQ=25)
                                          cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATg                                               <  1:2035916/101‑1 (MQ=25)
                                          cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CC.CG.GGGGGCTGGCCCCCCAGCCCCCGGGA.ATTTTTggggaaaag                                               >  1:241764/1‑94 (MQ=255)
                                           aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGGTATTTTTGGCAAGATGa                                              >  1:376223/1‑101 (MQ=25)
                                           aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGa                                              >  1:2119024/1‑101 (MQ=25)
                                           aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGa                                              <  2:3365205/101‑1 (MQ=25)
                                           aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGa                                              <  2:44732/101‑1 (MQ=25)
                                             ccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGC.GG.G..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGaaa                                            <  1:4077201/101‑1 (MQ=25)
                                             ccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGaaa                                            <  1:1484392/101‑1 (MQ=25)
                                             ccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGaaa                                            <  2:1753727/101‑1 (MQ=25)
                                             ccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGaaa                                            <  2:206095/101‑1 (MQ=25)
                                                ttCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACag                                         <  1:3950896/101‑1 (MQ=25)
                                                 tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                        <  2:2118887/101‑1 (MQ=25)
                                                 tCCGGGAAGGCATCAAGGGCCAGCCGGTGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                        <  1:442261/101‑1 (MQ=21)
                                                  ccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACagag                                       <  1:875770/101‑1 (MQ=25)
                                                     ggcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGt                                    <  1:3551019/101‑1 (MQ=25)
                                                      gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGtt                                   >  1:2201858/1‑101 (MQ=25)
                                                       caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTa                                  <  1:2993664/101‑1 (MQ=25)
                                                         ggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAgg                                >  2:1348539/1‑101 (MQ=25)
                                                          gcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAggg                               <  1:3595288/101‑1 (MQ=25)
                                                           caTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG.G..GGGGGCG.GCCCCCGGCCCCCCGGG.TATTTTGGGGAAGATGAAAAAGAGCGTTAgggg                              >  1:2620111/1‑100 (MQ=255)
                                                            aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGGGGC.GCCCCCCGGCCCCCCGGGGATTTTTGGCCAGATGAAACAGAGCGTTAggggt                             >  2:3275667/1‑99 (MQ=255)
                                                            aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAt                             >  2:506197/1‑101 (MQ=25)
                                                            aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAt                             <  1:1009476/101‑1 (MQ=25)
                                                             tCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAtt                            >  2:2478557/1‑101 (MQ=25)
                                                             tCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAtt                            <  2:1610110/101‑1 (MQ=25)
                                                              cAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAttt                           <  1:2142802/101‑1 (MQ=25)
                                                               aaGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTg                          <  1:3542324/101‑1 (MQ=25)
                                                                 gggCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGt                        <  2:861902/101‑1 (MQ=25)
                                                                 gccccACCCGGGGGTCG.GCT.TGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTaa                      <  1:4389627/98‑1 (MQ=17)
                                                                  ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTa                       <  1:2597893/101‑1 (MQ=25)
                                                                  ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTa                       <  2:4321209/101‑1 (MQ=25)
                                                                    ccAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAAc                     <  2:3150179/101‑1 (MQ=25)
                                                                        ccGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATg                 <  2:3339823/101‑1 (MQ=25)
                                                                         cGGGGGTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGa                <  1:2785271/101‑1 (MQ=25)
                                                                          gggggTCGAGTTCTGATCTC.AAG..GC.CGGGGGCGCT.GCCCCCGGGCCCCCCGGGTATTTGTGGCAAGATGAAAAAGAGCGGTAGGGGTTTGGTAACCATGAt               >  1:398646/1‑101 (MQ=21)
                                                                           ggggTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATc              <  2:3133631/101‑1 (MQ=25)
                                                                             ggTCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCgg            <  2:2729306/101‑1 (MQ=25)
                                                                             ggTCGAGTTCTGATCTC.AAG..GC.CGGGGGCGCG.GCCCCCGGACCCCCGGGGTATTTGTGGCAAGATGAAAGAGAGCGGTAGGGATTTGGGAACCATGATCgg            >  2:4371103/1‑101 (MQ=18)
                                                                               tCGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGc          <  2:1024317/101‑1 (MQ=25)
                                                                                cGAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCa         <  2:1835853/101‑1 (MQ=25)
                                                                                 gAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCaa        <  1:3368661/101‑1 (MQ=25)
                                                                                 gAGTTCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCaa        <  2:2106309/101‑1 (MQ=25)
                                                                                 gAGTTCTGATCTC.AAG..GC.CGGGGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCaa        >  1:593092/1‑101 (MQ=25)
                                                                                 gAGTTCTGATCTC.AAG..GC.CGGGGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCaa        >  1:1316381/1‑101 (MQ=25)
                                                                                    ttCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCt     <  2:362689/101‑1 (MQ=25)
                                                                                    ttCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCt     <  1:2848706/101‑1 (MQ=25)
                                                                                    ttCTGATCTCAAGG..CCGG..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCt     <  2:593024/101‑1 (MQ=25)
                                                                                       tGATCTCAAGG..CC.GG.GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCCAGCTGGc  >  2:2062942/1‑101 (MQ=21)
                                                                                       tGATCTCAAGG..CC.GG.GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGc  >  2:4134759/1‑101 (MQ=25)
                                                                                       tGATCTCAAGG..CC.GG.GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGc  <  1:163954/101‑1 (MQ=25)
                                                                                       tGATCTCAAGG..CC.GG.GGGCGCT.GCCCCCGGAACCCCGGGGTAATTTTGGCAAGATGAAACAGAGCGTTAGGGGTTCGGGAACCATGAACGGGCAAGCTGGc  >  2:3685757/1‑101 (MQ=19)
                                                                                                       |                                                                                         
ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAG..GC.C..GGGCGCT.GCCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGC  >  NC_014034/2209579‑2209765

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