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

GTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCG  >  Exported/2002356‑2002494
                                                                      |                                                                      
gTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACaa                                                                       >  1:990705/1‑71 (MQ=255)
gTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACaa                                                                       <  1:4667310/71‑1 (MQ=255)
    gtgtTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGc                                                                    <  1:914092/70‑1 (MQ=255)
    gtgtTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGc                                                                    <  1:3137013/70‑1 (MQ=255)
     tgtTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAAtaga                                                                   >  1:4184107/1‑67 (MQ=255)
       ttCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTg                                                                 <  1:1579360/70‑1 (MQ=255)
       ttCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGa                                                                <  1:3091002/71‑1 (MQ=255)
       ttCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGa                                                                <  1:2235435/71‑1 (MQ=255)
        tcgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCATGCCGCATCCGGCAATCGTGCACAATGCCTGAt                                                               >  1:2414298/1‑71 (MQ=255)
        tcgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGAt                                                               >  1:3347479/1‑71 (MQ=255)
         cgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATg                                                              >  1:3395797/1‑71 (MQ=255)
         cgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATg                                                              >  1:800646/1‑71 (MQ=255)
         cgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATg                                                              >  1:3300464/1‑71 (MQ=255)
          gtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGc                                                             <  1:2251105/71‑1 (MQ=255)
             gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGAc                                                          >  1:2601440/1‑71 (MQ=255)
             gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGAc                                                          >  1:4786456/1‑71 (MQ=255)
             gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGAc                                                          <  1:3554549/71‑1 (MQ=255)
             gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGAc                                                          >  1:3892702/1‑71 (MQ=255)
                gcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCt                                                      >  1:2718058/1‑71 (MQ=255)
                gcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCt                                                      >  1:1706087/1‑71 (MQ=255)
                 caAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTg                                                     >  1:549561/1‑71 (MQ=255)
                 caAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTg                                                     >  1:3460609/1‑71 (MQ=255)
                   aGCCATCCAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGcc                                                   >  1:300041/1‑71 (MQ=255)
                   aGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGc                                                    >  1:5213317/1‑70 (MQ=255)
                      cATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCcgcg                                                >  1:3412365/1‑71 (MQ=255)
                          taGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCt                                             >  1:1747449/2‑71 (MQ=11)
                          cAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCt                                             >  1:547717/1‑70 (MQ=255)
                           aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTa                                           <  1:2424158/71‑1 (MQ=255)
                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCtt                                            <  1:2369869/70‑1 (MQ=1)
                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTa                                           >  1:5575257/1‑71 (MQ=1)
                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTAt                                          <  1:398235/71‑1 (MQ=255)
                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTAt                                          <  1:1976584/71‑1 (MQ=255)
                             gCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTAt                                          <  1:366548/71‑1 (MQ=1)
                              ccGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCa                                        <  1:2557997/71‑1 (MQ=255)
                                 gATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAgg                                      >  1:103501/1‑71 (MQ=1)
                                  aTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAgac                                     <  1:2140569/71‑3 (MQ=1)
                                  aTAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGcc                                    >  1:299414/1‑71 (MQ=255)
                                    aaGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGc                                                        >  1:5086911/1‑50 (MQ=255)
                                    aaGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTa                                  >  1:502717/1‑71 (MQ=255)
                                     aGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCt                                   <  1:143691/70‑1 (MQ=11)
                                     aGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCt                                                                  <  1:4704908/39‑1 (MQ=255)
                                      ggCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTAc                                 <  1:4056090/71‑1 (MQ=1)
                                      ggCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGc                                                             <  1:4087008/43‑1 (MQ=255)
                                       gCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGa                               <  1:5425052/71‑1 (MQ=255)
                                        cGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGaa                              <  1:4972006/71‑1 (MQ=255)
                                        cGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGaa                              <  1:2040050/71‑1 (MQ=255)
                                           tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                           >  1:5073097/1‑71 (MQ=255)
                                           tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                           >  1:445070/1‑71 (MQ=255)
                                           tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                           >  1:2566937/1‑71 (MQ=255)
                                           tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                           >  1:2818161/1‑71 (MQ=255)
                                           tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                           >  1:1811793/1‑71 (MQ=255)
                                            cACGCCGCATCCGGCAATCGTTCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                           <  1:2863666/70‑1 (MQ=255)
                                            cACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGa                               >  1:448828/1‑66 (MQ=255)
                                            cACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                           >  1:1488499/1‑70 (MQ=255)
                                             aCGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCt                          <  1:270106/70‑1 (MQ=255)
                                             aCGCCGCATCCGGCAATCGAGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTa                         <  1:3990030/71‑1 (MQ=255)
                                                ccgcATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTaa                        >  1:3214259/1‑69 (MQ=255)
                                                ccgcATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatc                      <  1:5619053/71‑1 (MQ=255)
                                                  gcATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                    <  1:1862813/71‑1 (MQ=255)
                                                  gcATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                    <  1:1283253/71‑1 (MQ=255)
                                                      ccGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcg                >  1:853808/1‑71 (MQ=255)
                                                       cGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTaa                        >  1:1066588/1‑62 (MQ=255)
                                                       cGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgc               >  1:3899269/1‑71 (MQ=255)
                                                       cGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgc               >  1:2264226/1‑71 (MQ=255)
                                                         gCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACgcgc             >  1:313665/1‑71 (MQ=255)
                                                          cAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCt            >  1:4051214/1‑71 (MQ=255)
                                                          cAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCt            >  1:3221371/1‑71 (MQ=255)
                                                           aaTCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGaa                              <  1:4997144/52‑1 (MQ=255)
                                                           aaTCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTg           >  1:4606671/1‑71 (MQ=255)
                                                            aTCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTgg          <  1:1021144/71‑1 (MQ=255)
                                                             tCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatc                      >  1:5196056/1‑58 (MQ=255)
                                                             tCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTgg          >  1:4487556/1‑70 (MQ=255)
                                                             tCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGt         <  1:5235391/71‑1 (MQ=255)
                                                              cGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                    >  1:4082724/1‑59 (MQ=255)
                                                              cGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgcg              >  1:3916511/1‑65 (MQ=255)
                                                              cGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGt         <  1:5005017/70‑1 (MQ=255)
                                                               gTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatca                     <  1:307170/57‑1 (MQ=255)
                                                               gTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTaa       >  1:2209611/1‑71 (MQ=255)
                                                                tGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAAt      <  1:3547729/71‑1 (MQ=255)
                                                                  cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATg     <  1:2037908/70‑1 (MQ=255)
                                                                  cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATg     <  1:3390203/70‑1 (MQ=255)
                                                                  cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATg     <  1:3538600/70‑1 (MQ=255)
                                                                  cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt    <  1:2162253/71‑1 (MQ=255)
                                                                  cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt    <  1:887715/71‑1 (MQ=255)
                                                                  cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt    <  1:3471460/71‑1 (MQ=255)
                                                                  cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt    >  1:4828888/1‑71 (MQ=255)
                                                                   acaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt    <  1:5081828/70‑1 (MQ=255)
                                                                    caATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACgcgc             <  1:1787944/60‑1 (MQ=255)
                                                                    caATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTg           >  1:4726473/1‑62 (MQ=255)
                                                                    caATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCg  <  1:615478/71‑1 (MQ=255)
                                                                     aaTGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCg  >  1:3729234/1‑70 (MQ=255)
                                                                      |                                                                      
GTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCG  >  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: