Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 222,997 T→A 22.8% intergenic (+52/‑128) 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,9970TA22.8% 116.0 / 30.2 53intergenic (+52/‑128)lpcA/prfHD‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor
Reads supporting (aligned to +/- strand):  ref base T (17/23);  new base A (3/9);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

AAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAAT  >  minE/222927‑223064
                                                                      |                                                                   
aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                                     <  1:3471632/71‑1 (MQ=255)
aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                                     <  1:2535133/71‑1 (MQ=255)
aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                                     <  1:2110079/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:3590817/1‑71 (MQ=255)
       aTGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcct                                                              <  1:404501/71‑1 (MQ=255)
         ggTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatg                                                                 >  1:421964/1‑66 (MQ=255)
         ggTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatg                                                                 >  1:510023/1‑66 (MQ=255)
           ttAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatg                                                          >  1:2509624/1‑71 (MQ=255)
           ttAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatg                                                          >  1:1424715/1‑71 (MQ=255)
           ttAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatg                                                          >  1:430694/1‑71 (MQ=255)
              agtaagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcct                                                              <  1:1477924/64‑1 (MQ=255)
              agtaagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCg                                                        <  1:389116/70‑1 (MQ=255)
              agtaagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCg                                                        <  1:2909023/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)
                 aagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcct                                                              <  1:3887419/61‑1 (MQ=255)
                  agtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgc                                                                <  1:4090573/58‑1 (MQ=255)
                      tGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAcg                                               >  1:3623050/1‑71 (MQ=255)
                      tGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAcg                                               >  1:3618317/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:2568575/71‑1 (MQ=12)
                           tAGGCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTa                                        >  1:2702424/1‑71 (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)
                                 ggATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCAt                                    >  1:3875956/1‑71 (MQ=255)
                                    tAAGACGTTTACGCCGCATCCGGCATTTGTGCGctga                                                                   >  1:2564245/1‑37 (MQ=255)
                                      aGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctg                                                             <  1:3147794/41‑1 (MQ=255)
                                         cGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACa                             >  1:2348861/1‑70 (MQ=255)
                                         cGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACa                             >  1:2312175/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:3697971/1‑71 (MQ=255)
                                             tACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCt                        >  1:3760552/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:917430/71‑1 (MQ=255)
                                                         ggCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg            <  1:832762/71‑1 (MQ=255)
                                                         ggCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg            >  1:2075294/1‑71 (MQ=255)
                                                         ggCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg            <  1:2405781/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)
                                                              ttGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACaa                            <  1:1963291/50‑1 (MQ=11)
                                                              ttGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACaa                            <  1:2458922/50‑1 (MQ=11)
                                                                gtgtGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGc                                  >  1:2080797/1‑42 (MQ=14)
                                                                gtgCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGcc      <  1:3783617/70‑1 (MQ=255)
                                                                  gcgcTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCg     <  1:3882356/69‑1 (MQ=255)
                                                                   cgcTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaat  <  1:2928515/71‑1 (MQ=255)
                                                                   cgcTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaa   >  1:72883/1‑70 (MQ=255)
                                                                      tgatgcctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGt                      <  1:1416463/48‑1 (MQ=255)
                                                                      |                                                                   
AAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAAT  >  minE/222927‑223064

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: