Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,153,878	G→A	T246I (ACA→ATA)	astE ←	succinylglutamate desuccinylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,153,878	0	G	A	100.0%	225.9 / NA	67	T246I (ACA→ATA)	astE	succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/27);  total (40/27)

TATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTG  >  minE/1153813‑1153945                                                                  |                                                                    tatGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCATCCGATAACGGAGCGGCGGTATtct                                                                 >  1:1756262/1‑70 (MQ=255)         tCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATtct                                                                 <  1:1632178/62‑1 (MQ=255)         tCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATtc                                                                  <  1:2321052/61‑1 (MQ=255)          cGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATAcc                                                        >  1:1027113/1‑70 (MQ=255)          cGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATAcc                                                        >  1:2137574/1‑70 (MQ=255)             aGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGaca                                                    >  1:1855233/1‑71 (MQ=255)             aGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGaca                                                    >  1:1355040/1‑71 (MQ=255)              ggACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGacac                                                   >  1:335087/1‑71 (MQ=255)               gACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACAct                                                  <  1:313009/71‑1 (MQ=255)                   gcgagtggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGaca                                                    <  1:2697351/65‑1 (MQ=255)                   gcgagtggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGaca                                                    <  1:1459947/65‑1 (MQ=255)                   gcgagtggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGaca                                                    <  1:2687659/65‑1 (MQ=255)                     gagtggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACAct                                                  >  1:560735/1‑65 (MQ=255)                      agtggcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATtc                                                                  <  1:2072844/48‑1 (MQ=255)                          gcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATa                                                          >  1:977879/1‑52 (MQ=255)                          gcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATa                                                          >  1:122190/1‑52 (MQ=255)                          gcgagtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATa                                                          >  1:991650/1‑52 (MQ=255)                             agtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACActct                                                >  1:3070625/1‑59 (MQ=255)                             agtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATag                                      >  1:3053409/1‑69 (MQ=255)                              gtAATTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATagcag                                   >  1:513704/1‑71 (MQ=255)                                 aTTTGCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAgcagcg                                 <  1:2772068/70‑1 (MQ=255)                                     gCGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAgcagcgcaatt                            >  1:1265080/1‑68 (MQ=255)                                      cGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCa                                              <  1:1067502/52‑1 (MQ=255)                                      cGAAACCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCa                                              <  1:466718/52‑1 (MQ=255)                                        aaaCCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCAcc                                            >  1:1701952/1‑52 (MQ=255)                                        aaaCCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCAcc                                            >  1:1974194/1‑52 (MQ=255)                                        aaaCCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCAcc                                            >  1:2977920/1‑52 (MQ=255)                                        aaaCCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAgcagcg                                 <  1:69421/63‑1 (MQ=255)                                        aaaCCACCCGATAACGCAGCGGCGGTATTCTCACGATACCGACACTCTCAcc                                            >  1:2665975/1‑52 (MQ=255)                                         aaCCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTg                        >  1:634252/1‑71 (MQ=255)                                         aaCCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTg                        >  1:296446/1‑71 (MQ=255)                                         aaCCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTg                        <  1:707184/71‑1 (MQ=255)                                         aaCCACCCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTg                        >  1:762482/1‑71 (MQ=255)                                          accaccCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAgcagcgca                               >  1:3012147/1‑63 (MQ=255)                                          accaccCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAgcagcgca                               >  1:639997/1‑63 (MQ=255)                                          accaccCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTgc                       >  1:2182919/1‑71 (MQ=255)                                          accaccCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTgc                       >  1:1700676/1‑71 (MQ=255)                                          accaccCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTg                        >  1:2659713/1‑70 (MQ=255)                                          accaccCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTg                        >  1:2031266/1‑70 (MQ=255)                                           ccaccCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTgcg                      <  1:263876/71‑1 (MQ=255)                                           ccaccCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTgcg                      <  1:2817129/71‑1 (MQ=255)                                           ccaccCGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTgcg                      <  1:1113056/71‑1 (MQ=255)                                              cccGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCaa                           >  1:1542582/1‑63 (MQ=255)                                               ccGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCgcagca                            <  1:2683277/61‑1 (MQ=255)                                               ccGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCgcagca                            <  1:1748647/61‑1 (MQ=255)                                               ccGATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCgcagca                            <  1:2085188/61‑1 (MQ=255)                                                 gATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTgg                  >  1:845205/1‑69 (MQ=255)                                                 gATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTgg                  >  1:1331034/1‑69 (MQ=255)                                                 gATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTgg                  >  1:2217828/1‑69 (MQ=255)                                                 gATAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTg                   >  1:1691240/1‑68 (MQ=255)                                                 gATAACGGAGCGGCGCTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCaa                           >  1:570114/1‑60 (MQ=255)                                                  aTAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTgc                       <  1:3028699/63‑1 (MQ=255)                                                  aTAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTgc                       <  1:1722476/63‑1 (MQ=255)                                                  aTAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTgg                  >  1:716629/1‑68 (MQ=255)                                                  aTAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAg               >  1:1358312/1‑71 (MQ=255)                                                   tAACGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGt              >  1:1672498/1‑71 (MQ=255)                                                    aaCGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGt              >  1:2604037/1‑70 (MQ=255)                                                     aCGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGtt             <  1:1437072/70‑1 (MQ=255)                                                     aCGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTa            <  1:790771/71‑1 (MQ=255)                                                      cGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCgcagca                            <  1:996992/54‑1 (MQ=255)                                                      cGGAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTAc           <  1:1353649/71‑1 (MQ=255)                                                       ggAGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACt          <  1:760881/71‑1 (MQ=255)                                                         aGCGGCGGTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTg         >  1:309127/1‑70 (MQ=255)                                                             gcggTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAAc    <  1:2887083/71‑1 (MQ=255)                                                              cggTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACt   >  1:1445553/1‑71 (MQ=255)                                                              cggTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACt   >  1:2128085/1‑71 (MQ=255)                                                               ggTATTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTg  <  1:1733724/71‑1 (MQ=255)                                                                  |                                                                    TATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTG  >  minE/1153813‑1153945

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: ‑