Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,336 A→G 5.5% G322G (GGT→GGC ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,3360AG 5.5% 247.5 / 5.7 91G322G (GGT→GGCompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (39/47);  new base G (2/3);  total (41/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

GTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCA  >  minE/1439267‑1439406
                                                                     |                                                                      
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)
                 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: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:2967297/1‑71 (MQ=255)
                              aaCGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtcgt                                         >  1:2702276/1‑71 (MQ=255)
                              aaCGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtcgt                                         >  1:627788/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:1799226/1‑71 (MQ=255)
                                          cATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCAc                             >  1:1207108/1‑71 (MQ=255)
                                          cATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCAc                             <  1:118860/71‑1 (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:2397626/1‑71 (MQ=255)
                                              tttttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAcc                         >  1:4108720/1‑71 (MQ=255)
                                              tttttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAcc                         >  1:3367467/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:1277487/61‑1 (MQ=255)
                                                         gtagtaGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAccc                        <  1:3792926/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:2053667/70‑1 (MQ=255)
                                                               ggtagCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTa         <  1:1791809/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)
                                                                   gCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTa         >  1:660481/1‑66 (MQ=255)
                                                                   gCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTa         >  1:1535734/1‑66 (MQ=255)
                                                                   gCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTa         <  1:3269796/66‑1 (MQ=255)
                                                                     aCCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCa  <  1:2958476/71‑1 (MQ=255)
                                                                     |                                                                      
GTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCA  >  minE/1439267‑1439406

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: