Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 223,016 G→A 14.9% intergenic (+71/‑109) 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
*minE223,0160GA14.9% 188.0 / 20.1 67intergenic (+71/‑109)lpcA/prfHD‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base A (4/6);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCC  >  minE/222949‑223084
                                                                   |                                                                    
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)
   cGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGt                                                                >  1:238771/1‑71 (MQ=14)
    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)
        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)
                   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:3586635/1‑70 (MQ=255)
                       tACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATc                                             <  1:110868/70‑1 (MQ=255)
                       tACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCt                                            >  1:3760552/1‑71 (MQ=255)
                       tACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCt                                            >  1:3697971/1‑71 (MQ=255)
                       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)
                               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:2405781/71‑1 (MQ=255)
                                   ggCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg                                >  1:2075294/1‑71 (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:2458922/50‑1 (MQ=11)
                                        ttGTGCGCAGATGCCTGATGCGACGCTAGCGCGTCTTATCATGCCTACaa                                                <  1:1963291/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)
                                                  atgcctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaataat                   <  1:1202693/69‑1 (MQ=255)
                                                  atgcctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGc                 >  1:2373905/1‑71 (MQ=255)
                                                  atgcctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGc                 >  1:2193495/1‑71 (MQ=255)
                                                   tgcctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg                                <  1:1379248/55‑1 (MQ=255)
                                                   tgcctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCg                <  1:1940905/71‑1 (MQ=255)
                                                    gcctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGt               >  1:3771337/1‑71 (MQ=255)
                                                    gcctgatgcGACGCTAACGCGACTTATCATGCCTACa                                                 >  1:1237910/1‑37 (MQ=11)
                                                     cctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGtt              <  1:1022089/71‑1 (MQ=255)
                                                     cctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGtt              <  1:1202913/71‑1 (MQ=255)
                                                     cctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGtt              <  1:1570863/71‑1 (MQ=255)
                                                        gatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaataa                    >  1:1806158/1‑62 (MQ=255)
                                                        gatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaat                      <  1:1101179/60‑1 (MQ=255)
                                                         atgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCAcg          >  1:302684/1‑71 (MQ=255)
                                                          tgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCAcgc         >  1:1669211/1‑71 (MQ=255)
                                                                 cTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGt               >  1:1703369/1‑58 (MQ=255)
                                                                 cTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCAcgc         >  1:47448/1‑64 (MQ=255)
                                                                 cTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATcc  >  1:79132/1‑71 (MQ=255)
                                                                   gACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTa                              >  1:4079423/1‑41 (MQ=255)
                                                                   |                                                                    
TGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCC  >  minE/222949‑223084

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: