Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 228,316 G→C 6.7% F129L (TTC→TTG) ‡ phoE ← outer membrane phosphoporin protein E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported228,3160GC 6.7% 224.8 / 4.7 74F129L (TTC→TTG) ‡phoEouter membrane phosphoporin protein E
Reads supporting (aligned to +/- strand):  ref base G (34/35);  new base C (2/3);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCATAATCGAAAGAACCCAAA  >  Exported/228248‑228385
                                                                    |                                                                     
gCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGaa                                                                     <  1:4514319/71‑1 (MQ=255)
gCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGa                                                                      <  1:1480948/70‑1 (MQ=255)
gCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGa                                                                      <  1:2022930/70‑1 (MQ=255)
         cTGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCgg                                                             >  1:5018650/1‑70 (MQ=255)
         cTGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGt                                                            >  1:67333/1‑71 (MQ=255)
         cTGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGt                                                            >  1:1537934/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:2978170/54‑1 (MQ=255)
                 tttGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGaa                                                                     <  1:1959262/54‑1 (MQ=255)
                 tttGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGaa                                                                     <  1:1095715/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:2943974/1‑70 (MQ=255)
                        aTAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGt                                              >  1:4214848/1‑70 (MQ=255)
                        aTAAAGTTGTCGGTCTGCGACGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGt                                              >  1:5477579/1‑70 (MQ=255)
                         tAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGaa                                                                     <  1:2691582/46‑1 (MQ=255)
                         tAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGaa                                                                     <  1:1215067/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:1489484/71‑1 (MQ=255)
                                                   tCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCATa                  <  1:223126/71‑1 (MQ=255)
                                                    cGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGccc                                  >  1:1562158/1‑54 (MQ=255)
                                                    cGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGccc                                  >  1:5226986/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)
                                                                    |                                                                     
GCCAGACCGCTGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCCAGGTTACGACCATAATCGAAAGAACCCAAA  >  Exported/228248‑228385

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: