Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,324 G→A 6.7% Y326Y (TAC→TAT ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,3240GA 6.7% 235.1 / 5.4 74Y326Y (TAC→TATompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (36/33);  new base A (2/3);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTA  >  minE/1439254‑1439393
                                                                      |                                                                     
cACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAag                                                                       >  1:757978/1‑71 (MQ=255)
 aCGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagt                                                                      <  1:3106493/71‑1 (MQ=255)
 aCGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAag                                                                       >  1:3129531/1‑70 (MQ=255)
 aCGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAag                                                                       <  1:173355/70‑1 (MQ=255)
   gAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtag                                                                    <  1:1143437/71‑1 (MQ=255)
   gAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtag                                                                    <  1:1455832/71‑1 (MQ=255)
    aGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtag                                                                    <  1:2364407/70‑1 (MQ=255)
        aaCTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtagg                                                                >  1:1348598/1‑70 (MQ=255)
        aaCTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtagg                                                                >  1:2236668/1‑70 (MQ=255)
             gTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCAc                                                          >  1:643447/1‑71 (MQ=255)
               tGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCAcc                                                         >  1:2809950/1‑70 (MQ=255)
                gtcgtcCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtagg                                                                >  1:2449834/1‑62 (MQ=255)
                 tcgtcCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAag                                                                       <  1:1090898/54‑1 (MQ=255)
                  catcCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCcaat                                                      <  1:2924583/68‑2 (MQ=17)
                  catcCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCcaat                                                      <  1:1463257/68‑2 (MQ=17)
                   gtcCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACcaacat                                                    >  1:973159/1‑71 (MQ=255)
                   gtcCAGCAGGTTGATTTTGTAATCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTCGCACcaatat                                                    >  1:1123461/1‑71 (MQ=18)
                   gtcCAGCAGGTTGATTTTGTAATCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTCGCACcaatat                                                    >  1:1640258/1‑71 (MQ=18)
                   gtcCAGCAGGTTGATTTTGTAATCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTCGCACcaata                                                     <  1:2530049/70‑3 (MQ=18)
                   gtcCAGCAGGTTGATTTTGTAATCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTCGCACcaata                                                     <  1:2054752/70‑3 (MQ=18)
                      cagcagGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAag                                                                       >  1:1724234/1‑49 (MQ=255)
                      cagcagGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAag                                                                       >  1:1240559/1‑49 (MQ=255)
                        gcagGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtagt                                                                   <  1:259653/51‑1 (MQ=255)
                          agGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACata                                             >  1:3505148/1‑71 (MQ=255)
                            gTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtagt                                                                   <  1:717467/47‑1 (MQ=255)
                            gTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtag                                                                    <  1:2207987/46‑1 (MQ=255)
                            gTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtag                                                                    <  1:757823/46‑1 (MQ=255)
                            gTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACcaac                                                      <  1:3488174/60‑1 (MQ=255)
                            gTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATAtt                                           <  1:3667023/71‑1 (MQ=255)
                            gTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGcc                                                         <  1:2490528/57‑1 (MQ=14)
                             ttGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtagt                                                                   <  1:2067705/46‑1 (MQ=255)
                             ttGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtagt                                                                   <  1:1291669/46‑1 (MQ=255)
                              tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagta                                                                     >  1:3075523/1‑43 (MQ=255)
                              tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagta                                                                     >  1:3934270/1‑43 (MQ=255)
                              tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAgtagta                                                                  >  1:2249534/1‑46 (MQ=255)
                              tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAgtagta                                                                  >  1:2582508/1‑46 (MQ=255)
                              tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtaggta                                                              >  1:1740684/1‑50 (MQ=255)
                              tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATAttt                                          >  1:3326456/1‑70 (MQ=255)
                              tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATAttt                                          >  1:3466632/1‑70 (MQ=255)
                               gATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAag                                                                       <  1:2574066/40‑1 (MQ=255)
                               gATTTTGTAATCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTCGCACCAATATCAACGTAttta                                         <  1:2795632/70‑2 (MQ=17)
                                aTTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAg                                       <  1:3713999/71‑1 (MQ=255)
                                   ttGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGa                                     >  1:3307598/1‑70 (MQ=255)
                                   ttGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGa                                     >  1:989524/1‑70 (MQ=255)
                                   ttGTAATCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTCGCACCAATATCAACGTAtttaac                                       >  1:3859438/1‑65 (MQ=17)
                                     gTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGat                                    <  1:2152914/69‑1 (MQ=255)
                                      tAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGatat                                  <  1:757280/70‑1 (MQ=255)
                                        gTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCt                                <  1:2311225/70‑1 (MQ=255)
                                        gTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCt                                <  1:61310/70‑1 (MQ=255)
                                        gTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCt                                <  1:1263415/70‑1 (MQ=255)
                                         tCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATAtt                                           <  1:1438435/58‑1 (MQ=255)
                                         tCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtc                              <  1:3228435/71‑1 (MQ=255)
                                         tCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTCGCACcaatatc                                                   <  1:3065161/50‑1 (MQ=11)
                                          cAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCtt                               >  1:2219639/1‑69 (MQ=255)
                                          cAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtcg                             >  1:1811161/1‑71 (MQ=255)
                                           aaCGTAGGTGGACATGTTTTTTTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtc                              >  1:3473268/1‑69 (MQ=255)
                                           aaCGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGatat                                  >  1:3354198/1‑65 (MQ=255)
                                           aaCGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtcgt                            >  1:2702276/1‑71 (MQ=255)
                                           aaCGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtcgt                            >  1:627788/1‑71 (MQ=255)
                                           aaCGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtcgt                            >  1:2967297/1‑71 (MQ=255)
                                           aaCGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCcaa                                                       >  1:1669437/1‑44 (MQ=11)
                                           aaCGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCcaa                                                       >  1:3953001/1‑44 (MQ=11)
                                             cGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGat                                    <  1:553647/61‑1 (MQ=255)
                                              gTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtcgtcgt                         <  1:2832163/71‑1 (MQ=255)
                                               tAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTc                                         <  1:3108987/54‑1 (MQ=255)
                                                 ggTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTcgtcgtcgt                      <  1:2130629/71‑1 (MQ=255)
                                                  gTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAg                                       <  1:1505036/53‑1 (MQ=255)
                                                   tGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATAtt                                           <  1:1831647/48‑1 (MQ=255)
                                                    ggACATGTTTTTGTTGAAGTAGTAGGTAGCACcaacatcaa                                                 >  1:1448028/1‑41 (MQ=255)
                                                     gACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTcgtcgtcgt                      >  1:1593688/1‑67 (MQ=255)
                                                     gACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGc                   <  1:655539/70‑1 (MQ=255)
                                                     gACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGc                   >  1:865785/1‑70 (MQ=255)
                                                      aCATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGc                   <  1:417207/69‑1 (MQ=255)
                                                       cATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCAc                >  1:1207108/1‑71 (MQ=255)
                                                       cATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCAc                <  1:118860/71‑1 (MQ=255)
                                                       cATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCAc                >  1:1799226/1‑71 (MQ=255)
                                                        aTGTTTTTGTTGAAGTAGTAGGTCGCACcaatatcaacgt                                              <  1:3877550/40‑3 (MQ=18)
                                                         tGTTTTTGTTGAAGTAGTAGGTAGCACcaacatcaaca                                               <  1:2117872/38‑1 (MQ=255)
                                                         tGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAg                                       <  1:702818/46‑1 (MQ=255)
                                                         tGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAg                                       <  1:2132664/46‑1 (MQ=255)
                                                         tGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtcgtcgtc                        <  1:2729639/61‑1 (MQ=255)
                                                          gTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAc             >  1:3050398/1‑71 (MQ=255)
                                                           tttttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAcc            >  1:2397626/1‑71 (MQ=255)
                                                           tttttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAcc            >  1:3367467/1‑71 (MQ=255)
                                                           tttttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAcc            >  1:4108720/1‑71 (MQ=255)
                                                            ttttGTTGAAGTAGTAGGTCGCACCAATATCAACGTAttta                                         <  1:4015651/41‑2 (MQ=21)
                                                            ttttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGatat                                  >  1:768473/1‑48 (MQ=255)
                                                             tttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAccc           <  1:544190/70‑1 (MQ=255)
                                                             tttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAccc           <  1:1544930/70‑1 (MQ=255)
                                                             tttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAcc            >  1:2045107/1‑69 (MQ=255)
                                                             tttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAcc            >  1:2273403/1‑69 (MQ=255)
                                                              ttgttgAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGa                                     >  1:522397/1‑43 (MQ=255)
                                                              ttgttgAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGa                                     >  1:2302862/1‑43 (MQ=255)
                                                                    aaGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGttttt   <  1:2768811/71‑1 (MQ=255)
                                                                    aaGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGtttt    <  1:912838/70‑1 (MQ=255)
                                                                     agtagtAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGttttt   <  1:132807/70‑1 (MQ=255)
                                                                     agtagtAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTa  <  1:4101981/71‑1 (MQ=255)
                                                                      gtagtaGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAccc           <  1:1277487/61‑1 (MQ=255)
                                                                      gtagtaGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAccc           <  1:3792926/61‑1 (MQ=255)
                                                                      |                                                                     
CACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTA  >  minE/1439254‑1439393

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: