Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,921,568	+T	coding (455/1197 nt)	yjgN →	conserved inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,921,564	1	.	T	100.0%	87.8 / NA	27	I151F (ATT→TTT)	yjgN	conserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (17/10);  total (17/10)

TCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGA‑TTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAA  >  minE/2921496‑2921625                                                                      |                                                              tcAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGA‑tt                                                             >  1:826403/1‑71 (MQ=255)    aCGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGA‑ttttt                                                          <  1:1098955/71‑2 (MQ=255)    aCGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGA‑ttttt                                                          <  1:3007673/71‑2 (MQ=255)    aCGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGA‑tttt                                                           >  1:2277711/1‑70 (MQ=255)      ggCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGA‑tttttaa                                                        <  1:3061095/71‑4 (MQ=255)      ggCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGA‑tttttaa                                                        <  1:2499526/71‑4 (MQ=255)            aGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTTAATGGCCa                                                   <  1:1944660/71‑1 (MQ=255)                ttAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTTAATGGCCATTgg                                               >  1:1166962/1‑71 (MQ=255)                ttAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTTAATGGCCATTgg                                               >  1:414640/1‑71 (MQ=255)                ttAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTTAATGGCCATTgg                                               >  1:2879482/1‑71 (MQ=255)                ttAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTTAATGGCCATTgg                                               >  1:104402/1‑71 (MQ=255)                 tAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTTAATGGCCATTgg                                               <  1:2421223/70‑1 (MQ=255)                              tctATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGttt                                 >  1:2576099/1‑71 (MQ=255)                                tATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGA‑tttt                                                           <  1:1797478/42‑1 (MQ=255)                                  tGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCtt                             >  1:2226667/1‑71 (MQ=255)                                  tGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCtt                             >  1:775962/1‑71 (MQ=255)                                  tGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCtt                             >  1:2946472/1‑71 (MQ=255)                                  tGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCt                              <  1:656645/70‑1 (MQ=255)                                    aaaGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTa                           >  1:2318041/1‑71 (MQ=255)                                      aGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTAtc                         >  1:2820129/1‑71 (MQ=255)                                      aGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTAtc                         <  1:3170595/71‑1 (MQ=255)                                               gtgggtgACCTTTTTCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAg                <  1:532281/71‑1 (MQ=255)                                                    tgACCTTTTTCTTGCCGATTTTTAATGGCCATTGGGATGGGGa                                       >  1:3116811/1‑43 (MQ=38)                                                      aCCTTTTTCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTAc         >  1:2515472/1‑71 (MQ=255)                                                      aCCTTTTTCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTAc         >  1:2551116/1‑71 (MQ=255)                                                      aCCTTTTTCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTATTTATCTCGACAAAGATGCTAc         >  1:1855767/1‑71 (MQ=255)                                                        cTTTTTCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTAt                          <  1:354214/52‑1 (MQ=255)                                                          ttttCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTAt                          >  1:1922870/1‑50 (MQ=255)                                                          ttttCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTAt                          >  1:281606/1‑50 (MQ=255)                                                             tCTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCa   <  1:946328/70‑1 (MQ=255)                                                             tCTTGCCGATTTTTAATGGCCATTGGGATGGGCACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCa   <  1:2169648/70‑1 (MQ=255)                                                              cTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCat  <  1:2113491/70‑2 (MQ=255)                                                              cTTGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCat  <  1:2526987/70‑2 (MQ=255)                                                               ttGCCGATTTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCaa  >  1:704162/1‑69 (MQ=255)                                                                      |                                                              TCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGA‑TTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAA  >  minE/2921496‑2921625

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