Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,816 G→A 9.5% D162D (GAC→GAT) ‡ ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,8160GA 9.5% 364.8 / 20.1 115D162D (GAC→GAT) ‡ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (60/44);  new base A (5/6);  total (65/50)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 115 total aligned reads displayed.

ACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTG  >  minE/1439750‑1439885
                                                                  |                                                                     
aCTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTc                                                                     <  1:3982698/69‑1 (MQ=255)
aCTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                                   <  1:2590613/71‑1 (MQ=255)
 cTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                                                  <  1:1435957/71‑1 (MQ=255)
  tAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                                                  >  1:3863923/1‑70 (MQ=255)
       aaGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGa                                                              >  1:814612/1‑69 (MQ=255)
       aaGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGAcc                                                            >  1:2521302/1‑71 (MQ=255)
          ccTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAcc                                                                 >  1:3573862/1‑63 (MQ=255)
          ccTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAcc                                                                 >  1:2042108/1‑63 (MQ=255)
          ccTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                                         <  1:4038060/71‑1 (MQ=255)
          ccTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                                         <  1:939048/71‑1 (MQ=255)
              cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                                   >  1:3797908/1‑57 (MQ=255)
              cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                                   >  1:1057791/1‑57 (MQ=255)
              cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                                     >  1:4096180/1‑71 (MQ=255)
               aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                                    >  1:2897085/1‑71 (MQ=255)
               aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                                    >  1:2956021/1‑71 (MQ=255)
                  aGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                                                 >  1:375926/1‑71 (MQ=255)
                  aGATGGGTTGCCGTTTTTACCCTGGTACTGAACAAAAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                                                 <  1:1163407/71‑1 (MQ=255)
                     tGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                                              >  1:856898/1‑71 (MQ=255)
                      gggTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                                              >  1:1224135/1‑70 (MQ=255)
                       ggTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGtt                                             <  1:1277748/70‑1 (MQ=255)
                        gTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                                                 <  1:3481406/65‑1 (MQ=255)
                        gTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTAc                                           >  1:2736074/1‑71 (MQ=255)
                         ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCa                                                                <  1:956409/49‑1 (MQ=255)
                          cttcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACg                                          <  1:56719/67‑1 (MQ=2)
                           gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                         <  1:3717165/70‑1 (MQ=255)
                           gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                         >  1:3802273/1‑70 (MQ=255)
                           gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                         <  1:1424771/70‑1 (MQ=255)
                           gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                         >  1:829369/1‑70 (MQ=255)
                              gTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaagaa                                                      >  1:3198704/1‑54 (MQ=255)
                              gTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaagaa                                                      >  1:3146695/1‑54 (MQ=255)
                               tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                                                  <  1:1887122/41‑1 (MQ=255)
                               tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                                            <  1:1841386/63‑1 (MQ=255)
                               tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                         >  1:510827/1‑66 (MQ=255)
                               tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt                                    >  1:3220468/1‑71 (MQ=255)
                                ttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgt                                                >  1:2695949/1‑58 (MQ=255)
                                 tttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTc                                   >  1:214731/1‑70 (MQ=255)
                                  ttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc                                 <  1:2411618/71‑1 (MQ=9)
                                  ttACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgc                                 <  1:1804202/71‑1 (MQ=255)
                                  ttACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgc                                 <  1:3110226/71‑1 (MQ=255)
                                  ttACACTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgc                                 >  1:1130816/1‑71 (MQ=255)
                                   tACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgc                                 >  1:1937053/1‑70 (MQ=255)
                                   tACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgc                                 >  1:3849331/1‑70 (MQ=255)
                                     cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGaa                              <  1:3749021/71‑1 (MQ=255)
                                      ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg                             <  1:3279985/71‑1 (MQ=255)
                                       cTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                                            >  1:1374860/1‑55 (MQ=255)
                                       cTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                                            >  1:2445327/1‑55 (MQ=255)
                                       cTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGc                            <  1:2598011/71‑1 (MQ=255)
                                        tgatACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGt                                    >  1:3415493/4‑62 (MQ=2)
                                        tGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                                    >  1:3023380/1‑46 (MQ=255)
                                         ggTACTGAACAGCAAAGTTCAGGCCGTCAACCAGaa                                                             >  1:110833/1‑35 (MQ=255)
                                         gatACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGa                               <  1:695181/64‑1 (MQ=255)
                                         gatACTGCACAGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGa                               >  1:3624457/3‑66 (MQ=9)
                                          gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                         >  1:3478811/1‑55 (MQ=255)
                                          gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                         >  1:2716471/1‑55 (MQ=255)
                                          gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGGCGCGAAGcc                           <  1:2959033/69‑1 (MQ=255)
                                          gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGAGTTACGGTAGGTCGCGAAGCCGt                         >  1:1755088/1‑71 (MQ=255)
                                           tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                         <  1:835113/70‑1 (MQ=255)
                                            aCTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgt                                                >  1:3979937/1‑46 (MQ=255)
                                            aCTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                                                 >  1:1692942/1‑45 (MQ=255)
                                             cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                      >  1:3807269/1‑71 (MQ=255)
                                             cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                      <  1:1489981/71‑1 (MQ=255)
                                             cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                      <  1:3296023/71‑1 (MQ=255)
                                             cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                      <  1:3002901/71‑1 (MQ=255)
                                             cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                      <  1:2004825/71‑1 (MQ=255)
                                             cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                      >  1:3251623/1‑71 (MQ=255)
                                              tGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaat                                                        <  1:3753704/36‑2 (MQ=255)
                                              tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                                                 <  1:1233672/43‑1 (MQ=255)
                                              tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                         <  1:26319/51‑1 (MQ=255)
                                              tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                         <  1:2347834/51‑1 (MQ=255)
                                              tGAACAGCAAAGTTCAGCCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                         <  1:2017753/51‑1 (MQ=255)
                                               gAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                         <  1:1128244/50‑1 (MQ=255)
                                               gAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAg                                      >  1:391919/1‑53 (MQ=255)
                                               gAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGc                            >  1:539170/1‑63 (MQ=255)
                                                aaCAGCAAAGTTCAGGCCGTCACCCAGACCGAAGAAGTCAgtg                                               >  1:1301783/1‑43 (MQ=255)
                                                aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtg                                               >  1:33904/1‑43 (MQ=255)
                                                aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                                            >  1:896892/1‑46 (MQ=255)
                                                aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcg                                  >  1:2840096/1‑56 (MQ=255)
                                                aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcg                                  >  1:3684352/1‑56 (MQ=255)
                                                aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTa                       >  1:2715422/1‑67 (MQ=255)
                                                aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCa                    <  1:610035/70‑1 (MQ=255)
                                                aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCa                    <  1:3859286/70‑1 (MQ=255)
                                                aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCa                    >  1:791053/1‑70 (MQ=255)
                                                aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAc                   >  1:3444189/1‑71 (MQ=255)
                                                aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAc                   >  1:1326101/1‑71 (MQ=255)
                                                  cAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgc                 >  1:3674322/1‑71 (MQ=255)
                                                  cAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgc                 >  1:3726392/1‑71 (MQ=255)
                                                    gCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                                            >  1:795070/1‑42 (MQ=255)
                                                     cAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgctgc              <  1:3762712/71‑1 (MQ=255)
                                                      aaaGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTa                                       >  1:3439588/1‑45 (MQ=255)
                                                      aaaGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgctgct             >  1:2463292/1‑71 (MQ=255)
                                                         gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc                           <  1:4078005/54‑1 (MQ=255)
                                                         gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgct                <  1:3670763/65‑1 (MQ=255)
                                                         gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCa          >  1:1941883/1‑71 (MQ=255)
                                                         gTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACgat                                        <  1:154663/41‑3 (MQ=12)
                                                          ttCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                         <  1:3844695/39‑1 (MQ=255)
                                                          ttCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                         <  1:128481/39‑1 (MQ=255)
                                                          ttCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTa                                       >  1:455223/1‑41 (MQ=255)
                                                             aGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGa       <  1:3871500/70‑1 (MQ=255)
                                                             aGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGa       <  1:1720907/70‑1 (MQ=255)
                                                                 cGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTg  <  1:1330184/71‑1 (MQ=255)
                                                                  |                                                                     
ACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTG  >  minE/1439750‑1439885

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: