Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 258,840 A→G 6.7% V162A (GTT→GCT)  phoE ← outer membrane phosphoporin protein E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb258,8400AG 6.7% 162.5 / 3.5 60V162A (GTT→GCT) phoEouter membrane phosphoporin protein E
Reads supporting (aligned to +/- strand):  ref base A (20/36);  new base G (2/2);  total (22/38)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

TTGCTTTTTAACGTCGCGGTTTTCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGTCGGTCT  >  W3110S.gb/258773‑258909
                                                                   |                                                                     
ttGCTTTTTAACGTCGCGGTTTTCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACg                                                                     >  1:230387/1‑70 (MQ=255)
ttGCTTTTTAACGTCGCGGTTTTCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGc                                                                    >  1:376662/1‑71 (MQ=255)
 tGCTTTTTAACGTCGCGGTTTTCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATaa                                                                       >  1:4271730/1‑67 (MQ=255)
   cTTTTTAACGTCGCGGTTTTCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGcc                                                                   <  1:587124/69‑1 (MQ=255)
   cTTTTTAACGTCGCGGTTTTCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCg                                                                  <  1:123505/70‑1 (MQ=255)
            gTCGCGGTTTTCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCgg                                                        <  1:2634819/71‑1 (MQ=255)
              cgcgGTTTTCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAGCGCCGaagaag                                                            >  1:1000412/1‑65 (MQ=255)
              cgcgGTTTTCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGgtg                                                      >  1:1568231/1‑71 (MQ=255)
               gcgGTTTTCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGgtgt                                                     >  1:3254357/1‑71 (MQ=255)
                cgGTTTTCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCgaa                                                                <  1:2119046/59‑1 (MQ=255)
                  gTTTTCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACg                                                                     <  1:2588565/52‑1 (MQ=255)
                  gTTTTCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACg                                                                     <  1:2534614/52‑1 (MQ=255)
                  gTTTTCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTcc                                                  >  1:3715385/1‑71 (MQ=255)
                   ttttCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAGCGCCga                                                                 <  1:2690051/55‑1 (MQ=255)
                   ttttCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCgaagaa                                                             >  1:2599531/1‑59 (MQ=255)
                   ttttCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCg                                                 <  1:3131749/71‑1 (MQ=255)
                   ttttCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCAATCGATAACGCCGAAGAAGTCGGTGTTCCg                                                 <  1:1677286/71‑1 (MQ=255)
                    tttCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGa                                                >  1:426187/1‑71 (MQ=255)
                     ttCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGa                                                <  1:2323921/70‑1 (MQ=255)
                     ttCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGa                                                >  1:2618232/1‑70 (MQ=255)
                       cGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATaa                                                                       >  1:3577325/1‑45 (MQ=255)
                        gTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACg                                            >  1:1098605/1‑71 (MQ=255)
                         tttttCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGt                                           <  1:1601246/71‑1 (MQ=255)
                           tttCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTc                                          >  1:2278167/1‑70 (MQ=255)
                           tttCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCg                                         >  1:4065144/1‑71 (MQ=255)
                            ttCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCgaag                                                               <  1:2075106/48‑1 (MQ=255)
                             tCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGc                                        <  1:348167/70‑1 (MQ=255)
                             tCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGc                                        >  1:3340796/1‑70 (MQ=255)
                               ccTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCa                                      <  1:3270080/70‑1 (MQ=255)
                                cTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGa                                    <  1:987498/71‑1 (MQ=255)
                                 ttGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGAc                                   >  1:2660947/1‑71 (MQ=255)
                                  tGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGAcc                                  >  1:1848636/1‑71 (MQ=255)
                                   gATATTGCAGGGTTAAGTTCAGTCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCg                                 >  1:867549/1‑71 (MQ=255)
                                    atatTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCg                                 <  1:1741215/70‑1 (MQ=255)
                                     tatTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCt                               <  1:2775899/71‑1 (MQ=255)
                                     tatTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCt                               <  1:3105611/71‑1 (MQ=255)
                                       ttGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTcc                                                  <  1:1547462/50‑1 (MQ=255)
                                       ttGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGcc                                       <  1:1399461/61‑1 (MQ=255)
                                          cAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCAGATACGTCGCCAGACCGCTGgcg                           <  1:752136/70‑1 (MQ=255)
                                            gggTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGt                        >  1:1363006/1‑71 (MQ=255)
                                               ttAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATa                                              <  1:2053550/46‑1 (MQ=255)
                                               ttAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATa                                              <  1:1941705/46‑1 (MQ=255)
                                                 aaGTTCAGGCCATCGATAGCGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTg                     >  1:736714/1‑69 (MQ=255)
                                                 aaGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTgg                    >  1:3904167/1‑70 (MQ=255)
                                                   gTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGAt                                               <  1:3120031/41‑1 (MQ=255)
                                                    ttCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTTGTCAt                <  1:2186719/71‑1 (MQ=255)
                                                      cAGGCCATCGATAGCGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATa               <  1:4169959/70‑1 (MQ=255)
                                                      cAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCAtt               >  1:3142443/1‑69 (MQ=255)
                                                      cAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATa               <  1:42082/70‑1 (MQ=255)
                                                       aGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATaa              <  1:599733/70‑1 (MQ=255)
                                                       aGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATaa              <  1:2638379/70‑1 (MQ=255)
                                                       aGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATaa              <  1:3096537/70‑1 (MQ=255)
                                                        ggCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATaaa             <  1:2975654/70‑1 (MQ=255)
                                                         gCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATAAAGt           <  1:2338479/71‑1 (MQ=255)
                                                           cATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATAAAGtt          <  1:2560045/70‑1 (MQ=255)
                                                            aTCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGt        <  1:4281570/71‑1 (MQ=255)
                                                            aTCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGt        <  1:302432/71‑1 (MQ=255)
                                                            aTCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTAGCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGt        <  1:2231305/71‑1 (MQ=255)
                                                               gATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGTCgg     >  1:94732/1‑71 (MQ=255)
                                                               gATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGTCg      <  1:1082537/70‑1 (MQ=255)
                                                                aTAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATAAAg            <  1:830819/63‑1 (MQ=255)
                                                                aTAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGTCgg     >  1:167783/1‑70 (MQ=255)
                                                                   aCGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGTCGGTCt  <  1:2711709/70‑1 (MQ=255)
                                                                   |                                                                     
TTGCTTTTTAACGTCGCGGTTTTCGTTTTTCCCTTGATATTGCAGGGTTAAGTTCAGGCCATCGATAACGCCGAAGAAGTCGGTGTTCCGATACGTCGCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGTCGGTCT  >  W3110S.gb/258773‑258909

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: