Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,330 G→A 6.7% T324T (ACC→ACT ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,3300GA 6.7% 235.8 / 5.4 74T324T (ACC→ACTompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (32/37);  new base A (2/3);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

AACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTA  >  minE/1439262‑1439399
                                                                    |                                                                     
aaCTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtagg                                                                      >  1:2236668/1‑70 (MQ=255)
aaCTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtagg                                                                      >  1:1348598/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)
          catcCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCcaat                                                            <  1:1463257/68‑2 (MQ=17)
          catcCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCcaat                                                            <  1:2924583/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:3326456/1‑70 (MQ=255)
                      tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATAttt                                                >  1:3466632/1‑70 (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:1263415/70‑1 (MQ=255)
                                gTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCt                                      <  1:61310/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:627788/1‑71 (MQ=255)
                                   aaCGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtcgt                                  >  1:2967297/1‑71 (MQ=255)
                                   aaCGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtcgt                                  >  1:2702276/1‑71 (MQ=255)
                                   aaCGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCcaa                                                             >  1:3953001/1‑44 (MQ=11)
                                   aaCGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCcaa                                                             >  1:1669437/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: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:4108720/1‑71 (MQ=255)
                                                   tttttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAcc                  >  1:3367467/1‑71 (MQ=255)
                                                   tttttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAcc                  >  1:2397626/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: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:2845870/70‑1 (MQ=255)
                                                                    ggtagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTa  <  1:2053667/70‑1 (MQ=255)
                                                                    ggtagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTa  <  1:2863106/70‑1 (MQ=255)
                                                                    |                                                                     
AACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTA  >  minE/1439262‑1439399

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: