Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,523,658	C→T	Q153* (CAA→TAA)	flk →	predicted flagella assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,523,658	0	C	T	100.0%	121.1 / NA	41	Q153* (CAA→TAA)	flk	predicted flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/20);  total (21/20)

TCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGCAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAAC  >  minE/1523596‑1523721                                                               |                                                                tctcACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAAc                                                              >  1:137300/1‑66 (MQ=255)     aCCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCAcag                                                         >  1:1342609/1‑67 (MQ=255)        ggAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                                   >  1:2338691/1‑70 (MQ=255)         gAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAgcgc                                                   >  1:353904/1‑69 (MQ=255)               ttAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGcgac                                           >  1:1728354/1‑71 (MQ=255)                 aaaaaCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACcaca                                                          >  1:1195722/1‑54 (MQ=255)                  aaaaCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGcgcg                                                  >  1:1212847/1‑61 (MQ=255)                  aaaaCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGcgaccga                                        >  1:1721026/1‑71 (MQ=255)                    aaCCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCgaccgacc                                      <  1:1081211/71‑1 (MQ=255)                         ttCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGcgaccgac                                       >  1:2897596/1‑65 (MQ=255)                         ttCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGcgaccgac                                       >  1:734319/1‑65 (MQ=255)                           ctcACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCtt                                <  1:2942249/70‑1 (MQ=255)                            tcACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCtt                                <  1:1793158/69‑1 (MQ=255)                            tcACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTtta                              >  1:1837412/1‑71 (MQ=255)                             cACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGcg                                             <  1:311404/55‑1 (MQ=255)                              aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGcga                                            >  1:906562/1‑55 (MQ=255)                              aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTttatt                            <  1:2123507/71‑1 (MQ=255)                              aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTttatt                            <  1:565358/71‑1 (MQ=255)                              aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTttatt                            <  1:249617/71‑1 (MQ=255)                              aCCCTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTttat                             >  1:1914783/1‑70 (MQ=255)                                ccTGTTGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGcgaccga                                        <  1:1580500/57‑1 (MQ=255)                                 cTGTTGCAGGAAGGGAAGATGGTTATTACGTAACCACagcagc                                                     >  1:245807/1‑43 (MQ=38)                                    ttGCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGaccgaccg                                     <  1:2631268/56‑1 (MQ=255)                                      gCAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcggcg                     >  1:458139/1‑70 (MQ=255)                                       cAGGAAGGGAAGATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGa                   <  1:874977/71‑1 (MQ=255)                                        aGGAAGGGAAGATGGTTATTCC‑TAACCACAGCAGCGCGAGGCGACCGACCGTc                                   <  1:904012/53‑1 (MQ=38)                                           aagggaagATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGcac               <  1:1800988/71‑1 (MQ=255)                                            agggaagATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGcac               >  1:561679/1‑70 (MQ=255)                                               gaagATGGTTATTCCGTAACCACAGCAGCGCGAGGcgaccga                                        <  1:2699038/42‑1 (MQ=255)                                                aagATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGTGGAGCACAATgc          >  1:1178412/1‑71 (MQ=255)                                                aagATGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATgc          >  1:2730248/1‑71 (MQ=255)                                                    tGGTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCt       <  1:1943178/70‑1 (MQ=255)                                                     ggTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACcggcg                     >  1:2355925/1‑55 (MQ=255)                                                     ggTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGTCGGAGCACAATGCGCTa      >  1:1514500/1‑70 (MQ=255)                                                      gTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTaa     <  1:1741976/70‑1 (MQ=255)                                                      gTTATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTaa     <  1:1658462/70‑1 (MQ=255)                                                       ttATTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTaaaa   <  1:1161764/71‑1 (MQ=255)                                                         aTTCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAAc  <  1:204351/70‑1 (MQ=255)                                                           tCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAg                  >  1:1882588/1‑52 (MQ=255)                                                           tCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTaa     <  1:1686547/65‑1 (MQ=255)                                                           tCCGTAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTaa     <  1:9412/65‑1 (MQ=255)                                                               |                                                                TCTCACCGGAGCAATTAAAAACCATTCTCACCCTGTTGCAGGAAGGGAAGATGGTTATTCCGCAACCACAGCAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAAC  >  minE/1523596‑1523721

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