Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,209,678 +G 70.2% intergenic (+10/‑261) 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.→G70.2% 0.0 ‑4.6 ‑0.0326intergenic (+10/‑261)cobO/cobQ1cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand):  new base (87/142):  ref base (52/45):  total (139/187)
Fisher's exact test strand distribution p-value = 1.03e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.84e-01

Only 100 of 335 total aligned reads displayed.

ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAG..GC.C..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCA  >  NC_014034/2209579‑2209766
                                                                                                       |                                                                                             
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g                                                                                           <  2:4499857/101‑1 (MQ=255)
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g                                                                                           >  1:783061/1‑101 (MQ=255)
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g                                                                                           >  2:4304121/1‑101 (MQ=255)
aTCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGGGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGGGTTCGGGTTTC.AAG..GC.C..g                                                                                           >  2:5710569/1‑101 (MQ=255)
  cGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g..g                                                                                        >  1:4906921/1‑101 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..GC.G..g..gg                                                                                       >  1:2746270/1‑101 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g..ggg                                                                                      >  2:2411304/1‑100 (MQ=255)
   gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g..ggg                                                                                      >  1:943513/1‑100 (MQ=255)
    aGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g..gggg                                                                                     >  2:673280/1‑99 (MQ=255)
     ggcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..G..Ggggc                                                                                    >  2:3044539/1‑101 (MQ=255)
      gcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..G..Ggggcg                                                                                   >  2:1739604/1‑100 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..G..Ggggcgct                                                                                 >  2:6234477/1‑98 (MQ=255)
        ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..G..Ggggcgct                                                                                 >  1:1246710/1‑98 (MQ=255)
         ggctaTCCAGTGACCAATATCGGTCATGTAAATCC.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..G..Ggggcgctg                                                                                <  1:5695582/97‑5 (MQ=255)
         gcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..G..Ggggcgctg                                                                                <  1:5776816/101‑5 (MQ=255)
         gcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..G..Ggggcgctg                                                                                <  2:5581156/101‑5 (MQ=255)
             aTCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGcccc                                                                              >  1:5653442/1‑101 (MQ=255)
              tCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGG..GGCGCTGccccc                                                                             <  2:5763018/101‑1 (MQ=255)
              tCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..GGTGGCGCTGccccc                                                                             <  1:6309315/101‑1 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCg                                                                            >  1:6389082/1‑101 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCg                                                                            >  1:2235769/1‑101 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGG..GGCGCTGCCCCCg                                                                            <  1:1696547/101‑1 (MQ=255)
               cTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..G..Ggggcgctgccagat                                                                          >  1:6777581/1‑91 (MQ=255)
                  gTGACCGCGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGG..GGCGCTGCCCCCGGAc                                                                         <  2:6059931/101‑1 (MQ=255)
                   tGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGG..GGCGCTGCCCCCGGAcc                                                                        <  1:33577/101‑1 (MQ=255)
                      ccGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGAccccc                                                                     <  1:52967/101‑1 (MQ=38)
                       cGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCg                                                                    >  1:3086254/1‑101 (MQ=38)
                       cGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCg                                                                    <  2:4142508/101‑1 (MQ=38)
                       cGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCg                                                                    >  1:2084234/1‑101 (MQ=38)
                        gagaTGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCgg                                                                   >  2:7135954/1‑101 (MQ=38)
                          gaTGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGa                                                                 <  2:1805137/101‑1 (MQ=38)
                           aTGGGGCAGGTGAAGCG.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGat                                                                <  2:1033298/101‑1 (MQ=37)
                               gTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATAttt                                                            <  2:4873518/101‑1 (MQ=38)
                               gTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATAttt                                                            >  2:3924250/1‑101 (MQ=38)
                                  aGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG.G..G..GGGGCGGCCCCCGGCCCCCCggggatttttg                                                          >  1:2353083/1‑93 (MQ=255)
                                  aGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTtgtg                                                         >  1:6636308/1‑101 (MQ=37)
                                    gTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGGTATTTGTGGc                                                       >  1:3032725/1‑101 (MQ=37)
                                    gTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGc                                                       <  1:1347302/101‑1 (MQ=37)
                                      gAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCaa                                                     >  1:5227729/1‑101 (MQ=37)
                                      gAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCaa                                                     >  2:1105581/1‑101 (MQ=37)
                                       aaGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAg                                                    <  1:1865664/101‑1 (MQ=25)
                                        aGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGa                                                   >  2:5921954/1‑101 (MQ=25)
                                         gCA.CCCGTTCCGGGCAGGCATCAACGGCCAGCCGGGCGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAt                                                  <  2:7031303/101‑1 (MQ=21)
                                           cccccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.g                                                <  1:6947705/99‑1 (MQ=25)
                                           a.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.Ga                                               >  2:6371564/1‑101 (MQ=25)
                                           a.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.Ga                                               <  1:3397579/101‑1 (MQ=25)
                                            ccccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.Ga                                               <  2:663408/100‑1 (MQ=25)
                                             cccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.Gaa                                              <  1:1312667/101‑1 (MQ=25)
                                              ccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.Gaaa                                             <  1:3941619/101‑1 (MQ=25)
                                               cGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAAc                                            <  2:2833274/101‑1 (MQ=25)
                                                gTTCCGTGCAGGCAGCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACa                                           <  1:1508512/101‑1 (MQ=21)
                                                gTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACa                                           <  2:5496448/101‑1 (MQ=25)
                                                gTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACa                                           >  2:5270858/1‑101 (MQ=25)
                                                gTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACa                                           >  2:3984300/1‑101 (MQ=25)
                                                 ttCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACag                                          >  2:5394878/1‑101 (MQ=25)
                                                 ttCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACag                                          <  2:3177781/101‑1 (MQ=25)
                                                   ccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACagag                                        >  1:1772358/1‑101 (MQ=25)
                                                     gggCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCg                                      >  2:3813012/1‑101 (MQ=25)
                                                       gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGtt                                    <  1:5720842/101‑1 (MQ=25)
                                                        caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTa                                   >  2:3637725/1‑101 (MQ=25)
                                                         aggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAg                                  <  2:1337460/101‑1 (MQ=25)
                                                         aggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAg                                  >  2:2433990/1‑101 (MQ=25)
                                                          ggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAgg                                 >  1:4104707/1‑101 (MQ=25)
                                                          ggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAgg                                 <  2:1531788/101‑1 (MQ=25)
                                                           gcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAggg                                >  1:7073988/1‑101 (MQ=25)
                                                             aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGAt                              >  1:4970995/1‑101 (MQ=25)
                                                             aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGAt                              >  1:5954865/1‑101 (MQ=25)
                                                             aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGAt                              <  1:2535598/101‑1 (MQ=25)
                                                             aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGAt                              <  1:4374270/101‑1 (MQ=25)
                                                              tCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGAtt                             <  1:1019408/101‑1 (MQ=25)
                                                               cAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGAttt                            <  1:4614136/101‑1 (MQ=25)
                                                               cAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGAttt                            >  2:4967050/1‑101 (MQ=25)
                                                                 aGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCCGGGAATTTGTGGCAAGAA.GCAACAGAAAGTTAGGGGCTTgg                          >  2:1666101/1‑101 (MQ=16)
                                                                 aGGGCCAGCCGGAGGTCGAGTTCTGATCTC.AAG..GC.CGGG..GGCGCTGCCCCCGGACCCCCGGGGTATTTGTGGCCAGAT.GAAACAGAGCGTTAGGGGTTTgg                          >  1:5091996/1‑101 (MQ=21)
                                                                  gggCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGt                         >  1:2858792/1‑101 (MQ=25)
                                                                  gggCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGG..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGt                         >  2:4086838/1‑101 (MQ=25)
                                                                   ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTa                        <  1:4573249/101‑1 (MQ=25)
                                                                   ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTa                        <  1:673300/101‑1 (MQ=25)
                                                                   ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTa                        <  2:4408864/101‑1 (MQ=25)
                                                                   ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGCGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTa                        <  1:3716603/101‑1 (MQ=25)
                                                                     ccAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAAc                      <  1:6240316/101‑1 (MQ=25)
                                                                     ccAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAAc                      <  1:1652291/101‑1 (MQ=25)
                                                                      cAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAAcc                     <  1:7136156/101‑1 (MQ=25)
                                                                        gCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAACCAt                   <  2:2772527/101‑1 (MQ=25)
                                                                         ccGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTACCCATg                  <  2:3512081/101‑1 (MQ=21)
                                                                         ccGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAACCATg                  <  2:1971847/101‑1 (MQ=25)
                                                                         ccGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAACCATg                  <  1:219110/101‑1 (MQ=25)
                                                                         ccGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGG..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAACCATg                  >  2:264809/1‑101 (MQ=25)
                                                                          cGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGG..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAACCATGa                 >  2:4761424/1‑101 (MQ=25)
                                                                             gggTCGAGTTCTGATCTC.AAG..GC.CGGG..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAACCATGATCg              >  2:503362/1‑101 (MQ=25)
                                                                             gggTCGAGTTCTGATCTC.AAG..GC.CGGG..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGGAACCCTGATCg              >  2:2501945/1‑101 (MQ=21)
                                                                              ggTCGAGTTCTGCTCTC.AAG..GC.CGGG..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAACCATGATCgg             >  1:5170860/1‑101 (MQ=25)
                                                                                  gAGTTCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCaa         <  2:341956/101‑1 (MQ=25)
                                                                                   aGTTCTGATCTC.AAG..GC.CGGG..GGCGCTGCCCCCGGACCCCCGGGGTATTTGTGGGAAGAT.GAAACAGAGCGTTAGGGGTTTGGTAACCATGATCGGGCAAg        >  1:120592/1‑101 (MQ=25)
                                                                                   aGTTCTGATCTC.AAG..GC.CGGG..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAg        >  2:3698325/1‑101 (MQ=25)
                                                                                     ttCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCt      <  2:585837/101‑1 (MQ=25)
                                                                                     ttCTGATCTCAAGG..CCGG..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCt      >  1:825269/1‑101 (MQ=25)
                                                                                      tCTGATCTC.AAG..GC.CGGG..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTg     >  1:2017389/1‑101 (MQ=25)
                                                                                        tGATCTC.AAG..GC.CGGG..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GCAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGc   >  1:3301702/1‑101 (MQ=25)
                                                                                        tGATCTC.AAG..GC.CGGG..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGC.AGATGGAAACA.AGCGGTAAGGGTTTGGGAAACATGATCGGGCAAGCTggga  >  2:760472/1‑99 (MQ=9)
                                                                                                       |                                                                                             
ATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCA.CCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAG..GC.C..G..GGCGCTGCCCCCGGACCCCCGGGATATTTGTGGCAAGAT.GAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGCA  >  NC_014034/2209579‑2209766

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG