Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,974,881	T→C	E148E (GAA→GAG)	cheA ←	chemotaxis protein CheA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,974,881	0	T	C	100.0%	111.3 / NA	32	E148E (GAA→GAG)	cheA	chemotaxis protein CheA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/18);  total (14/18)

TTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCG  >  NC_000913/1974755‑1975016                                                                                                                               |                                                                                                                                        tttCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCaa                                                                                                                            <  1:109893/140‑1 (MQ=255)  ttCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAAc                                                                                                                           <  2:696840/140‑1 (MQ=255)   tCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAAcc                                                                                                                          <  2:23843/140‑1 (MQ=255)       cacGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCacac                                                                                                                      <  2:144210/140‑1 (MQ=255)           tCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTaa                                                                                                                  >  1:796419/1‑140 (MQ=255)                    gTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCAc                                                                                                         >  2:570080/1‑140 (MQ=255)                    gTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCAc                                                                                                         >  2:29562/1‑140 (MQ=255)                        tCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCg                                                                                                     >  1:655744/1‑140 (MQ=255)                        tCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCg                                                                                                     >  1:7312/1‑140 (MQ=255)                           gATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGAt                                                                                                  <  2:546739/140‑1 (MQ=255)                              gTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGc                                                                                               >  1:791984/1‑140 (MQ=255)                               tCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCg                                                                                              >  2:661454/1‑140 (MQ=255)                                 ccAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGtt                                                                                            <  1:135576/140‑1 (MQ=255)                                  cAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGttt                                                                                           >  2:150680/1‑140 (MQ=255)                                         tcttcCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACATTTGGCAACCACACTTCATCGGGTCTCTGCGCATGGCGTATCGCCttt                                                                                    >  2:488128/1‑140 (MQ=255)                                               agcagGTCGACTTCCGCGGCCTCCACGCGGGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTc                                                                              <  1:488128/140‑1 (MQ=255)                                                          ttCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGtt                                                                   <  1:570080/140‑1 (MQ=255)                                                           tCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTg                                                                  >  2:224544/1‑140 (MQ=255)                                                                ggCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCa                                                             <  2:959398/140‑1 (MQ=255)                                                                      cAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCt                                                       >  2:659106/1‑140 (MQ=255)                                                                      cAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCt                                                       >  1:127388/1‑140 (MQ=255)                                                                        ggCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTgg                                                     <  1:768591/140‑1 (MQ=255)                                                                        ggCGCGAAAGGATAACTCGTCGGGGCGACTGACTGCGTCTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTgg                                                     <  1:499388/140‑1 (MQ=255)                                                                                       ttCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAg                                      <  1:150680/140‑1 (MQ=255)                                                                                            cgcgGCGACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTggc                                 >  1:581170/1‑140 (MQ=255)                                                                                             gcggcgACTGACTGCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTggcg                                <  2:581170/140‑1 (MQ=255)                                                                                                      gactgCGACTCTGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGc                       <  2:501757/140‑1 (MQ=255)                                                                                                              ctctGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGa               <  2:127388/140‑1 (MQ=255)                                                                                                              ctctGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGa               <  1:772793/140‑1 (MQ=255)                                                                                                              ctctGCTCATCTTGCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGa               <  1:29562/140‑1 (MQ=255)                                                                                                                      aTCTTGTGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTtat       >  2:840682/1‑140 (MQ=255)                                                                                                                           gCGGCTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCg  <  2:707777/140‑1 (MQ=255)                                                                                                                               |                                                                                                                                        TTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTGGCAGATATAATCGAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCG  >  NC_000913/1974755‑1975016

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑