Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,013,916 G→A 37.9% L30L (CTG→TTG)  exbD ← membrane spanning protein in TonB‑ExbB‑ExbD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,013,9160GA37.9% 64.3 / 70.3 111L30L (CTG→TTG) exbDmembrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base G (44/24);  new base A (29/13);  total (73/38)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 116 total aligned reads displayed.

TTCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAG  >  minE/2013867‑2013979
                                                 |                                                               
ttCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCATAACCAACATCACGTCGATa                                              <  1:1793564/69‑1 (MQ=255)
ttCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATa                                              >  1:26431/1‑69 (MQ=255)
 tCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAAAACCAACATCACGTCGATaa                                             <  1:2263678/69‑1 (MQ=255)
 tCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAAAACCAACATCACGTCGATaa                                             <  1:1784603/69‑1 (MQ=255)
 tCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAAAACCAACATCACGTCGATaa                                             <  1:865787/69‑1 (MQ=255)
 tCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAAAACCAACATCACGTCGATaa                                             <  1:1550176/69‑1 (MQ=255)
 tCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAAAACCAACATCACGTCGATaa                                             <  1:1511207/69‑1 (MQ=255)
 tCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAAAACCAACATCACGTCGATaa                                             <  1:133603/69‑1 (MQ=255)
 tCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAAAACCAACATCACGTCGATaa                                             <  1:2777696/69‑1 (MQ=255)
 tCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAAAACCAACATCACGTCGATaa                                             <  1:91713/69‑1 (MQ=255)
 tCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAAAACCAACATCACGTCGATaa                                             <  1:3028179/69‑1 (MQ=255)
 tCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAAAACCAACATCACGTCGATaa                                             <  1:3260230/69‑1 (MQ=255)
 tCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAAAACCAACATCACGTCGATaa                                             <  1:3270723/69‑1 (MQ=255)
 tCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAAAACCAACATCACGTCGATaa                                             <  1:1059838/69‑1 (MQ=255)
     aTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACgg                                         >  1:805747/1‑69 (MQ=255)
     aTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACgg                                         >  1:1071416/1‑69 (MQ=255)
       cTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGc                                        >  1:580595/1‑68 (MQ=255)
       cTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGc                                        >  1:2785482/1‑68 (MQ=255)
       cTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGc                                        >  1:803266/1‑68 (MQ=255)
       cTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGc                                        >  1:1817670/1‑68 (MQ=255)
       cTACAGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACgg                                         >  1:235644/1‑67 (MQ=255)
        tACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGt                                      >  1:3334026/1‑69 (MQ=255)
        tACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGt                                      >  1:530625/1‑69 (MQ=255)
        tACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGt                                      >  1:1126651/1‑69 (MQ=255)
        tACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGt                                      >  1:3256049/1‑69 (MQ=255)
        tACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGt                                      >  1:2786469/1‑69 (MQ=255)
        tACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGt                                      >  1:635538/1‑69 (MQ=255)
        tACCGTCGCTAACGGTGCCGCCACCAAAAAGATAATCAGCAGAACAAACATCACGTCGATAAACGGCGt                                      >  1:3261370/1‑69 (MQ=255)
               gCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGtt                                >  1:908584/1‑68 (MQ=255)
                   aCGGTGCCGCCACCATAAAGATAATCAGCAGAACCaa                                                           <  1:544696/37‑1 (MQ=255)
                   aCGGTGCCGCCACCATAAAGATAATCAGCAGAACCaa                                                           <  1:229910/37‑1 (MQ=255)
                   aCGGTGCCGCCACCATAAAGATAATCAGCAGAACCaa                                                           <  1:2673117/37‑1 (MQ=255)
                   aCGGTGCCGCCACCATAAAGATAATCAGCAGAACCaa                                                           <  1:1274283/37‑1 (MQ=255)
                   aCGGTGCCGCCACCATAAAGATAATCAGCAGAACCaa                                                           <  1:1364664/37‑1 (MQ=255)
                   aCGGTGCCGCCACCATAAAGATAATCAGCAGAACCaa                                                           <  1:927491/37‑1 (MQ=255)
                   aCGGTGCCGCCACCATAAAGATAATCAGCAGAACCaa                                                           <  1:1607560/37‑1 (MQ=255)
                     ggTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGatat                           <  1:440072/67‑1 (MQ=255)
                     ggTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGatat                           <  1:528399/67‑1 (MQ=255)
                     ggTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGatat                           <  1:1778886/67‑1 (MQ=255)
                     ggTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGatat                           <  1:2789771/67‑1 (MQ=255)
                     ggTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGatat                           <  1:891323/67‑1 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCTTTTCGCCGtt           >  1:2615357/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCGTTTCGCCGt            >  1:3107785/1‑68 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGcc              >  1:3285700/1‑66 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:2794823/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:1786783/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:862233/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:1880118/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:3115471/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:3083668/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:1144559/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:3004517/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:3419967/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:542514/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:2794145/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:1894835/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:1897024/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:1987806/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:672170/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:407209/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:2195766/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:1287022/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:1303400/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:2296897/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:2355804/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGtt           >  1:2580781/1‑69 (MQ=255)
                                   aaaGATAATCAGCAAAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGt            >  1:1875812/1‑68 (MQ=255)
                                    aaGATAATCAGCAGGACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTa          >  1:1101923/1‑69 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGatat                           <  1:371297/52‑1 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGatat                           <  1:3396807/52‑1 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGatat                           <  1:967997/52‑1 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGatat                           <  1:3309034/52‑1 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGatat                           <  1:1212414/52‑1 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGatat                           <  1:1825146/52‑1 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGatat                           <  1:1995047/52‑1 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCAt                    >  1:91342/1‑59 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTa          >  1:3033367/1‑69 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTa          >  1:1048576/1‑69 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTa          >  1:160671/1‑69 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTa          >  1:1645702/1‑69 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTa          >  1:660358/1‑69 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTa          >  1:634625/1‑69 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTa          >  1:225656/1‑69 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTa          >  1:558039/1‑69 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTa          >  1:2959829/1‑69 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTa          >  1:2651942/1‑69 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTa          >  1:2623832/1‑69 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTa          >  1:2614786/1‑69 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCAAGCATTTCGCCGTTa          >  1:2412114/1‑69 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTAACGTTGATATCATGCATTTCGCCGTTa          >  1:2522679/1‑69 (MQ=255)
                                    aaGATAATCAGCAGAACCAACATAACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTa          >  1:1947178/1‑69 (MQ=255)
                                        tAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTAtcg       >  1:2870107/1‑68 (MQ=255)
                                        tAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTAtcg       >  1:342891/1‑68 (MQ=255)
                                         aaTCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTAtcg       >  1:1665022/1‑67 (MQ=255)
                                         aaTCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTAtcg       >  1:1272791/1‑67 (MQ=255)
                                         aaTCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTAtcg       >  1:2232050/1‑67 (MQ=255)
                                            cagcagAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAg  <  1:790964/69‑1 (MQ=255)
                                            cagcagAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAg  <  1:838533/69‑1 (MQ=255)
                                            cagcagAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAg  <  1:2172198/69‑1 (MQ=255)
                                            cagcagAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAg  <  1:2741221/69‑1 (MQ=255)
                                                 |                                                               
TTCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAG  >  minE/2013867‑2013979

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: