Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,973,522	A→G	G601G (GGT→GGC)	cheA ←	chemotaxis protein CheA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,973,522	0	A	G	100.0%	60.7 / NA	19	G601G (GGT→GGC)	cheA	chemotaxis protein CheA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/6);  total (13/6)

TTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGACCAATTAATTG‑ATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATT  >  NC_000913/1973387‑1973655                                                                                                                                         |                                                                                                                                       ttCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCaa                                                                                                                                    >  1:418761/1‑140 (MQ=255)    gcgGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCaatta                                                                                                                                 <  2:598604/140‑1 (MQ=255)    gcgGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCaatta                                                                                                                                 <  2:1036135/140‑1 (MQ=255)       gTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCaattaatt                                                                                                                              >  2:573367/1‑140 (MQ=255)         ttATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCAATTAATTG‑a                                                                                                                           <  1:172883/140‑1 (MQ=255)           aTCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCAATTAATTG‑ATc                                                                                                                         >  2:458725/1‑140 (MQ=255)              gCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCAATTAATTG‑ATccac                                                                                                                      <  1:573367/140‑1 (MQ=255)                cTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCAATTAATTG‑ATccacca                                                                                                                    <  2:418761/140‑1 (MQ=255)                  gCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCAATTAATTG‑ATCCACCAGc                                                                                                                  >  1:738198/1‑140 (MQ=255)                  gCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCAATTAATTG‑ATCCACCAGc                                                                                                                  >  2:715838/1‑140 (MQ=255)                                       tCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCAATTAATTG‑ATCCACCAGCAAGGCGTAGCGGCGACCGCCa                                                                                             >  1:46802/1‑140 (MQ=255)                                               aCGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACACCCTGGTGTTGGCAAATTAATTG‑ATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTaa                                                                                     <  1:447170/140‑1 (MQ=255)                                                     ccGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCAATTAATTG‑ATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATcac                                                                               >  2:884407/1‑140 (MQ=255)                                                        tCGCAAAGAATGGTCGCAGAAGAAACGCCGGGGACTTTGCGATATTTTTTTACAAGTTTTTT‑AACCACACCCTGGTGTTGGCAAATTAATTG‑ATCCACCAGCAGGGCGTAGCGGCGACCGCCACTTTGTAAGATcaccac                                                                            <  2:738198/140‑1 (MQ=255)                                                               gAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCAATTAATTG‑ATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTccc                                                                     >  1:321025/1‑140 (MQ=255)                                                                   ggTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCAATTAATTG‑ATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTggg                                                                 >  1:1031592/1‑140 (MQ=255)                                                                      cGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCAATTAATTG‑ATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTgg                                                              >  2:818242/1‑140 (MQ=255)                                                                                          ttttGCGATAGTTACTTT‑AAGGTTTTTAAACCACACCCTGGTG‑TGGCCAATTTATTGAATCCCCCAGCAGGGCGGGGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCcgcg                                          <  2:494378/139‑1 (MQ=255)                                                                                                     ttACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCAATTAATTG‑ATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAacac                               >  2:106127/1‑140 (MQ=255)                                                                                                       aCTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGGCCAATTAATTG‑ATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACtt                             >  2:950622/1‑140 (MQ=255)                                                                                                                   ttt‑AACCACACCCGGGTGTGGCCAAATTAATTG‑ACCCCCCAGCAGGGCGTGGGGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTcga                 <  2:46802/140‑1 (MQ=255)                                                                                                                      tacccACAACCTGGTGTTGCCCATTTAATTG‑ACCCCCCAGCAGGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTcgacg               <  1:162669/137‑1 (MQ=255)                                                                                                                       acccACACCCTGGTGTCGGCAAATTAATTG‑ACCCACCAGCAGGGCGTGGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTcgacga              <  2:1031592/138‑1 (MQ=255)                                                                                                                         ccACACCCTGGTGTTGGCAAATTAATTG‑ACCCCCCAGCAGGGCGGGGGGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATg            <  1:715838/140‑1 (MQ=255)                                                                                                                         ccACACCCTGGTGTTGGCAAATTAATTG‑ACCCCCCAGCAAGGCGTGGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATg            <  1:884407/140‑1 (MQ=255)                                                                                                                                   gtgtTGGCCAATTAATTG‑ATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATAtt  >  2:331602/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                       TTCGCGGTTTATCGCCTGCAAGGCGGAGACATCAACAATCAGTGCCACGCTGCCGTCGCCAAGAATGGTCGCAGCAGAAATGCCGGGGACTTTGCGATAGTTACTTTCAAGGTTTTT‑AACCACAACCTGGTGTTGACCAATTAATTG‑ATCCACCAGCAAGGCGTAGCGGCGACCGCCACTTTGTAAGATCACCACAATTCCCTGGGTGGCTTCGGTTTTCGCGCCCGCGACGTTGAACACTTTCCACAGTTCGACGATGGGCAGATATT  >  NC_000913/1973387‑1973655

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