Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,002,425:1 +G 9.1% intergenic (+72/+51) yghU → / ← hybG predicted S‑transferase/hydrogenase 2 accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,002,4251.G 9.1% 295.2 / 28.6 77intergenic (+72/+51)yghU/hybGpredicted S‑transferase/hydrogenase 2 accessory protein
Reads supporting (aligned to +/- strand):  ref base . (41/29);  new base G (3/4);  total (44/33)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

GTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCG  >  Exported/2002356‑2002494
                                                                      |                                                                      
gTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACa‑a                                                                       >  1:990705/1‑71 (MQ=255)
gTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACa‑a                                                                       <  1:4667310/71‑1 (MQ=255)
    gtgtTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGc                                                                    <  1:914092/70‑1 (MQ=255)
    gtgtTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGc                                                                    <  1:3137013/70‑1 (MQ=255)
     tgtTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑Ataga                                                                   >  1:4184107/1‑67 (MQ=255)
       ttCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTg                                                                 <  1:1579360/70‑1 (MQ=255)
       ttCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGa                                                                <  1:3091002/71‑1 (MQ=255)
       ttCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGa                                                                <  1:2235435/71‑1 (MQ=255)
        tcgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCATGCCGCATCCGGCAATCGTGCACA‑ATGCCTGAt                                                               >  1:2414298/1‑71 (MQ=255)
        tcgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGAt                                                               >  1:3347479/1‑71 (MQ=255)
         cgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATg                                                              >  1:3395797/1‑71 (MQ=255)
         cgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATg                                                              >  1:800646/1‑71 (MQ=255)
         cgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATg                                                              >  1:3300464/1‑71 (MQ=255)
          gtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGc                                                             <  1:2251105/71‑1 (MQ=255)
             gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGAc                                                          >  1:2601440/1‑71 (MQ=255)
             gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGAc                                                          >  1:4786456/1‑71 (MQ=255)
             gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGAc                                                          <  1:3554549/71‑1 (MQ=255)
             gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGAc                                                          >  1:3892702/1‑71 (MQ=255)
                gcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑t                                                      >  1:2718058/1‑71 (MQ=255)
                gcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑t                                                      >  1:1706087/1‑71 (MQ=255)
                 caAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑Tg                                                     >  1:549561/1‑71 (MQ=255)
                 caAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑Tg                                                     >  1:3460609/1‑71 (MQ=255)
                   aGCCATCCAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGcc                                                   >  1:300041/1‑71 (MQ=255)
                   aGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGc                                                    >  1:5213317/1‑70 (MQ=255)
                      cATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCcgcg                                                >  1:3412365/1‑71 (MQ=255)
                          taGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGC‑TAACGCGTCt                                             >  1:1747449/2‑71 (MQ=11)
                          cAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCt                                             >  1:547717/1‑70 (MQ=255)
                           aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTa                                           <  1:2424158/71‑1 (MQ=255)
                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGC‑TAACGCGTCtt                                            <  1:2369869/70‑1 (MQ=1)
                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGC‑TAACGCGTCTTa                                           >  1:5575257/1‑71 (MQ=1)
                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTAt                                          <  1:398235/71‑1 (MQ=255)
                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTAt                                          <  1:1976584/71‑1 (MQ=255)
                             gCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGC‑TAACGCGTCTTAt                                          <  1:366548/71‑1 (MQ=1)
                              ccGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCa                                        <  1:2557997/71‑1 (MQ=255)
                                 gATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGC‑TAACGCGTCTTATCAgg                                      >  1:103501/1‑71 (MQ=1)
                                  aTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGC‑TAACGCGTCTTATCAgac                                     <  1:2140569/71‑3 (MQ=1)
                                  aTAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGcc                                    >  1:299414/1‑71 (MQ=255)
                                    aaGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGc                                                        >  1:5086911/1‑50 (MQ=255)
                                    aaGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTa                                  >  1:502717/1‑71 (MQ=255)
                                     aGGCGTTCACGCCGCATCCGGCAGTCGTGCGCC‑ATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCt                                   <  1:143691/70‑1 (MQ=11)
                                     aGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCt                                                                  <  1:4704908/39‑1 (MQ=255)
                                      ggCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGC‑TAACGCGTCTTATCAGGCCTAc                                 <  1:4056090/71‑1 (MQ=1)
                                      ggCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGc                                                             <  1:4087008/43‑1 (MQ=255)
                                       gCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGa                               <  1:5425052/71‑1 (MQ=255)
                                        cGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGaa                              <  1:4972006/71‑1 (MQ=255)
                                        cGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGaa                              <  1:2040050/71‑1 (MQ=255)
                                           tCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGc                           >  1:5073097/1‑71 (MQ=255)
                                           tCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGc                           >  1:445070/1‑71 (MQ=255)
                                           tCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGc                           >  1:2566937/1‑71 (MQ=255)
                                           tCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGc                           >  1:2818161/1‑71 (MQ=255)
                                           tCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGc                           >  1:1811793/1‑71 (MQ=255)
                                            cACGCCGCATCCGGCAATCGTTCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGc                           <  1:2863666/70‑1 (MQ=255)
                                            cACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGa                               >  1:448828/1‑66 (MQ=255)
                                            cACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGc                           >  1:1488499/1‑70 (MQ=255)
                                             aCGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCt                          <  1:270106/70‑1 (MQ=255)
                                             aCGCCGCATCCGGCAATCGAGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTa                         <  1:3990030/71‑1 (MQ=255)
                                                ccgcATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTaa                        >  1:3214259/1‑69 (MQ=255)
                                                ccgcATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAatc                      <  1:5619053/71‑1 (MQ=255)
                                                  gcATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                    <  1:1862813/71‑1 (MQ=255)
                                                  gcATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                    <  1:1283253/71‑1 (MQ=255)
                                                      ccGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcg                >  1:853808/1‑71 (MQ=255)
                                                       cGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTaa                        >  1:1066588/1‑62 (MQ=255)
                                                       cGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgc               >  1:3899269/1‑71 (MQ=255)
                                                       cGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgc               >  1:2264226/1‑71 (MQ=255)
                                                         gCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACgcgc             >  1:313665/1‑71 (MQ=255)
                                                          cAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCt            >  1:4051214/1‑71 (MQ=255)
                                                          cAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCt            >  1:3221371/1‑71 (MQ=255)
                                                           aaTCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGaa                              <  1:4997144/52‑1 (MQ=255)
                                                           aaTCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTg           >  1:4606671/1‑71 (MQ=255)
                                                            aTCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTgg          <  1:1021144/71‑1 (MQ=255)
                                                             tCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAatc                      >  1:5196056/1‑58 (MQ=255)
                                                             tCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTgg          >  1:4487556/1‑70 (MQ=255)
                                                             tCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGt         <  1:5235391/71‑1 (MQ=255)
                                                              cGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                    >  1:4082724/1‑59 (MQ=255)
                                                              cGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgcg              >  1:3916511/1‑65 (MQ=255)
                                                              cGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGt         <  1:5005017/70‑1 (MQ=255)
                                                               gTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAatca                     <  1:307170/57‑1 (MQ=255)
                                                               gTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTaa       >  1:2209611/1‑71 (MQ=255)
                                                                tGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAAt      <  1:3547729/71‑1 (MQ=255)
                                                                  caca‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATg     <  1:2037908/70‑1 (MQ=255)
                                                                  caca‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATg     <  1:3390203/70‑1 (MQ=255)
                                                                  caca‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATg     <  1:3538600/70‑1 (MQ=255)
                                                                  caca‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt    <  1:2162253/71‑1 (MQ=255)
                                                                  caca‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt    <  1:887715/71‑1 (MQ=255)
                                                                  caca‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt    <  1:3471460/71‑1 (MQ=255)
                                                                  caca‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt    >  1:4828888/1‑71 (MQ=255)
                                                                   aca‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt    <  1:5081828/70‑1 (MQ=255)
                                                                    ca‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACgcgc             <  1:1787944/60‑1 (MQ=255)
                                                                    ca‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTg           >  1:4726473/1‑62 (MQ=255)
                                                                    ca‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCg  <  1:615478/71‑1 (MQ=255)
                                                                     a‑aTGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCg  >  1:3729234/1‑70 (MQ=255)
                                                                      |                                                                      
GTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCG  >  Exported/2002356‑2002494

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑