Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 222,990 T→C 12.4% intergenic (+45/‑135) 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,9900TC12.4% 191.2 / 17.8 74intergenic (+45/‑135)lpcA/prfHD‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor
Reads supporting (aligned to +/- strand):  ref base T (27/37);  new base C (3/6);  total (30/44)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGG  >  minE/222921‑223058
                                                                     |                                                                    
gattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTtg                                                                     >  1:3933237/1‑71 (MQ=255)
gattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTtg                                                                     >  1:1142411/1‑71 (MQ=255)
gattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTtg                                                                     >  1:421712/1‑71 (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:2535133/71‑1 (MQ=255)
      aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                               <  1:3471632/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:404501/71‑1 (MQ=255)
             aTGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcct                                                        >  1:3590817/1‑71 (MQ=255)
              tGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTtgt                                                                    <  1:2221802/58‑1 (MQ=255)
              tGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTtgt                                                                    <  1:2990156/58‑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:430694/1‑71 (MQ=255)
                 ttAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatg                                                    >  1:1424715/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)
                        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)
                                       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: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: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: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)
                                                                     |                                                                    
GATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGG  >  minE/222921‑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: