Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,023,066	G→A	P281S (CCA→TCA)	ygiQ ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,023,066	0	G	A	100.0%	274.2 / NA	85	P281S (CCA→TCA)	ygiQ	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/45);  total (40/45)

CCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGG  >  minE/2023002‑2023129                                                                 |                                                                cccTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCAc                                                           <  1:846427/71‑1 (MQ=255)      cACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGtt                                                       <  1:4734229/70‑1 (MQ=255)      cACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGtt                                                       <  1:2790779/70‑1 (MQ=255)      cACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTa                                                      >  1:4894154/1‑71 (MQ=255)      cACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTa                                                      >  1:2530656/1‑71 (MQ=255)      cACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTa                                                      >  1:855058/1‑71 (MQ=255)      cACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTa                                                      >  1:1464673/1‑71 (MQ=255)       aCTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTAc                                                     <  1:1461604/71‑1 (MQ=255)       aCTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTAc                                                     <  1:2000048/71‑1 (MQ=255)       aCTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTAc                                                     <  1:2080246/71‑1 (MQ=255)       aCTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTAc                                                     <  1:3486325/71‑1 (MQ=255)       aCTTTCTAGAAATAAGGCAGCAACACTTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTAc                                                     <  1:3847077/71‑1 (MQ=255)        cTTTCTCGAAAGAAGGCAGCGACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTAc                                                     <  1:733204/70‑1 (MQ=255)          ttCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGc                                                  >  1:897353/1‑71 (MQ=255)          ttCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGc                                                  >  1:75155/1‑71 (MQ=255)          ttCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGc                                                  >  1:2563966/1‑71 (MQ=255)          ttCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGc                                                  >  1:4546264/1‑71 (MQ=255)           tctcGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGc                                                  >  1:2153393/1‑70 (MQ=255)           tctcGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCt                                                 >  1:3886636/1‑71 (MQ=255)           tctcGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCt                                                 >  1:3107764/1‑71 (MQ=255)           tctcGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCt                                                 >  1:362729/1‑71 (MQ=255)           tctcGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCt                                                 >  1:2347694/1‑71 (MQ=255)           tctcGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCt                                                 >  1:3468177/1‑71 (MQ=255)           tctcGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCt                                                 >  1:1230581/1‑71 (MQ=255)                 aagaagGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCt                                           <  1:930353/71‑1 (MQ=255)                 aagaagGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCt                                           <  1:1116117/71‑1 (MQ=255)                 aagaagGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCt                                           <  1:4240050/71‑1 (MQ=255)                  agaagGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCGCGGTTACGGCTTTGGCtt                                          >  1:4443202/1‑71 (MQ=255)                      ggCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTg                                      <  1:4478278/71‑1 (MQ=255)                      ggCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTg                                      <  1:749781/71‑1 (MQ=255)                      ggCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTg                                      <  1:1736434/71‑1 (MQ=255)                      ggCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTg                                      <  1:2761782/71‑1 (MQ=255)                      ggCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTg                                      <  1:1851233/71‑1 (MQ=255)                      ggCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTg                                      <  1:917862/71‑1 (MQ=255)                      ggCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTg                                      <  1:3034460/71‑1 (MQ=255)                      ggCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTg                                      <  1:3075254/71‑1 (MQ=255)                      ggCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTg                                      <  1:557731/71‑1 (MQ=255)                       gcagcaACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTcc                                        >  1:1317430/1‑68 (MQ=255)                       gcagcaACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTg                                      >  1:4614386/1‑70 (MQ=255)                          gcaACAGGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCttt                                  <  1:3653496/71‑1 (MQ=255)                          gcaACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCttt                                  <  1:415231/71‑1 (MQ=255)                          gcaACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCttt                                  <  1:1481864/71‑1 (MQ=255)                          gcaACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGGGACTGCACGGTTACGGCTTTGGCTTCCTGCttt                                  <  1:39727/71‑1 (MQ=255)                           caACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCttt                                  <  1:529705/70‑1 (MQ=255)                            aaCACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCttttt                                >  1:3856936/1‑71 (MQ=255)                            aaCACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCttttt                                >  1:3293800/1‑71 (MQ=255)                            aaCACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCttttt                                >  1:3344262/1‑71 (MQ=255)                                cGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGt                            >  1:2320137/1‑71 (MQ=255)                                  tAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGc                          >  1:824603/1‑71 (MQ=255)                                  tAGGTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGc                          >  1:226477/1‑71 (MQ=255)                                    ggTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCa                        >  1:167893/1‑71 (MQ=255)                                    ggTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCa                        >  1:2419183/1‑71 (MQ=255)                                    ggTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCa                        >  1:317274/1‑71 (MQ=255)                                    ggTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCa                        >  1:3047532/1‑71 (MQ=255)                                    ggTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCa                        >  1:3038919/1‑71 (MQ=255)                                    ggTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCa                        >  1:602140/1‑71 (MQ=255)                                    ggTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCa                        >  1:2111526/1‑71 (MQ=255)                                    ggTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCa                        >  1:2150865/1‑71 (MQ=255)                                    ggTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCa                        >  1:4362882/1‑71 (MQ=255)                                    ggTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCa                        >  1:2218001/1‑71 (MQ=255)                                    ggTTTTTTCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCa                        >  1:4843118/1‑71 (MQ=255)                                        ttttCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTg                                              >  1:4621285/1‑45 (MQ=255)                                        ttttCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTg                                              >  1:3716772/1‑45 (MQ=255)                                        ttttCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTg                                              >  1:614074/1‑45 (MQ=255)                                        ttttCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTg                                              >  1:1675286/1‑45 (MQ=255)                                        ttttCCCACGGTTTCGGGCGCGGTGACTGCACGGTTACGGCTTTg                                              >  1:373994/1‑45 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGc                          <  1:4601509/43‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGc                          <  1:1085928/43‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGc                          <  1:1107077/43‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGc                          <  1:778125/43‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGc                          <  1:1691550/43‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGc                          <  1:1943283/43‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGc                          <  1:3220771/43‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGc                          <  1:2882101/43‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGc                          <  1:2827278/43‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGc                          <  1:4476431/43‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGc                          <  1:3407328/43‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACgg  <  1:60299/67‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACgg  <  1:1754191/67‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACgg  <  1:4256097/67‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACgg  <  1:2436253/67‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACgg  <  1:3816206/67‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACgg  <  1:1028176/67‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACgg  <  1:242082/67‑1 (MQ=255)                                                              gTGACTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACgg  <  1:1037743/67‑1 (MQ=255)                                                                 |                                                                CCCTTCACTTTCTCGAAAGAAGGCAGCAACACGTAGGTTTTTTCCCACGGTTTCGGGCGCGGTGGCTGCACGGTTACGGCTTTGGCTTCCTGCTTTTTCGGTGCCACCGGTTTGTTATCCGCGCACGG  >  minE/2023002‑2023129

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