Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,231,014	C→T	G581D (GGT→GAT)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,231,014	0	C	T	100.0%	150.3 / NA	50	G581D (GGT→GAT)	ptrB	protease II
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/29);  total (21/29)

TGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAA  >  minE/1230947‑1231080                                                                    |                                                                   tgCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTAtcc                                                                  <  1:1119284/70‑1 (MQ=255)   cGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCAc                                                               >  1:2210289/1‑71 (MQ=255)      tGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCActct                                                            <  1:2451602/71‑1 (MQ=255)      tGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCActct                                                            <  1:1097293/71‑1 (MQ=255)         cGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCa                                                         <  1:2675826/71‑1 (MQ=255)            tGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCaaa                                                       >  1:179844/1‑70 (MQ=255)            tGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAAc                                                      >  1:2617599/1‑71 (MQ=255)            tGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAAc                                                      >  1:1263846/1‑71 (MQ=255)             gTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCtt                                                           >  1:478809/1‑65 (MQ=255)             gTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACt                                                     >  1:2949214/1‑71 (MQ=255)                aTACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCAc                                                  <  1:156173/71‑1 (MQ=255)                     ggCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTg                                             <  1:1734370/71‑1 (MQ=255)                            gCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACt                                                     >  1:2864884/1‑56 (MQ=255)                            gCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAagag                                      >  1:2904837/1‑71 (MQ=255)                            gCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAagag                                      >  1:443777/1‑71 (MQ=255)                            gCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAagag                                      <  1:2292313/71‑1 (MQ=255)                            gCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAaga                                       <  1:1916449/70‑1 (MQ=255)                             cTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAgg                                     <  1:959895/71‑1 (MQ=255)                                  cATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAAtt                                 <  1:560344/70‑1 (MQ=255)                                  cATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAAtt                                 <  1:223972/70‑1 (MQ=255)                                  cATGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAAtt                                 >  1:2723222/1‑70 (MQ=255)                                   aTGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGa                               <  1:2164927/71‑1 (MQ=255)                                   aTGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGa                               <  1:2726978/71‑1 (MQ=255)                                    tGTACTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGAt                              <  1:1057251/71‑1 (MQ=255)                                        cTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATtca                           <  1:2161720/70‑1 (MQ=255)                                        cTCGTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATtca                           <  1:620318/70‑1 (MQ=255)                                           gTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATtcatc                         >  1:2460207/1‑69 (MQ=255)                                           gTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATtcatc                         <  1:650787/69‑1 (MQ=255)                                           gTAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATtcatc                         <  1:2071247/69‑1 (MQ=255)                                            tAATATTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAg                      >  1:9002/1‑71 (MQ=255)                                              atatTGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGc                     <  1:2116513/70‑1 (MQ=255)                                                 ttGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCa                                                >  1:2781559/1‑40 (MQ=255)                                                 ttGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGTAATTGATTCATCAAGCATc                  <  1:2795913/70‑1 (MQ=255)                                                 ttGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATtca                           <  1:958700/61‑1 (MQ=255)                                                 ttGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATc                  <  1:2109378/70‑1 (MQ=255)                                                 ttGCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCg                 >  1:1554721/1‑71 (MQ=255)                                                   gCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAAtt                                 <  1:1063160/53‑1 (MQ=255)                                                   gCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCgt                <  1:441489/70‑1 (MQ=255)                                                   gCGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCgt                <  1:954064/70‑1 (MQ=255)                                                    cGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCa                                                >  1:884176/1‑37 (MQ=255)                                                    cGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCgttg              <  1:1180885/71‑1 (MQ=255)                                                    cGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCgttg              >  1:2635757/1‑71 (MQ=255)                                                    cGGATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCgttg              >  1:2632225/1‑71 (MQ=255)                                                     ggATCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCgttgt             >  1:2429800/1‑71 (MQ=255)                                                        tCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGAt                              <  1:951598/51‑1 (MQ=255)                                                        tCCTGCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTac          <  1:946034/71‑1 (MQ=255)                                                            gCGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTacaaca      >  1:2887252/1‑71 (MQ=255)                                                             cGGGTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTacaaca      <  1:870506/70‑1 (MQ=255)                                                               ggTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATca   <  1:2542824/71‑1 (MQ=255)                                                               ggTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATc    >  1:1072048/1‑70 (MQ=255)                                                                gTTATCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATcaa  >  1:227334/1‑71 (MQ=255)                                                                    |                                                                   TGCGGTGACGTTGTCATACGGGCTGTAGCTTTTCATGTACTCGTAATATTGCGGATCCTGCGGGTTACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAA  >  minE/1230947‑1231080

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