Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,333 A→C 18.2% A323A (GCT→GCG ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,3330AC18.2% 157.4 / 38.9 77A323A (GCT→GCGompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base C (5/9);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACT  >  minE/1439267‑1439403
                                                                  |                                                                      
gTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCAc                                                                    >  1:643447/1‑71 (MQ=255)
  tGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCAcc                                                                   >  1:2809950/1‑70 (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:2054752/70‑3 (MQ=18)
      gtcCAGCAGGTTGATTTTGTAATCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTCGCACcaata                                                               <  1:2530049/70‑3 (MQ=18)
             agGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACata                                                       >  1:3505148/1‑71 (MQ=255)
               gTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACcaac                                                                <  1:3488174/60‑1 (MQ=255)
               gTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATAtt                                                     <  1:3667023/71‑1 (MQ=255)
               gTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGcc                                                                   <  1:2490528/57‑1 (MQ=14)
                 tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtaggta                                                                        >  1:1740684/1‑50 (MQ=255)
                 tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATAttt                                                    >  1:3466632/1‑70 (MQ=255)
                 tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATAttt                                                    >  1:3326456/1‑70 (MQ=255)
                  gATTTTGTAATCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTCGCACCAATATCAACGTAttta                                                   <  1:2795632/70‑2 (MQ=17)
                   aTTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAg                                                 <  1:3713999/71‑1 (MQ=255)
                      ttGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGa                                               >  1:989524/1‑70 (MQ=255)
                      ttGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGa                                               >  1:3307598/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:61310/70‑1 (MQ=255)
                           gTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCt                                          <  1:1263415/70‑1 (MQ=255)
                           gTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCt                                          <  1:2311225/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:2967297/1‑71 (MQ=255)
                              aaCGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtcgt                                      >  1:627788/1‑71 (MQ=255)
                              aaCGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtcgt                                      >  1:2702276/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:865785/1‑70 (MQ=255)
                                        gACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGc                             <  1:655539/70‑1 (MQ=255)
                                         aCATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGc                             <  1:417207/69‑1 (MQ=255)
                                          cATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCAc                          <  1:118860/71‑1 (MQ=255)
                                          cATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCAc                          >  1:1207108/1‑71 (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:2132664/46‑1 (MQ=255)
                                            tGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAg                                                 <  1:702818/46‑1 (MQ=255)
                                            tGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtcgtcgtc                                  <  1:2729639/61‑1 (MQ=255)
                                             gTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAc                       >  1:3050398/1‑71 (MQ=255)
                                              tttttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAcc                      >  1:3367467/1‑71 (MQ=255)
                                              tttttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAcc                      >  1:2397626/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:1544930/70‑1 (MQ=255)
                                                tttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAccc                     <  1:544190/70‑1 (MQ=255)
                                                tttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAcc                      >  1:2273403/1‑69 (MQ=255)
                                                tttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAcc                      >  1:2045107/1‑69 (MQ=255)
                                                 ttgttgAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGa                                               >  1:2302862/1‑43 (MQ=255)
                                                 ttgttgAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGa                                               >  1:522397/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:3792926/61‑1 (MQ=255)
                                                         gtagtaGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAccc                     <  1:1277487/61‑1 (MQ=255)
                                                          tagtagGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAg                   <  1:358433/62‑1 (MQ=255)
                                                          tagtagGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTa            <  1:3307129/69‑1 (MQ=255)
                                                           agtagGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCt         >  1:1883901/1‑71 (MQ=255)
                                                            gtaggtagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGGCACGACCCAGGt                 >  1:1870377/1‑62 (MQ=255)
                                                            gtaggtagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCtt        <  1:1794312/71‑1 (MQ=255)
                                                             taggtagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCttt       <  1:932219/71‑1 (MQ=255)
                                                               ggtagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCtt        >  1:1468403/1‑68 (MQ=255)
                                                               ggtagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTa      <  1:1791809/70‑1 (MQ=255)
                                                               ggtagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTa      <  1:2053667/70‑1 (MQ=255)
                                                               ggtagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTa      <  1:2845870/70‑1 (MQ=255)
                                                               ggtagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTa      <  1:2863106/70‑1 (MQ=255)
                                                                gtagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTa      <  1:3580281/69‑1 (MQ=255)
                                                                 tagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGa    >  1:3166963/1‑70 (MQ=255)
                                                                  agCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACt  >  1:1707152/1‑71 (MQ=255)
                                                                  |                                                                      
GTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACT  >  minE/1439267‑1439403

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