Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,798 A→C 20.5% V168V (GTT→GTG ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,7980AC20.5% 244.5 / 71.8 127V168V (GTT→GTGompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (44/57);  new base C (16/10);  total (60/67)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 134 total aligned reads displayed.

GACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAC  >  minE/1439731‑1439868
                                                                   |                                                                      
gACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGaa                                                                       >  1:5640637/1‑69 (MQ=255)
    gttgttAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCa                                                                  <  1:2955158/70‑1 (MQ=255)
     ttgttAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCa                                                                  >  1:2907751/1‑69 (MQ=255)
       gttagttaCGCCACTATTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAg                                                               <  1:4691544/70‑1 (MQ=255)
       gttagttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAg                                                               <  1:867747/70‑1 (MQ=255)
              aCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGcc                                                       >  1:2776364/1‑71 (MQ=255)
                    cTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                                 >  1:901969/1‑71 (MQ=255)
                    cTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                                 >  1:4028265/1‑71 (MQ=255)
                     tAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACa                                                                     >  1:5563721/1‑50 (MQ=255)
                        tAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAg                                                               <  1:946638/53‑1 (MQ=255)
                             ccTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                  >  1:2616335/1‑61 (MQ=255)
                               ttCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaag                                       >  1:3829570/1‑70 (MQ=255)
                                 cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCa                                                           >  1:3144134/1‑48 (MQ=255)
                                 cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                  >  1:2568228/1‑57 (MQ=255)
                                 cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                  >  1:2246951/1‑57 (MQ=255)
                                 cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                  >  1:4472787/1‑57 (MQ=255)
                                 cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                    >  1:4820904/1‑71 (MQ=255)
                                 cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                    >  1:4628469/1‑71 (MQ=255)
                                 cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                    >  1:4258669/1‑71 (MQ=255)
                                 cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                    >  1:1775763/1‑71 (MQ=255)
                                  aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                   <  1:4362760/71‑1 (MQ=255)
                                  aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                   <  1:5050667/71‑1 (MQ=255)
                                   ccAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                   >  1:4219703/1‑70 (MQ=255)
                                     aGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCa                                                                  <  1:3519884/37‑1 (MQ=255)
                                      gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                                 >  1:1593116/1‑53 (MQ=255)
                                      gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaag                                       <  1:4076419/63‑1 (MQ=255)
                                      gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaag                                       <  1:2264004/63‑1 (MQ=255)
                                      gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgt                               >  1:346052/1‑71 (MQ=255)
                                      gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgt                               >  1:1500305/1‑71 (MQ=255)
                                        tGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                             <  1:594962/71‑1 (MQ=255)
                                          ggTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                           <  1:2939076/71‑1 (MQ=255)
                                            ttGCCGTTTTTACCCTTGTACTGAACAGCAAAGtt                                                             <  1:3747349/35‑1 (MQ=255)
                                            ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                             <  1:2788405/35‑1 (MQ=255)
                                            ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                             <  1:3148972/35‑1 (MQ=255)
                                            ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                             <  1:194325/35‑1 (MQ=255)
                                            ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCa                                 <  1:5095300/63‑1 (MQ=255)
                                            ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCa                                 <  1:2691678/63‑1 (MQ=255)
                                              ttcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAgt                               <  1:1222257/61‑1 (MQ=2)
                                              ttcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACga                        >  1:245171/3‑69 (MQ=255)
                                              gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                        <  1:2365396/70‑1 (MQ=255)
                                               ccGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                        <  1:3448231/53‑1 (MQ=255)
                                                 gTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAgg                    >  1:3481802/1‑71 (MQ=255)
                                                 gTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAgg                    >  1:4051572/1‑71 (MQ=255)
                                                  tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                  >  1:4609455/1‑40 (MQ=255)
                                                  tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAgg                    >  1:3235629/1‑70 (MQ=255)
                                                  tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt                   <  1:5418776/71‑1 (MQ=255)
                                                  ttgttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGt                   <  1:4163262/68‑1 (MQ=9)
                                                  ttattACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGt                   >  1:944137/4‑71 (MQ=255)
                                                  ttattACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGt                   >  1:5490431/4‑71 (MQ=255)
                                                   tattACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGtt                  >  1:2472725/3‑70 (MQ=255)
                                                    tttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCa                                               >  1:2124127/1‑41 (MQ=255)
                                                    tttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCa                                               >  1:5206816/1‑41 (MQ=255)
                                                    tttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcg                 >  1:2293748/1‑71 (MQ=255)
                                                    tttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcg                 >  1:5234665/1‑71 (MQ=255)
                                                    gttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcg                 >  1:2122834/2‑71 (MQ=255)
                                                    gttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTc                  <  1:1903393/69‑1 (MQ=255)
                                                    gttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTc                  <  1:4359271/69‑1 (MQ=255)
                                                     ttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc                <  1:3486189/71‑1 (MQ=9)
                                                     ttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                        <  1:4326007/47‑1 (MQ=255)
                                                     ttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                        <  1:755752/47‑1 (MQ=255)
                                                      tACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc                >  1:1925733/1‑70 (MQ=9)
                                                      tACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc                >  1:1432098/1‑70 (MQ=9)
                                                       aCCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg               >  1:77650/1‑70 (MQ=255)
                                                       aCCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg               >  1:2255134/1‑70 (MQ=255)
                                                       aCCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgcaa              >  1:4607544/1‑69 (MQ=9)
                                                        cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGCAGAAGTCAGTGTTACGGTAGGTCGCGaa             <  1:4488359/71‑1 (MQ=255)
                                                        cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGtt                            <  1:1886572/56‑1 (MQ=255)
                                                        cccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgcaa              <  1:446008/70‑3 (MQ=255)
                                                           tgatACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGc           >  1:4504168/4‑70 (MQ=9)
                                                           tgatACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGc           >  1:1718842/4‑70 (MQ=9)
                                                            ggTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                           >  1:2497059/1‑53 (MQ=255)
                                                            ggTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCg         >  1:1261773/1‑71 (MQ=255)
                                                            gatACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAGTCAGTATTACGATAAGTCGCGACGcca         >  1:3367110/3‑70 (MQ=9)
                                                            gatACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGcca         >  1:5310415/3‑70 (MQ=9)
                                                             gTACTGAACAGCAAAGTTCAGGCCTTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt        <  1:2621710/71‑1 (MQ=255)
                                                             gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc          <  1:888923/69‑1 (MQ=255)
                                                             gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc          <  1:1149380/69‑1 (MQ=255)
                                                             gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCg         <  1:4788670/70‑1 (MQ=255)
                                                             gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt        <  1:798438/71‑1 (MQ=255)
                                                             gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt        >  1:1097808/1‑71 (MQ=255)
                                                             gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt        <  1:4839570/71‑1 (MQ=255)
                                                             gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt        >  1:2312503/1‑71 (MQ=255)
                                                             gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt        >  1:3983863/1‑71 (MQ=255)
                                                              tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCGGTGTTACGGTAGGTCGCGAAGc           <  1:5030844/67‑1 (MQ=255)
                                                              tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt       <  1:5506495/71‑1 (MQ=255)
                                                              tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt       <  1:2037482/71‑1 (MQ=255)
                                                               aCTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGaa             >  1:5538641/1‑64 (MQ=255)
                                                                cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                        <  1:4804220/52‑1 (MQ=255)
                                                                cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                        <  1:2206403/52‑1 (MQ=255)
                                                                 tGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                        <  1:5518907/35‑1 (MQ=255)
                                                                 tGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                        <  1:5617523/35‑1 (MQ=255)
                                                                 tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTc                                  <  1:1269071/41‑1 (MQ=255)
                                                                 tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTc                                  <  1:3644829/41‑1 (MQ=255)
                                                                 tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                                <  1:1015214/43‑1 (MQ=255)
                                                                 tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt       <  1:1523523/68‑1 (MQ=255)
                                                                  gAACAGCAAAGTTCAGGCCGTCAACCAGACCgaag                                       <  1:3798414/35‑1 (MQ=255)
                                                                  gAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCa                                 <  1:4521908/41‑1 (MQ=255)
                                                                   aaCAGCAAAGTTCAGGCCGTGAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGaa             >  1:4705293/1‑60 (MQ=255)
                                                                   aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGaa             >  1:1589268/1‑60 (MQ=255)
                                                                   aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAc  >  1:4472473/1‑71 (MQ=255)
                                                                   |                                                                      
GACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAC  >  minE/1439731‑1439868

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: