Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 222,991 G→A 10.0% intergenic (+46/‑134) lpcA → / → prfH D‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE222,9910GA10.0% 153.1 / 8.6 70intergenic (+46/‑134)lpcA/prfHD‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor
Reads supporting (aligned to +/- strand):  ref base G (35/28);  new base A (3/4);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCC  >  minE/222927‑223060
                                                                |                                                                     
aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                                 <  1:3333644/71‑1 (MQ=255)
aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                                 >  1:155804/1‑71 (MQ=255)
aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                                 <  1:1090087/71‑1 (MQ=255)
 aaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctg                                                                <  1:764669/71‑1 (MQ=255)
   agagATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgat                                                              >  1:1296964/1‑71 (MQ=255)
   agagATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgat                                                              >  1:3852053/1‑71 (MQ=255)
     agaTGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcg                                                                   >  1:1077818/1‑64 (MQ=255)
        tGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctg                                                         <  1:1196182/71‑1 (MQ=255)
              agtaagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCg                                                    <  1:3505789/70‑1 (MQ=255)
                  agtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCt                                               <  1:1655151/71‑1 (MQ=255)
                     cTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAc                                            >  1:434380/1‑71 (MQ=255)
                         tgTAGGCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGcagatgcct                                                          <  1:2467046/50‑1 (MQ=9)
                         tgTAGGCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCt                                      <  1:4325091/70‑1 (MQ=9)
                         tgTAGGCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCt                                      <  1:3991590/70‑1 (MQ=9)
                             ggCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcgc                                                                  >  1:1101933/1‑41 (MQ=255)
                             ggCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTa                                    >  1:4039942/1‑71 (MQ=255)
                             ggCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTa                                    <  1:3968490/71‑1 (MQ=255)
                             ggCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTAt                                   >  1:2766014/1‑70 (MQ=9)
                             ggCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTAt                                   <  1:4582333/70‑1 (MQ=9)
                             ggCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTAt                                   <  1:4683535/70‑1 (MQ=9)
                             ggCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTAt                                   <  1:3118110/70‑1 (MQ=9)
                             ggCCGGACAAGACGTTTA‑‑TCGAATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTAt                                   <  1:834604/70‑1 (MQ=2)
                              gCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTa                                    >  1:4619478/1‑70 (MQ=11)
                                cGGATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGt                                        <  1:2226868/64‑1 (MQ=14)
                                 ggATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGcagatgcctgatgc                                                     <  1:1756742/50‑1 (MQ=11)
                                 ggATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGcagatgcctgatgc                                                     <  1:422928/50‑1 (MQ=11)
                                 ggATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTAt                                   >  1:1063915/1‑68 (MQ=14)
                                 ggACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTATc                                  <  1:3573976/67‑1 (MQ=9)
                                  gATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCagg                               >  1:4350560/1‑69 (MQ=11)
                                  gATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCagg                               >  1:5554197/1‑69 (MQ=11)
                                  gATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCagg                               >  1:824945/1‑69 (MQ=11)
                                  gATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcct                                                          <  1:1701772/44‑1 (MQ=255)
                                   aTAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctga                                                        >  1:1805980/1‑45 (MQ=255)
                                   aTAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctga                                                        >  1:2609340/1‑45 (MQ=255)
                                   aTAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctga                                                        >  1:3616525/1‑45 (MQ=255)
                                    tAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGcc                             <  1:3814930/71‑1 (MQ=11)
                                    tAAGACGTTTACGCCGCATCCGGCATTTGTGCGctga                                                               >  1:4290447/1‑37 (MQ=255)
                                      aggcGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTa                           <  1:1875605/68‑1 (MQ=11)
                                      aggcGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCt                            <  1:75157/67‑1 (MQ=11)
                                      aGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTa                           >  1:3666333/1‑71 (MQ=255)
                                      aGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTa                           >  1:989136/1‑71 (MQ=255)
                                       ggcGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTAc                          <  1:5436686/69‑1 (MQ=11)
                                          gTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACaaa                       >  1:1943052/1‑71 (MQ=14)
                                          gTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATg                               >  1:1076548/1‑63 (MQ=255)
                                          gTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGc                              >  1:2237443/1‑64 (MQ=255)
                                          gTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACaaa                       >  1:2742367/1‑71 (MQ=255)
                                           tttATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATc                                  <  1:3095830/59‑1 (MQ=11)
                                           ttcaCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGACCTACAAAt                      <  1:316021/68‑1 (MQ=11)
                                            ttATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACaaact                     >  1:187324/1‑69 (MQ=11)
                                            ttACGCCGCATCCGGCATTTGTGCGATGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATc                     >  1:2157744/1‑71 (MQ=255)
                                             tACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCt                    >  1:3034541/1‑71 (MQ=255)
                                             tACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCt                    >  1:5551855/1‑71 (MQ=255)
                                             tACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGc                                                >  1:5012285/1‑43 (MQ=14)
                                             caCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACAAAtt                     >  1:3545990/2‑69 (MQ=11)
                                             caCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGACCTACAAAttc                    >  1:2380082/2‑69 (MQ=11)
                                             caCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGACCTACAAAttc                    >  1:4035262/2‑69 (MQ=11)
                                               cgccgcATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACAAAttct                   <  1:3513868/70‑4 (MQ=11)
                                                 ccgcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAc                >  1:4984833/1‑71 (MQ=255)
                                                 ccgcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAc                <  1:5228399/71‑1 (MQ=255)
                                                  cgcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAcg               <  1:4007667/71‑1 (MQ=255)
                                                  cgcATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGt                                        >  1:557805/1‑46 (MQ=14)
                                                   gcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAcgc              >  1:3970778/1‑71 (MQ=255)
                                                   gcATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACa                         <  1:1360600/60‑1 (MQ=11)
                                                    cATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTAcg                         <  1:1850183/58‑2 (MQ=14)
                                                      tCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACAAAttc                    >  1:1944324/1‑60 (MQ=11)
                                                           cATTTGTGCGCTGATGCCTGATGCGACGCTGCCGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg        >  1:1816658/1‑69 (MQ=255)
                                                           cATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg        >  1:2460913/1‑69 (MQ=255)
                                                           cATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg        >  1:374981/1‑69 (MQ=255)
                                                           cATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTa      <  1:1381827/71‑1 (MQ=255)
                                                             tttGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACaa                        <  1:1505484/51‑1 (MQ=255)
                                                              ctgtgtgCCGATGCCTGATGCGACGCTGACGCGTTTTATCAt                                <  1:1633189/41‑1 (MQ=14)
                                                              ctgtgtgCCGATGCCTGATGCGACGCTGACGCGTTTTATCAt                                <  1:3307470/41‑1 (MQ=14)
                                                               tgtgtgCCGATGCCTGATGCGACGCTGACGCGTTTTAt                                   <  1:4718543/38‑1 (MQ=14)
                                                               tgtgtgCCGATGCCTGATGCGACGCTGACGCGTTTTAt                                   <  1:599282/38‑1 (MQ=14)
                                                               tgtgCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAg     >  1:3118635/1‑68 (MQ=255)
                                                               tgtgCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAg     >  1:2547384/1‑68 (MQ=255)
                                                               tgtgCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGcc  >  1:2958194/1‑71 (MQ=255)
                                                                gtgtGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGc                              >  1:1834998/1‑42 (MQ=14)
                                                                gtgtGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGc                              >  1:1661702/1‑42 (MQ=14)
                                                                gtgtGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGc                              >  1:4043461/1‑42 (MQ=14)
                                                                gtgtGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGc                              >  1:2721273/1‑42 (MQ=14)
                                                                gtgtGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTa                           >  1:3473134/1‑45 (MQ=11)
                                                                gtgtGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTa                           >  1:3616852/1‑45 (MQ=11)
                                                                gtgCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACAAAtt                     >  1:4971846/1‑50 (MQ=14)
                                                                |                                                                     
AAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCC  >  minE/222927‑223060

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: