Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,318 G→A 6.7% N328N (AAC→AAT ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,3180GA 6.7% 335.4 / 9.7 105N328N (AAC→AATompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (47/50);  new base A (3/4);  total (51/54)
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.66e-01

Only 100 of 120 total aligned reads displayed.

CAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCC  >  minE/1439248‑1439384
                                                                      |                                                                  
cAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttg                                                                    <  1:3748841/71‑1 (MQ=255)
   cGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgttg                                                                 <  1:2411768/71‑1 (MQ=255)
    gTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgttg                                                                 >  1:4099008/1‑70 (MQ=255)
    gTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgttg                                                                 >  1:3870815/1‑70 (MQ=255)
     tCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgttg                                                                 <  1:2295035/69‑1 (MQ=255)
     tCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGaa                                                               <  1:3072708/71‑1 (MQ=255)
       aCGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgtt                                                                  >  1:3223453/1‑66 (MQ=255)
       aCGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGaa                                                               >  1:1843073/1‑69 (MQ=255)
       aCGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGaa                                                               >  1:1853380/1‑69 (MQ=255)
       aCGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGaa                                                               >  1:2192452/1‑69 (MQ=255)
       aCGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGaa                                                               >  1:2331440/1‑69 (MQ=255)
       aCGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAag                                                              >  1:3129531/1‑70 (MQ=255)
       aCGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAag                                                              <  1:173355/70‑1 (MQ=255)
         gAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtag                                                           <  1:1455832/71‑1 (MQ=255)
         gAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtag                                                           <  1:1143437/71‑1 (MQ=255)
          aGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGa                                                                <  1:2846807/65‑1 (MQ=255)
              aaCTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtagg                                                       >  1:1348598/1‑70 (MQ=255)
              aaCTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtagg                                                       >  1:2236668/1‑70 (MQ=255)
                cTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGa                                                                <  1:2689910/59‑1 (MQ=255)
                cTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGa                                                                <  1:982141/59‑1 (MQ=255)
                   gTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGa                                                                <  1:3693229/56‑1 (MQ=255)
                   gTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCAc                                                 >  1:643447/1‑71 (MQ=255)
                     tGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCAcc                                                >  1:2809950/1‑70 (MQ=255)
                      gtcgtcCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtagg                                                       >  1:2449834/1‑62 (MQ=255)
                       tcgtcCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgttg                                                                 <  1:132814/51‑1 (MQ=255)
                       tcgtcCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAag                                                              <  1:1090898/54‑1 (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)
                            cagcagGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAag                                                              >  1:1724234/1‑49 (MQ=255)
                              gcagGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtagt                                                          <  1:259653/51‑1 (MQ=255)
                                agGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACata                                    >  1:3505148/1‑71 (MQ=255)
                                agGTTGATTTTATAGTCAACGTAGGTGGACATGTTTttatt                                                                  >  1:2867153/1‑41 (MQ=18)
                                agGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGa                                                                >  1:820667/1‑43 (MQ=18)
                                  gTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgttg                                                                 <  1:3174141/40‑1 (MQ=255)
                                  gTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtagt                                                          <  1:717467/47‑1 (MQ=255)
                                  gTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtag                                                           <  1:757823/46‑1 (MQ=255)
                                  gTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtag                                                           <  1:2207987/46‑1 (MQ=255)
                                  gTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATAtt                                  <  1:3667023/71‑1 (MQ=255)
                                  gTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGcc                                                <  1:2490528/57‑1 (MQ=14)
                                   ttGATTTTGTAGTCAACGTAGGTGGACATGTTTttgttg                                                                 <  1:3935708/39‑1 (MQ=255)
                                   ttGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtagt                                                          <  1:1291669/46‑1 (MQ=255)
                                   ttGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagtagt                                                          <  1:2067705/46‑1 (MQ=255)
                                    tGATTTTGTAGTCAACGTAGGTGGACATGTTTttgt                                                                   >  1:725818/1‑36 (MQ=255)
                                    tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagta                                                            >  1:3934270/1‑43 (MQ=255)
                                    tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagta                                                            >  1:3075523/1‑43 (MQ=255)
                                    tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAgtagta                                                         >  1:2582508/1‑46 (MQ=255)
                                    tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATAttt                                 >  1:3466632/1‑70 (MQ=255)
                                    tGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATAttt                                 >  1:3326456/1‑70 (MQ=255)
                                    tGATTTTGTAATCAACATAGGTGGACATGTTTttgtt                                                                  >  1:1525115/1‑37 (MQ=14)
                                     gATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGaa                                                               >  1:706619/1‑39 (MQ=255)
                                     gATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGaa                                                               >  1:4142818/1‑39 (MQ=255)
                                      aTTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAg                              <  1:3713999/71‑1 (MQ=255)
                                           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:2967297/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)
                                                           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: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)
                                                               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)
                                                                  ttttGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGatat                         >  1:768473/1‑48 (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)
                                                                      |                                                                  
CAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCC  >  minE/1439248‑1439384

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: