Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,732,742 T→A 54.3% D30E (GAT→GAA yheN → predicted intracellular sulfur oxidation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,732,7420TA54.3% 11.2 / 90.9 81D30E (GAT→GAAyheNpredicted intracellular sulfur oxidation protein
Reads supporting (aligned to +/- strand):  ref base T (21/16);  new base A (32/12);  total (53/28)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

GGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAGGGGTCTATAACGCTAACCA  >  minE/2732677‑2732798
                                                                 |                                                        
gggCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGc                                                       <  1:1764119/69‑1 (MQ=255)
gggCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGc                                                       <  1:1679358/69‑1 (MQ=255)
gggCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGc                                                       <  1:678747/69‑1 (MQ=255)
gggCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGc                                                       <  1:443551/69‑1 (MQ=255)
 ggCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGc                                                       <  1:1788753/68‑1 (MQ=255)
        aTACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGa                                                  >  1:2459832/1‑66 (MQ=255)
        aTACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGa                                                  >  1:1812264/1‑66 (MQ=255)
        aTACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAg                                                 >  1:3247148/1‑67 (MQ=255)
        aTACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAg                                                 >  1:2398075/1‑67 (MQ=255)
             gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:586331/69‑1 (MQ=255)
             gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2450051/69‑1 (MQ=255)
             gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2955477/69‑1 (MQ=255)
             gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2775411/69‑1 (MQ=255)
             gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2693239/69‑1 (MQ=255)
             gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:1733794/69‑1 (MQ=255)
             gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2470664/69‑1 (MQ=255)
             gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2074620/69‑1 (MQ=255)
             gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2100599/69‑1 (MQ=255)
             gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:1147932/69‑1 (MQ=255)
             gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:1827375/69‑1 (MQ=255)
             gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2186523/69‑1 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCAt                                                    >  1:3408190/1‑50 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          >  1:96400/1‑60 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:349938/1‑69 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:2444399/1‑69 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:689858/1‑69 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:2332599/1‑69 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:847621/1‑69 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:2271545/1‑69 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:889350/1‑69 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:1128421/1‑69 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:2156328/1‑69 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:2150394/1‑69 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:1662117/1‑69 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:1056466/1‑69 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:1937723/1‑69 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttt                                   >  1:2469163/1‑67 (MQ=255)
                      aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGAAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:896153/1‑69 (MQ=255)
                         cGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:1562301/69‑1 (MQ=255)
                         cGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:669122/69‑1 (MQ=255)
                         cGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:1412962/69‑1 (MQ=255)
                         cGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:738306/69‑1 (MQ=255)
                         cGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:825245/69‑1 (MQ=255)
                         cGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:2059184/69‑1 (MQ=255)
                         cGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:2214592/69‑1 (MQ=255)
                          gAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:1738073/68‑1 (MQ=255)
                             tagtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:83268/1‑68 (MQ=255)
                             tagtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:1444800/1‑68 (MQ=255)
                             tagtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:3361876/1‑68 (MQ=255)
                               gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:896478/1‑66 (MQ=255)
                               gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:736403/1‑66 (MQ=255)
                               gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:2961415/1‑66 (MQ=255)
                               gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:604878/1‑66 (MQ=255)
                               gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:473088/1‑66 (MQ=255)
                               gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:1494815/1‑66 (MQ=255)
                               gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:3345374/1‑66 (MQ=255)
                               gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:1533597/1‑66 (MQ=255)
                                 gCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGaa                      >  1:1998609/1‑69 (MQ=255)
                                  cTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAg                     >  1:1503190/1‑69 (MQ=255)
                                  cTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAg                     >  1:3002264/1‑69 (MQ=255)
                                         gTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAGGGGTc                >  1:975685/1‑67 (MQ=255)
                                         gTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAGGGGTCt               >  1:1521455/1‑68 (MQ=255)
                                          tttGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAGGGGTCtat             <  1:2814800/69‑1 (MQ=255)
                                          tttGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTCTCTATCGGGAAGGGGTCtat             <  1:320116/69‑1 (MQ=255)
                                            tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:921910/1‑42 (MQ=255)
                                            tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:1664600/1‑42 (MQ=255)
                                            tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:209644/1‑42 (MQ=255)
                                            tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:2288623/1‑42 (MQ=255)
                                            tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:1165384/1‑42 (MQ=255)
                                            tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:1246866/1‑42 (MQ=255)
                                            tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:698278/1‑42 (MQ=255)
                                            tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:130830/1‑42 (MQ=255)
                                            tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:1926366/1‑42 (MQ=255)
                                            tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:449442/1‑42 (MQ=255)
                                            tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:406741/1‑42 (MQ=255)
                                            tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:2156090/1‑42 (MQ=255)
                                            tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:3286169/1‑42 (MQ=255)
                                            tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:2993392/1‑42 (MQ=255)
                                            tGCGCAGGCGCTAATAACAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:1999057/1‑42 (MQ=255)
                                                     gcTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAGGGGTCTATAACGCTAACCa  <  1:903156/69‑1 (MQ=255)
                                                        gATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAGGGGTCTATaa           >  1:2427476/1‑57 (MQ=255)
                                                                 |                                                        
GGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAGGGGTCTATAACGCTAACCA  >  minE/2732677‑2732798

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