Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,967,179	T→C	T87A (ACT→GCT)	cheY ←	chemotaxis protein CheY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,967,179	0	T	C	100.0%	152.5 / NA	44	T87A (ACT→GCT)	cheY	chemotaxis protein CheY
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/24);  total (20/24)

CATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTA  >  NC_000913/1967043‑1967306                                                                                                                                         |                                                                                                                                cATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCAc                                                                                                                              >  1:862787/1‑140 (MQ=255)    ccTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCAt                                                                                                                           <  2:225006/140‑1 (MQ=255)     cTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCAtt                                                                                                                          >  1:27023/1‑140 (MQ=255)        acaTGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAac                                                                                                                       >  1:798402/1‑140 (MQ=255)            gCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTg                                                                                                                   >  2:965841/1‑140 (MQ=255)              ccAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGc                                                                                                                 >  1:551525/1‑140 (MQ=255)              ccAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGc                                                                                                                 <  1:763513/140‑1 (MQ=255)                   ttCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGc                                                                                                            <  2:63918/140‑1 (MQ=255)                       cAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGAc                                                                                                        <  2:165063/140‑1 (MQ=255)                            aTTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATcgc                                                                                                   >  1:414015/1‑140 (MQ=255)                                     aGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGc                                                                                          <  2:798402/140‑1 (MQ=255)                                               ccAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGtt                                                                                <  2:27023/140‑1 (MQ=255)                                                   cgacgcgGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGc                                                                                          <  2:620740/123‑1 (MQ=255)                                                   cGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGc                                                                                          >  1:736970/1‑126 (MQ=255)                                                   cGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGc                                                                                          >  1:620740/1‑126 (MQ=255)                                                   cGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGc                                                                                          <  2:736970/126‑1 (MQ=255)                                                    gTCGCGGCGGTAAATGGCTTCACCACATAGCCACTTGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAg                                                                           <  2:364352/140‑1 (MQ=255)                                                    gTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAg                                                                           <  2:551525/140‑1 (MQ=255)                                                     tCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAg                                                                           <  1:739961/139‑1 (MQ=255)                                                     tCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAg                                                                           >  2:739961/1‑139 (MQ=255)                                                              tAAATGGCTTCACCTCATAGCCACTGGCCCTCGCTTGCGCCGCAGCAATGAGGCTCTCTATCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAgg                                                                 <  2:80869/140‑1 (MQ=255)                                                              tAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCGCCATTAACACGGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAgg                                                                 <  1:89344/140‑1 (MQ=255)                                                                  tGGCTTCACCGCATAGCCACTGGCCCCCGGTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGccat                                                             <  1:447975/140‑1 (MQ=255)                                                                  tGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGccat                                                             <  1:784215/140‑1 (MQ=255)                                                                      ttCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGTCATCGCGCCATCCGCAAGAATTGTTGTCAGCAATTCCAGGccatccat                                                         >  2:394712/1‑140 (MQ=255)                                                                      ttCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGccatccat                                                         >  2:829862/1‑140 (MQ=255)                                                                       tCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCata                                                        <  1:965841/140‑1 (MQ=255)                                                                        caccacATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCatat                                                       >  1:785149/1‑140 (MQ=255)                                                                               tAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCAt                                                <  2:157145/140‑1 (MQ=255)                                                                                       ggCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGt                                        <  2:129464/140‑1 (MQ=255)                                                                                       ggCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGt                                        >  2:581478/1‑140 (MQ=255)                                                                                        gCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCGGCAATTCCAGGCCAGCCATATGGGGCATGTTCCAGTc                                       >  1:259497/1‑140 (MQ=255)                                                                                           cccGCTTGCGCCGCAGCGATGATGTTCTCTTCCTTCGCTACTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGga                                    <  2:259497/140‑1 (MQ=255)                                                                                            ccGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGgag                                   <  2:862787/140‑1 (MQ=255)                                                                                            ccGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGgag                                   <  1:37093/140‑1 (MQ=255)                                                                                             cGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGTTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGgaga                                  <  2:414015/140‑1 (MQ=255)                                                                                                 tGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCGTGTTCCAGTCGGAGATAAc                              <  2:824834/140‑1 (MQ=255)                                                                                                        cagcaATGATGTTCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCa                       >  2:133517/1‑140 (MQ=255)                                                                                                         agcaATGATGTTCTTTTCCTTCGCTTCTGCAGCCATCATTAACACTGGTAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCAt                      <  2:319735/140‑1 (MQ=255)                                                                                                                   ttCTCTTTCTTCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCgcctgc            >  2:678321/1‑140 (MQ=255)                                                                                                                       ttttcttCGCTTCTGCAGCCCCCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCctgcctgc        <  2:216728/139‑1 (MQ=255)                                                                                                                        tttcttCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCa       >  1:173541/1‑140 (MQ=255)                                                                                                                        tttcttCGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCa       >  1:193378/1‑140 (MQ=255)                                                                                                                             tcGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCgcctgcct          <  2:86119/132‑1 (MQ=255)                                                                                                                             tcGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCgcctgcct          >  1:86119/1‑132 (MQ=255)                                                                                                                             tcGCTTCTGCAGCCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTa  >  2:25533/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                CATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTA  >  NC_000913/1967043‑1967306

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