Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,967,290	C→T	G50S (GGT→AGT)	cheY ←	chemotaxis regulator transmitting signal to flagellar motor component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,967,290	0	C	T	100.0%	164.8 / NA	44	G50S (GGT→AGT)	cheY	chemotaxis regulator transmitting signal to flagellar motor component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/27);  total (17/27)

GCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCAT  >  NC_000913/1967142‑1967437                                                                                                                                                     |                                                                                                                                                    gccgcAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTgc                                                                                                                                                   <  1:275802/151‑1 (MQ=255)  ccgcAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTgcc                                                                                                                                                  <  2:145375/151‑1 (MQ=255)          atgatgTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAActgcctgcctgc                                                                                                                                           >  1:24718/1‑150 (MQ=255)          atgatgTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCa                                                                                                                                          <  2:362118/151‑1 (MQ=255)                      ttcttcGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTgaga                                                                                                                               <  2:191049/150‑1 (MQ=255)                      ttcttcGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTga                                                                                                                                 <  1:583609/148‑1 (MQ=255)                      ttcttcGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTga                                                                                                                                 <  2:499972/148‑1 (MQ=255)                          tcGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCg                                                                                                                            <  1:507619/149‑1 (MQ=255)                          tcGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGt                                                                                                                           <  1:332546/150‑1 (MQ=255)                           cGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTc                                                                                                                          <  1:559298/150‑1 (MQ=255)                            gCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCg                                                                                                                         >  2:287576/1‑150 (MQ=255)                              ttCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGAc                                                                                                                       >  1:377196/1‑150 (MQ=255)                               tCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACg                                                                                                                      <  2:371624/150‑1 (MQ=255)                                 tGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGcc                                                                                                                    >  1:237570/1‑150 (MQ=255)                                 tGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGcc                                                                                                                    >  2:168573/1‑150 (MQ=255)                                  gCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCa                                                                                                                   <  2:93493/150‑1 (MQ=255)                                   cAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCGGCCTGCAACTTATTGAGAGCGTCGACGCCATc                                                                                                                 <  1:28972/151‑1 (MQ=255)                                       cACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTcc                                                                                                             <  1:459824/151‑1 (MQ=255)                                                    ggCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACat                                                                                                >  2:73293/1‑151 (MQ=255)                                                    ggCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACa                                                                                                 <  1:230393/150‑1 (MQ=255)                                                                     gcCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGctct                                                                               >  1:72696/1‑151 (MQ=255)                                                                       cATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCttt                                                                             >  1:66724/1‑151 (MQ=255)                                                                            gCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTcagca                                                                        <  2:236910/151‑1 (MQ=255)                                                                              aCGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAgg                                                                      >  2:533583/1‑151 (MQ=255)                                                                                gAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGtt                                                                    >  1:445561/1‑151 (MQ=255)                                                                                     gTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGc                                                                <  1:666849/150‑1 (MQ=255)                                                                                            gCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCg                                                        <  1:473639/151‑1 (MQ=255)                                                                                             cAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGt                                                       <  2:666078/151‑1 (MQ=255)                                                                                              aaTTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTc                                                      >  2:24636/1‑151 (MQ=255)                                                                                                    aGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCAtgg                                                >  1:87134/1‑151 (MQ=255)                                                                                                       ccatccatATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCAtggtgg                                             >  2:643286/1‑151 (MQ=255)                                                                                                          tccatATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGaa                                           >  2:397186/1‑150 (MQ=255)                                                                                                                  gggCATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATcc                                  <  2:24717/151‑1 (MQ=255)                                                                                                                     cATGTTCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCcaca                               <  1:299408/151‑1 (MQ=255)                                                                                                                         ttCCAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAAcc                            >  1:401885/1‑150 (MQ=255)                                                                                                                           ccAGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCaa                          <  2:435769/150‑1 (MQ=255)                                                                                                                             aGTCGGAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCaaaaa                       <  2:56221/151‑1 (MQ=255)                                                                                                                                 ggAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTa                   <  2:445566/151‑1 (MQ=255)                                                                                                                                 ggAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTa                   <  1:643276/151‑1 (MQ=255)                                                                                                                                 ggAGATAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTa                   <  1:24637/151‑1 (MQ=255)                                                                                                                                   agaTAACAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTa                   >  2:514132/1‑149 (MQ=255)                                                                                                                                       aaCAAATCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCt             >  1:178289/1‑151 (MQ=255)                                                                                                                                           aaTCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTAt         <  1:513881/151‑1 (MQ=255)                                                                                                                                            aTCCATAACTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATc        <  2:592781/151‑1 (MQ=255)                                                                                                                                                  aaCTGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCAt  >  2:73354/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTACGCACTATGCGTCGCATGGTGGAAAAGTCATCCACAACCAAAAATTTAAGTTCTTTATCCGCCAT  >  NC_000913/1967142‑1967437

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑