Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 228,318 A→G 8.9% F129L (TTC→CTC) ‡ phoE ← outer membrane phosphoporin protein E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE228,3180AG 8.9% 175.9 / 7.3 67F129L (TTC→CTC) ‡phoEouter membrane phosphoporin protein E
Reads supporting (aligned to +/- strand):  ref base A (34/27);  new base G (2/4);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

GCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCATAATCGAAAGAACCCAAATCT  >  minE/228248‑228388
                                                                      |                                                                      
gCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGaa                                                                        <  1:4514319/71‑1 (MQ=255)
         cTGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCgg                                                                >  1:5018650/1‑70 (MQ=255)
         cTGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGt                                                               >  1:1537934/1‑71 (MQ=255)
         cTGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGt                                                               >  1:67333/1‑71 (MQ=255)
            gcgcGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTcc                                                             <  1:5148096/70‑1 (MQ=255)
            gcgcGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTcc                                                             <  1:2699299/70‑1 (MQ=255)
            gcgcGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTcc                                                             <  1:2705790/70‑1 (MQ=255)
            gcgcGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTcc                                                             <  1:3800690/70‑1 (MQ=255)
            gcgcGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCa                                                            <  1:2094244/71‑1 (MQ=255)
                gTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCt                                                        >  1:4802152/1‑71 (MQ=255)
                gTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATAGCCACCAAATTCCGGGAACATATCGGTCAAGGCt                                                        >  1:1168160/1‑71 (MQ=255)
                 tttGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGaa                                                                        <  1:1095715/54‑1 (MQ=255)
                 tttGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGaa                                                                        <  1:2978170/54‑1 (MQ=255)
                 tttGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGaa                                                                        <  1:1959262/54‑1 (MQ=255)
                   tGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTc                                                      >  1:228010/1‑70 (MQ=255)
                   tGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTc                                                      >  1:4382516/1‑70 (MQ=255)
                   tGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTc                                                      >  1:473025/1‑70 (MQ=255)
                       cATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCGGGGAACATATCGGTCCAGGCTTcc                                                     <  1:3334581/67‑1 (MQ=255)
                       cATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTcc                                                     <  1:3826850/67‑1 (MQ=255)
                        aTAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGt                                                 >  1:3327623/1‑70 (MQ=255)
                        aTAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGt                                                 >  1:4214848/1‑70 (MQ=255)
                        aTAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGt                                                 >  1:2943974/1‑70 (MQ=255)
                        aTAAAGTTGTCGGTCTGCGACGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGt                                                 >  1:5477579/1‑70 (MQ=255)
                         tAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGaa                                                                        <  1:1215067/46‑1 (MQ=255)
                         tAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGaa                                                                        <  1:2691582/46‑1 (MQ=255)
                         tAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGaa                                                                        <  1:1212664/46‑1 (MQ=255)
                         tAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCa                                               <  1:2178694/71‑1 (MQ=255)
                          aaaGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACa                                                                      >  1:1730454/1‑47 (MQ=255)
                           aaGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGc                                                         >  1:2091362/1‑59 (MQ=255)
                            aGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATAc                                            <  1:410217/71‑1 (MQ=255)
                                gTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATa                                             >  1:681500/1‑66 (MQ=255)
                                 tCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGTGAACATATCGGTCCAg                                                           <  1:200812/51‑1 (MQ=255)
                                 tCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACg                                        <  1:469558/70‑1 (MQ=255)
                                 tCGGTCTGCGCCGAGGAATCGCCAACAAATTCCGGGAACATATCGGTCCAg                                                           <  1:1776079/51‑1 (MQ=255)
                                    gTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACatat                                                                   <  1:3615088/40‑1 (MQ=255)
                                     tCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATAc                                            <  1:4283894/62‑1 (MQ=255)
                                       tGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAg                                 >  1:3546922/1‑71 (MQ=255)
                                         cgcCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAg                                 <  1:3694876/69‑1 (MQ=255)
                                         cgcCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAg                                 <  1:1326057/69‑1 (MQ=255)
                                         cgcCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGt                               >  1:3631384/1‑71 (MQ=255)
                                           ccGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTa                             >  1:3812058/1‑71 (MQ=255)
                                              aGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGa                          >  1:3286098/1‑71 (MQ=255)
                                              aGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGa                          >  1:4810085/1‑71 (MQ=255)
                                                gAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGccccc                                   <  1:4421983/60‑1 (MQ=255)
                                                 aaTCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGTCCa                       <  1:1473738/71‑1 (MQ=255)
                                                  aTCGCCACCAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATa                                             >  1:5650054/1‑48 (MQ=14)
                                                  aTCGCCACCAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATa                                             >  1:3918596/1‑48 (MQ=14)
                                                  aTCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCa                       >  1:5461541/1‑70 (MQ=255)
                                                   tCGCCACCAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATa                     <  1:3091891/71‑1 (MQ=2)
                                                   tCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCATa                     <  1:223126/71‑1 (MQ=255)
                                                   tCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCATa                     <  1:1489484/71‑1 (MQ=255)
                                                    cGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGccc                                     >  1:5226986/1‑54 (MQ=255)
                                                    cGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGccc                                     >  1:1562158/1‑54 (MQ=255)
                                                     gCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACaa                                          >  1:368172/1‑48 (MQ=255)
                                                     gCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACaa                                          >  1:4533659/1‑48 (MQ=255)
                                                      ccaccaAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCa                                               >  1:4851033/1‑42 (MQ=255)
                                                       caccaAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCATAATc                  >  1:2613759/1‑70 (MQ=255)
                                                        acctAATTCCGGGAACTTTTCGGTACAGGCTTCCTCGTCATAc                                            <  1:5594220/43‑1 (MQ=19)
                                                        accaAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCATAATCGa                <  1:5584933/71‑1 (MQ=255)
                                                         ccaAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGt                               <  1:1350133/55‑1 (MQ=255)
                                                         ccaAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCATAATCGaa               <  1:23206/71‑1 (MQ=255)
                                                          caAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACaa                                          >  1:4896614/1‑43 (MQ=255)
                                                          caAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACaa                                          >  1:5352103/1‑43 (MQ=255)
                                                           aaaTTCCGGGAACATATCGGTCCAGGCTTCCACGTCATa                                             >  1:3223565/1‑39 (MQ=255)
                                                           aaaTTCCGGGAACATATCGGTCCAGGCTTCCACGTCATa                                             >  1:2031837/1‑39 (MQ=255)
                                                           aaaTTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCATAATCGa                >  1:5473907/1‑68 (MQ=255)
                                                           aaaTTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAg             <  1:886751/71‑1 (MQ=255)
                                                           aaaTTCAGGCAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAg             <  1:4996157/71‑1 (MQ=255)
                                                            aaTTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCATAATCGAAAg             <  1:2273794/70‑1 (MQ=255)
                                                             aTTCAGGGAGCATATCGGTCCAGCCTTCGATGTCGTACATTACGCCATAATTACGTCCGTAATCGAAAGa            <  1:1140555/70‑1 (MQ=255)
                                                               tCCGGGAACATATCGGTCCAGGCTTCCACGTCATa                                             >  1:4964631/1‑35 (MQ=255)
                                                               tCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTAc                            >  1:694195/1‑52 (MQ=255)
                                                                   ggAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCATAATCGAAAGAACCCaaa     <  1:3284920/71‑1 (MQ=255)
                                                                   ggAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCATAATCGAAAGAACCCaaa     <  1:1674353/71‑1 (MQ=255)
                                                                      aCATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCATAATCGAAAGAACCCAAATCt  >  1:1659261/1‑71 (MQ=255)
                                                                      |                                                                      
GCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCATAATCGAAAGAACCCAAATCT  >  minE/228248‑228388

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: