Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,222,865 G→A 42.0% G480G (GGG→GGA yebT → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,222,8650GA42.0% 24.5 / 48.8 69G480G (GGG→GGAyebTconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/27);  new base A (14/15);  total (27/42)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

GTGAAGGTGGGGCTGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTAT  >  minE/1222800‑1222919
                                                                 |                                                      
gtgAAGGTGGGGCTGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATa                                                     >  1:3053920/1‑69 (MQ=255)
gtgAAGGTGGGGCTGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATa                                                     >  1:127789/1‑69 (MQ=255)
gtgAAGGTGGGGCTGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATa                                                     >  1:2898727/1‑69 (MQ=255)
   aaGGTGGGGCTGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACg                                                   <  1:2310408/68‑1 (MQ=255)
      gTGGGGCTGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTAtt                                               >  1:1623469/1‑69 (MQ=255)
                     gTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATaaa                                <  1:3374479/69‑1 (MQ=255)
                     gTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATaaa                                <  1:2264613/69‑1 (MQ=255)
                     gTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATaaa                                <  1:625599/69‑1 (MQ=255)
                     gTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATaaa                                <  1:704830/69‑1 (MQ=255)
                     gTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATaaa                                <  1:813976/69‑1 (MQ=255)
                           tttCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTgag                          >  1:2124943/1‑69 (MQ=255)
                           tttCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTgag                          >  1:1692588/1‑69 (MQ=255)
                           tttCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTga                           >  1:2008975/1‑68 (MQ=255)
                                  gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                    >  1:1441117/1‑68 (MQ=255)
                                  gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                    >  1:1075892/1‑68 (MQ=255)
                                  gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                    >  1:1170015/1‑68 (MQ=255)
                                  gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                    >  1:1905203/1‑68 (MQ=255)
                                  gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                    >  1:2393636/1‑68 (MQ=255)
                                  gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                    >  1:2988393/1‑68 (MQ=255)
                                  gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                    >  1:2619194/1‑68 (MQ=255)
                                  gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                    >  1:732136/1‑68 (MQ=255)
                                  gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                    >  1:2927065/1‑68 (MQ=255)
                                  gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                    >  1:804513/1‑68 (MQ=255)
                                  gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                    >  1:157121/1‑68 (MQ=255)
                                  gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                    >  1:3327919/1‑68 (MQ=255)
                                  gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                    >  1:2912283/1‑68 (MQ=255)
                                  gTGCCAGCGCCTCAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGaaa                    >  1:144889/1‑68 (MQ=255)
                                       aGCGCCTCAGAATGGATCAATGGCGGGATACGTAt                                                >  1:1938027/1‑35 (MQ=255)
                                       aGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGc              >  1:1289915/1‑69 (MQ=255)
                                       aGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGc              >  1:388387/1‑69 (MQ=255)
                                       aGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGc              >  1:2502480/1‑69 (MQ=255)
                                              cataaTGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                <  1:1060533/57‑1 (MQ=255)
                                              cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGTCGATAAAGGCGAGATGAAAGCCa                <  1:1990151/60‑1 (MQ=255)
                                              cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                <  1:2014779/60‑1 (MQ=255)
                                              cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                <  1:2992223/60‑1 (MQ=255)
                                              cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                <  1:1553681/60‑1 (MQ=255)
                                              cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                <  1:2418095/60‑1 (MQ=255)
                                              cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                <  1:1442937/60‑1 (MQ=255)
                                              cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                <  1:3319001/60‑1 (MQ=255)
                                              cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                <  1:3182979/60‑1 (MQ=255)
                                              cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                <  1:2971013/60‑1 (MQ=255)
                                              cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                <  1:1340909/60‑1 (MQ=255)
                                              cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                <  1:1307341/60‑1 (MQ=255)
                                              cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCGATAAACGCGAGATGAAAGCCa                <  1:122040/60‑1 (MQ=255)
                                              cAGAATGGATCAATGGCGGAATACGTATTCTGCCGGGCCATAAAGGCGAGATGAAAGCCa                <  1:2620653/60‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAGGCCAGCTATCCACTg     <  1:1398691/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:992538/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:3390389/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:308074/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:844736/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:587139/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:2257479/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:556369/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:2954673/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:2821099/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:2621188/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:2129709/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:2100379/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:2053302/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:1897147/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:1685291/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:1574888/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:1547393/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:1396870/68‑1 (MQ=255)
                                                 aaTGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTg     <  1:1075718/68‑1 (MQ=255)
                                                   tGGATCAATGGCGGGATACGTATTCTGCCGGGCGATaaa                                <  1:2400822/39‑1 (MQ=255)
                                                    ggATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTAt  >  1:1759572/1‑68 (MQ=255)
                                                        cAATGGCGGAATACGTATTCTGCCGGGCGATAAAGGCGAGATGAAAGCCa                <  1:2158091/50‑1 (MQ=39)
                                                           tGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCAAGCTATCCACTGTAt  >  1:1865580/1‑61 (MQ=255)
                                                                 |                                                      
GTGAAGGTGGGGCTGGATGGCGTTGAGTTTCTCGGTGCCAGCGCCTCAGAATGGATCAATGGCGGGATACGTATTCTGCCGGGCGATAAAGGTGAGATGAAAGCCAGCTATCCACTGTAT  >  minE/1222800‑1222919

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: