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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,8130GA 6.0% 379.6 / 9.3 116G163G (GGC→GGTompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (64/44);  new base A (3/4);  total (67/49)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 117 total aligned reads displayed.

ACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGA  >  minE/1439745‑1439880
                                                                    |                                                                   
aCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGcc                                                                   >  1:2028957/1‑71 (MQ=255)
 cGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGcc                                                                   <  1:2429539/70‑1 (MQ=255)
     aCTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTc                                                                <  1:3982698/69‑1 (MQ=255)
     aCTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                              <  1:2590613/71‑1 (MQ=255)
      cTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                                             <  1:1435957/71‑1 (MQ=255)
      cTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                                             >  1:511295/1‑71 (MQ=255)
       tAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                                             >  1:3321781/1‑70 (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:2042108/1‑63 (MQ=255)
               ccTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAcc                                                            >  1:3573862/1‑63 (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:2022812/71‑1 (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:3487694/1‑71 (MQ=255)
                          tGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                                         >  1:856898/1‑71 (MQ=255)
                           gggTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAgggc                                                                   >  1:1191602/1‑42 (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)
                                gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                    <  1:3717165/70‑1 (MQ=255)
                                gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                    >  1:3802273/1‑70 (MQ=255)
                                gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                    <  1:1424771/70‑1 (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)
                                      attACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGttg                             >  1:3924542/2‑69 (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)
                                          cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                                         <  1:382139/55‑1 (MQ=255)
                                          cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGaa                         <  1:3749021/71‑1 (MQ=255)
                                           ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg                        <  1:3279985/71‑1 (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:3092619/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:2716471/1‑55 (MQ=255)
                                               gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                    >  1:3478811/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:2004825/71‑1 (MQ=255)
                                                  cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                 <  1:1489981/71‑1 (MQ=255)
                                                  cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                 >  1:1429344/1‑71 (MQ=255)
                                                  cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                 >  1:3807269/1‑71 (MQ=255)
                                                  cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                 <  1:3296023/71‑1 (MQ=255)
                                                  cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                 <  1:3002901/71‑1 (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:1760466/1‑53 (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:3684352/1‑56 (MQ=255)
                                                     aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcg                             >  1:2840096/1‑56 (MQ=255)
                                                     aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTa                  >  1:2715422/1‑67 (MQ=255)
                                                     aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCa               <  1:610035/70‑1 (MQ=255)
                                                     aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCa               >  1:791053/1‑70 (MQ=255)
                                                     aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAc              >  1:1326101/1‑71 (MQ=255)
                                                       cAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgc            >  1:3726392/1‑71 (MQ=255)
                                                         gCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                                       >  1:795070/1‑42 (MQ=255)
                                                          cAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGc                       <  1:1241078/57‑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)
                                                              gTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACgat                                   <  1:154663/41‑3 (MQ=12)
                                                               ttCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                    <  1:128481/39‑1 (MQ=255)
                                                               ttCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                    <  1:3844695/39‑1 (MQ=255)
                                                               ttCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTa                                  >  1:455223/1‑41 (MQ=255)
                                                                  aGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt                               <  1:81904/41‑1 (MQ=255)
                                                                  aGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGa  <  1:3871500/70‑1 (MQ=255)
                                                                    |                                                                   
ACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGA  >  minE/1439745‑1439880

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: