Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,356,770 A→G 58.1% S784S (AGT→AGC yehB ← predicted outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,356,7700AG58.1% 47.1 / 87.3 119S784S (AGT→AGCyehBpredicted outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (46/4);  new base G (25/44);  total (71/48)
Fisher's exact test for biased strand distribution p-value = 3.21e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

TCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAAT  >  minE/1356704‑1356836
                                                                  |                                                                  
tCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGAc                                                                   >  1:2489067/1‑68 (MQ=255)
               taTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:673896/67‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:3336464/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:2135437/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:219480/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:2386439/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:2696206/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:2801936/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:3397103/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:3423745/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:1337719/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:1001352/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:1206098/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:1199046/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:92319/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:905013/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:119900/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACaa                                                    <  1:437786/68‑1 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAAc                                                   >  1:2660363/1‑69 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAAc                                                   >  1:1920962/1‑69 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAAc                                                   >  1:1150498/1‑69 (MQ=255)
               aaTTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAAc                                                   >  1:2977441/1‑69 (MQ=255)
                   gCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                >  1:1196243/1‑68 (MQ=255)
                   gCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                >  1:2511311/1‑68 (MQ=255)
                   gCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                >  1:3141242/1‑68 (MQ=255)
                   gCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                >  1:2305260/1‑68 (MQ=255)
                   gCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                >  1:1528708/1‑68 (MQ=255)
                   gCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                >  1:1901040/1‑68 (MQ=255)
                   gCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                >  1:2194755/1‑68 (MQ=255)
                   gCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                >  1:1221167/1‑68 (MQ=255)
                   gCCACATTAACCGATGGCGGTAATTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGcc                                                >  1:1080600/1‑68 (MQ=255)
                              cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                    >  1:2693085/1‑69 (MQ=255)
                              cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                    >  1:3361085/1‑69 (MQ=255)
                              cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                    >  1:2956026/1‑69 (MQ=255)
                              cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                    >  1:1912687/1‑69 (MQ=255)
                              cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                    >  1:3314248/1‑69 (MQ=255)
                              cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                    >  1:2659810/1‑69 (MQ=255)
                              cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                    >  1:2638099/1‑69 (MQ=255)
                              cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                    <  1:3098399/69‑1 (MQ=255)
                              cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                    >  1:804839/1‑69 (MQ=255)
                              cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                    >  1:764061/1‑69 (MQ=255)
                              cGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGAcc                                    <  1:1840741/69‑1 (MQ=255)
                                        aCTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCa                          >  1:315718/1‑69 (MQ=255)
                                        aCTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCa                          >  1:1617332/1‑69 (MQ=255)
                                        aCTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCa                          >  1:3152147/1‑69 (MQ=255)
                                        aCTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCa                          >  1:465554/1‑69 (MQ=255)
                                        aCTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCa                          >  1:89243/1‑69 (MQ=255)
                                        aCTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCa                          >  1:1878420/1‑69 (MQ=255)
                                         cTTCATTGGTGCTAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:539308/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGAGAACACCAATATTATGACCATGGATATCa                          >  1:3073885/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAAca                                                  >  1:1774278/1‑43 (MQ=39)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:2945617/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:1097336/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:47849/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:2066708/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:423863/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:3351704/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:1392046/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:1535949/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:636197/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:2104761/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:740612/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:2599338/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:2527842/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:765475/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:850853/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:1974189/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:874709/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:222785/47‑1 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATc                           >  1:784994/1‑67 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:15306/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:623078/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:1388273/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:348941/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:1308377/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:662811/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:635818/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:1876992/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:1869890/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:2315661/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:2490050/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:2674093/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:1569791/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:2902950/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:1531107/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:3075841/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:3164143/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:3227897/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:2270864/1‑68 (MQ=255)
                                         cTTCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCAATATTATGACCATGGATATCa                          >  1:1411655/1‑68 (MQ=255)
                                          ttCATTGGTGCGAATAAATAACTGGCTTCCCTGGCCGACAACACCa                                               <  1:2472329/46‑1 (MQ=255)
                                                    cGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCAtt                        >  1:2785965/1‑59 (MQ=255)
                                                    cGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCAtt                        >  1:1160694/1‑59 (MQ=255)
                                                        taaataaCTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGatat                            >  1:3306812/1‑51 (MQ=255)
                                                        taaataaCTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGatat                            >  1:1112243/1‑51 (MQ=255)
                                                        taaataaCTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGatat                            >  1:2831606/1‑51 (MQ=255)
                                                        taaataaCTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGatat                            >  1:3255215/1‑51 (MQ=255)
                                                        taaataaCTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGatat                            >  1:408426/1‑51 (MQ=255)
                                                         aaataaCTGACTTCCCTGGCCGACAACGCCAatat                                           <  1:2099650/35‑1 (MQ=255)
                                                                tGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAAt  <  1:2463426/69‑1 (MQ=255)
                                                                  |                                                                  
TCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAACTGACTTCCCTGGCCGACAACGCCAATATTATGACCATGGATATCATTGACTTCATAACCAAACGTTAAT  >  minE/1356704‑1356836

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: