Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,002,416 A→G 11.7% intergenic (+63/+60) 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,4160AG11.7% 180.9 / 12.4 68intergenic (+63/+60)yghU/hybGpredicted S‑transferase/hydrogenase 2 accessory protein
Reads supporting (aligned to +/- strand):  ref base A (36/24);  new base G (3/5);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

AGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGG  >  Exported/2002346‑2002486
                                                                      |                                                                        
aGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                          >  1:628797/1‑71 (MQ=255)
aGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                          <  1:5217679/71‑1 (MQ=255)
          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:3300464/1‑71 (MQ=255)
                   cgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATg                                                      >  1:800646/1‑71 (MQ=255)
                    gtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGc                                                     <  1:2251105/71‑1 (MQ=255)
                       gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGAc                                                  <  1:3554549/71‑1 (MQ=255)
                       gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGAc                                                  >  1:3892702/1‑71 (MQ=255)
                       gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGAc                                                  >  1:4786456/1‑71 (MQ=255)
                       gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGAc                                                  >  1:2601440/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)
                                           gATAAGGCGTTCACGCCGCATCCGGCAATCGTGcac                                                                  <  1:1027287/36‑1 (MQ=255)
                                            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:2040050/71‑1 (MQ=255)
                                                  cGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGaa                      <  1:4972006/71‑1 (MQ=255)
                                                     tCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGc                   >  1:2818161/1‑71 (MQ=255)
                                                     tCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGc                   >  1:1811793/1‑71 (MQ=255)
                                                     tCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGc                   >  1:2566937/1‑71 (MQ=255)
                                                     tCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGc                   >  1:5073097/1‑71 (MQ=255)
                                                     tCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGc                   >  1:445070/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:1283253/71‑1 (MQ=255)
                                                            gcATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat            <  1:1862813/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:3221371/1‑71 (MQ=255)
                                                                    cAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCt    >  1:4051214/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)
                                                                      |                                                                        
AGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACA‑ATGCCTGATGCGACGC‑TGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGG  >  Exported/2002346‑2002486

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: ‑