Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,268,188 T→G 53.7% A139A (GCA→GCC pgsA ← phosphatidylglycerophosphate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,268,1880TG53.7% 24.5 / 86.3 82A139A (GCA→GCCpgsAphosphatidylglycerophosphate synthetase
Reads supporting (aligned to +/- strand):  ref base T (36/2);  new base G (32/12);  total (68/14)
Fisher's exact test for biased strand distribution p-value = 8.97e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGACGGCCACGCTACTGCGTT  >  minE/1268126‑1268255
                                                              |                                                                   
gCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAAc                                                               <  1:370670/69‑1 (MQ=255)
      aagaagAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGccacc                                                         >  1:1393411/1‑69 (MQ=255)
      aaGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGccacc                                                         >  1:812970/1‑69 (MQ=255)
      aaGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGccacc                                                         >  1:756212/1‑69 (MQ=255)
      aaGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGccacc                                                         >  1:453783/1‑69 (MQ=255)
      aaGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGccacc                                                         >  1:2602654/1‑69 (MQ=255)
      aaGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGccacc                                                         >  1:1243722/1‑69 (MQ=255)
      aaGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGccacc                                                         >  1:1337783/1‑69 (MQ=255)
      aaGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGccacc                                                         >  1:1634636/1‑69 (MQ=255)
      aaGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGccacc                                                         >  1:1461076/1‑69 (MQ=255)
      aaGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGccacc                                                         >  1:3175254/1‑69 (MQ=255)
      aaGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGccacc                                                         >  1:1551607/1‑69 (MQ=255)
                cAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAg                                               >  1:2629523/1‑69 (MQ=255)
                cAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAg                                               >  1:343754/1‑69 (MQ=255)
                cAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAg                                               >  1:2203007/1‑69 (MQ=255)
                cAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAg                                               >  1:3357684/1‑69 (MQ=255)
                cAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAg                                               >  1:3309944/1‑69 (MQ=255)
                cAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAg                                               >  1:3000442/1‑69 (MQ=255)
                cAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAg                                               >  1:2903608/1‑69 (MQ=255)
                cAATACCGGCGTACTCAACCCAAATGTTCGGACCCCACAGCAGCCATGCCAACGCCACCATCTGGGCAg                                               >  1:1676471/1‑69 (MQ=255)
                     ccGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGtt                                          >  1:1834352/1‑69 (MQ=255)
                     ccGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGtt                                          >  1:2207105/1‑69 (MQ=255)
                     ccGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGtt                                          >  1:2594632/1‑69 (MQ=255)
                     ccGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGtt                                          >  1:2920648/1‑69 (MQ=255)
                     ccGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGtt                                          >  1:2694483/1‑69 (MQ=255)
                     ccGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGtt                                          >  1:857167/1‑69 (MQ=255)
                     ccGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGtt                                          >  1:471524/1‑69 (MQ=255)
                     ccGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGtt                                          >  1:470502/1‑69 (MQ=255)
                     ccGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGt                                           >  1:578630/1‑68 (MQ=255)
                     ccGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGt                                           >  1:2350496/1‑68 (MQ=255)
                     ccGGCGTACTCAACCCAAATGTTAGGACGCCGCAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGt                                           >  1:464065/1‑68 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGTTTTTCACttt                                  >  1:1389678/1‑69 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACttt                                  >  1:3395194/1‑69 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACttt                                  >  1:45516/1‑69 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACttt                                  >  1:1664920/1‑69 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACttt                                  >  1:2347407/1‑69 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACttt                                  >  1:2962497/1‑69 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGccacca                                                        >  1:1069198/1‑47 (MQ=39)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGc                                                 >  1:2992945/1‑54 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGc                                                 >  1:470161/1‑54 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:2966251/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:3372684/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:591875/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:910723/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:874526/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:359364/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:1890306/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:2681795/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:107491/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:2486107/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:1811507/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:2400996/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:1183374/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:1202202/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:1676622/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:15233/1‑55 (MQ=255)
                             cTCAACCCAAATGTTCGGACGCCACAGCAGCCAGCCTAACGCCACCATCTgg                                                   >  1:3206286/1‑52 (MQ=255)
                             cTCAACCCAAATGTTAGGACGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCa                                                >  1:1310950/1‑55 (MQ=255)
                                            cGGACGCCATAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAc                   >  1:2851054/1‑69 (MQ=255)
                                            cGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAc                   >  1:1586206/1‑69 (MQ=255)
                                            cGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAc                   >  1:585792/1‑69 (MQ=255)
                                            cGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAc                   >  1:1107532/1‑69 (MQ=255)
                                            cGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAc                   >  1:1897691/1‑69 (MQ=255)
                                            cGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAc                   >  1:2131043/1‑69 (MQ=255)
                                            cGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGAc                   >  1:21420/1‑69 (MQ=255)
                                                cGCCACAGCAGCCAGGCTTACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATc                           <  1:1237328/57‑1 (MQ=255)
                                                cGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATc                           <  1:2162867/57‑1 (MQ=255)
                                                cGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATc                           <  1:426156/57‑1 (MQ=255)
                                                cGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATc                           <  1:1036112/57‑1 (MQ=255)
                                                cGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATc                           <  1:3336816/57‑1 (MQ=255)
                                                cGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATc                           <  1:1634229/57‑1 (MQ=255)
                                                cGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATc                           <  1:2411797/57‑1 (MQ=255)
                                                cGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATc                           <  1:230622/57‑1 (MQ=255)
                                                cGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATc                           <  1:2390480/57‑1 (MQ=255)
                                                cGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATc                           <  1:3238690/57‑1 (MQ=255)
                                                cGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATc                           <  1:2338734/57‑1 (MQ=255)
                                                cGCCACAGCAGCCAGGCTAACGCCACCATCTGGGCAGTGGTTCTCACTTTCCCAATc                           <  1:3243663/57‑1 (MQ=255)
                                                         agcCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGACGGCCACGCTACt       >  1:2847506/1‑68 (MQ=255)
                                                         agcCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGACGGCCACGCTACTg      >  1:1135422/1‑69 (MQ=255)
                                                         agcCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGACGGCCACGCTACTg      >  1:2222606/1‑69 (MQ=255)
                                                            cATGCCAACGCCACCATCTGGGCAGTTGTTTTCACTTTCCCAATCCAGGAGACGGCCACGCTACTGCg    >  1:565634/1‑68 (MQ=255)
                                                             aTGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGACGGCCACGCTACTGCGtt  <  1:2626592/69‑1 (MQ=255)
                                                              tGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGACGGCCACGCTACTGCGtt  <  1:809355/68‑1 (MQ=255)
                                                              tGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGACGGCCACGCTACTGCGtt  <  1:3316855/68‑1 (MQ=255)
                                                              |                                                                   
GCCACAAAGAAAAGTGCAATACCGGCGTACTCAACCCAAATGTTCGGACGCCACAGCAGCCATGCCAACGCCACCATCTGGGCAGTGGTTTTCACTTTCCCAATCCAGGAGACGGCCACGCTACTGCGTT  >  minE/1268126‑1268255

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: