Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 222,988 T→G 15.0% intergenic (+43/‑137) 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,9880TG15.0% 186.6 / 33.4 80intergenic (+43/‑137)lpcA/prfHD‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor
Reads supporting (aligned to +/- strand):  ref base T (30/38);  new base G (4/8);  total (34/46)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGG  >  minE/222919‑223058
                                                                     |                                                                      
ttgattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAtt                                                                       <  1:3093897/71‑1 (MQ=255)
ttgattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAtt                                                                       >  1:2127581/1‑71 (MQ=255)
ttgattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAtt                                                                       <  1:3735513/71‑1 (MQ=255)
ttgattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAtt                                                                       <  1:3834809/71‑1 (MQ=255)
ttgattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAtt                                                                       >  1:1733439/1‑71 (MQ=255)
  gattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAt                                                                        <  1:1617161/68‑1 (MQ=255)
  gattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTtg                                                                     >  1:1142411/1‑71 (MQ=255)
  gattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTtg                                                                     >  1:3933237/1‑71 (MQ=255)
  gattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTtg                                                                     >  1:421712/1‑71 (MQ=255)
    ttgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAt                                                                        <  1:1052485/66‑1 (MQ=255)
    ttgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTtgtg                                                                   >  1:3354382/1‑71 (MQ=255)
      gaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcg                                                                 >  1:67990/1‑71 (MQ=255)
       aaaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcg                                                                 <  1:544005/70‑1 (MQ=255)
        aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAttt                                                                      >  1:1431471/1‑64 (MQ=255)
        aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                               <  1:2110079/71‑1 (MQ=255)
        aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                               <  1:2535133/71‑1 (MQ=255)
        aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                               <  1:3471632/71‑1 (MQ=255)
         aaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctg                                                              <  1:298310/71‑1 (MQ=255)
         aaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                               >  1:1579794/1‑70 (MQ=255)
           agagATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgat                                                            >  1:1190847/1‑71 (MQ=255)
               aTGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcct                                                        <  1:404501/71‑1 (MQ=255)
               aTGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcct                                                        >  1:3590817/1‑71 (MQ=255)
                tGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTtgt                                                                    <  1:2990156/58‑1 (MQ=255)
                tGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTtgt                                                                    <  1:2221802/58‑1 (MQ=255)
                 ggTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatg                                                           >  1:510023/1‑66 (MQ=255)
                 ggTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatg                                                           >  1:421964/1‑66 (MQ=255)
                   ttAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatg                                                    >  1:430694/1‑71 (MQ=255)
                   ttAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatg                                                    >  1:1424715/1‑71 (MQ=255)
                   ttAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatg                                                    >  1:2509624/1‑71 (MQ=255)
                      agtaagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcct                                                        <  1:1477924/64‑1 (MQ=255)
                      agtaagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCg                                                  <  1:2909023/70‑1 (MQ=255)
                      agtaagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCg                                                  <  1:389116/70‑1 (MQ=255)
                       gtaagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatg                                                           <  1:2619229/60‑1 (MQ=255)
                       gtaagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                               <  1:3301149/56‑1 (MQ=255)
                       gtaagtCTGGCCTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                               <  1:1945690/56‑1 (MQ=255)
                        taagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcg                                                                 <  1:1414670/53‑1 (MQ=255)
                         aagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcct                                                        <  1:3887419/61‑1 (MQ=255)
                          agtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAtt                                                                       >  1:1508497/1‑45 (MQ=255)
                          agtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgc                                                          <  1:4090573/58‑1 (MQ=255)
                              tGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAcg                                         >  1:3618317/1‑71 (MQ=255)
                              tGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAcg                                         >  1:3623050/1‑71 (MQ=255)
                                gCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGa                                           >  1:1322626/1‑67 (MQ=255)
                                gCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAcgc                                        <  1:413923/70‑1 (MQ=255)
                                gCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAc                                          <  1:3192193/68‑1 (MQ=255)
                                 tgTAGGCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCt                                             <  1:2882678/61‑1 (MQ=9)
                                 cGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGt                                      >  1:238771/1‑71 (MQ=14)
                                 cGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                               <  1:1250631/46‑1 (MQ=255)
                                  gTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                               <  1:2129690/45‑1 (MQ=255)
                                  gTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTc                                     <  1:1687108/71‑1 (MQ=255)
                                   tAGGCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTa                                  >  1:2702424/1‑71 (MQ=12)
                                   tAGGCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTa                                  <  1:2568575/71‑1 (MQ=12)
                                    aGGCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTa                                  <  1:398636/70‑1 (MQ=12)
                                     ggCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTa                                  <  1:2526099/71‑1 (MQ=11)
                                     ggCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCtt                                   <  1:3840587/70‑1 (MQ=14)
                                     ggCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgat                                                     <  1:575863/52‑1 (MQ=255)
                                      gCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcct                                                        <  1:733193/48‑1 (MQ=255)
                                      gCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctga                                                             >  1:1327297/1‑43 (MQ=255)
                                      gCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCt                                    <  1:1453856/68‑1 (MQ=255)
                                      gCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTAt                                 <  1:3147032/71‑1 (MQ=255)
                                      gCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTAt                                 <  1:719785/71‑1 (MQ=255)
                                         ggATAAGACGTTTACGCCGCATCCGGCATTTGTgc                                                                  <  1:2301882/35‑1 (MQ=255)
                                         ggATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCAt                              >  1:3875956/1‑71 (MQ=255)
                                            tAAGACGTTTACGCCGCATCCGGCATTTGTGCGctga                                                             >  1:2564245/1‑37 (MQ=255)
                                              aGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctg                                                       <  1:3147794/41‑1 (MQ=255)
                                                 cGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACa                       >  1:2312175/1‑70 (MQ=255)
                                                 cGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACa                       >  1:2348861/1‑70 (MQ=255)
                                                   tttATGCCGCATCCGGCAGTAATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTAt                                 <  1:3129315/58‑1 (MQ=11)
                                                   tttACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACaaa                     <  1:1937058/70‑1 (MQ=255)
                                                   tttACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAAt                    <  1:1172589/71‑1 (MQ=255)
                                                     tACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATc                   <  1:680451/70‑1 (MQ=255)
                                                     tACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATc                   <  1:110868/70‑1 (MQ=255)
                                                     tACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATc                   >  1:3586635/1‑70 (MQ=255)
                                                     tACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCt                  >  1:3760552/1‑71 (MQ=255)
                                                     tACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCt                  >  1:3697971/1‑71 (MQ=255)
                                                     tACGCCGCATCCGGCAGTTGTACGcaggtgcctga                                                      >  1:2082527/1‑35 (MQ=9)
                                                     tACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACAAAttc                  >  1:1879168/1‑69 (MQ=11)
                                                       cgccgcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTg                 <  1:20578/70‑1 (MQ=255)
                                                         ccgcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTa               <  1:1513807/70‑1 (MQ=255)
                                                           gcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAcgc            >  1:249184/1‑71 (MQ=255)
                                                            cATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAcgc            >  1:912032/1‑70 (MQ=255)
                                                            cATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCt                                             <  1:505567/37‑1 (MQ=11)
                                                             aTCCGGCAGTCGTGTGCTGATGCCTGATGCGACGCTAACGCGTCTTAt                                 <  1:37824/48‑1 (MQ=14)
                                                                cGGCATTTGTGCGCTGATGCCTGATGCGACGCTGa                                           >  1:1050908/1‑35 (MQ=255)
                                                                cGGCATTTGTGCGCTGATGCCTGATGCGACGCTGCCGCGTCTTATCATGCCTACAAATCTGTACGCGAAcc       >  1:2331005/1‑71 (MQ=255)
                                                                cGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAAcc       >  1:3142606/1‑71 (MQ=255)
                                                                 ggCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg      >  1:2075294/1‑71 (MQ=255)
                                                                 ggCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg      <  1:2405781/71‑1 (MQ=255)
                                                                 ggCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg      <  1:832762/71‑1 (MQ=255)
                                                                 ggCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg      <  1:917430/71‑1 (MQ=255)
                                                                     tttGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATc                   <  1:418304/54‑1 (MQ=255)
                                                                     tttGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTg                 <  1:2126657/56‑1 (MQ=255)
                                                                     tttGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCATACAAATCTGTACGCGAACCGTAgg  >  1:1900840/1‑71 (MQ=255)
                                                                     |                                                                      
TTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGG  >  minE/222919‑223058

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: