Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,971,513	Δ1 bp	coding (1179/1662 nt)	tar ←	methyl‑accepting chemotaxis protein Tar

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,971,510	0	A	.	100.0%	222.5 / NA	52	coding (1182/1662 nt)	tar	methyl‑accepting chemotaxis protein Tar
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (25/27);  total (25/27)

TCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGC  >  NC_000913/1971371‑1971645                                                                                                                                            |                                                                                                                                        tCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCa                                                                                                                                         >  2:772112/1‑140 (MQ=255)     tGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCaaacc                                                                                                                                     >  1:769061/1‑138 (MQ=255)     tGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCaaacc                                                                                                                                     >  1:517866/1‑138 (MQ=255)      gTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCaaacca                                                                                                                                    <  1:1001343/140‑4 (MQ=255)      gTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCaaacca                                                                                                                                    <  1:457157/140‑4 (MQ=255)         tCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACgg                                                                                                                                >  1:346868/1‑140 (MQ=255)          ccccGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGc                                                                                                                               >  2:187048/1‑140 (MQ=255)           cccGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTCCGCACTTCCCCCGCCACCACGGC‑AAACCACGGcc                                                                                                                              <  2:211059/140‑1 (MQ=255)           cccGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGcc                                                                                                                              <  2:637503/140‑1 (MQ=255)           cccGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGcc                                                                                                                              <  2:525502/140‑1 (MQ=255)               gcgcTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTCTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGt                                                                                                                          >  1:830811/1‑140 (MQ=255)                 gcTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTc                                                                                                                        >  1:1079990/1‑140 (MQ=255)                            gCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTAAATGGGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGAAAGATTCCGCACTTCCCCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCAcgcg                                                                                                             <  2:535225/140‑1 (MQ=255)                                  aaCCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCtt                                                                                                       >  2:187355/1‑140 (MQ=255)                                     cGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCaa                                                                                                    <  2:769061/140‑1 (MQ=255)                                               gcgTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCAtt                                                                                          >  2:306753/1‑140 (MQ=255)                                                    agaCGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAgcg                                                                                     >  2:803482/1‑140 (MQ=255)                                                     gaCGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCCCTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAgcgc                                                                                    <  2:1055028/140‑1 (MQ=255)                                                        ggAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAg                                                                                 <  2:574189/140‑1 (MQ=255)                                                          aGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGa                                                                               >  1:1015408/1‑140 (MQ=255)                                                           gTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGat                                                                              >  2:548566/1‑140 (MQ=255)                                                            tCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGata                                                                             <  1:187355/140‑1 (MQ=255)                                                            tCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGata                                                                             <  2:1015408/140‑1 (MQ=255)                                                               tCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTa                                                                          >  1:108183/1‑140 (MQ=255)                                                                cAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCCCCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAg                                                                         <  1:941611/140‑1 (MQ=255)                                                                 aaTGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCTCTTCGCCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGt                                                                        <  1:1053206/140‑1 (MQ=255)                                                                    gAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTg                                                                     <  1:803482/140‑1 (MQ=255)                                                                    gAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTg                                                                     <  1:366498/140‑1 (MQ=255)                                                                    gAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTg                                                                     <  2:346868/140‑1 (MQ=255)                                                                    gAGGGCTTTGATCTCTTTTGCCCCCTGGGCGCAGGGACTGGCAAGTTTACGCCCTTCCCCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTg                                                                     <  2:443386/140‑1 (MQ=255)                                                                      gggCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGa                                                                                  <  2:507524/125‑1 (MQ=255)                                                                      gggCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGa                                                                                  >  1:507524/1‑125 (MQ=255)                                                                      gggCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGa                                                                   >  2:861794/1‑140 (MQ=255)                                                                       ggCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGaa                                                                  >  1:526186/1‑140 (MQ=255)                                                                        gCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAg                                                                 >  2:991912/1‑140 (MQ=255)                                                                        gCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAg                                                                 >  1:656021/1‑140 (MQ=255)                                                                           ttGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCa                                                              <  2:606484/140‑1 (MQ=255)                                                                              ttctcTTTTGCCGCCGGGGCGCGGCGACGGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATa                                                           <  2:517866/139‑1 (MQ=255)                                                                                   ttttGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTc                                                      >  2:1157/1‑140 (MQ=255)                                                                                      tGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGAt                                                   >  2:96677/1‑140 (MQ=255)                                                                                        ccgccTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATaa                                                 <  1:772112/140‑1 (MQ=255)                                                                                          gccTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACg                                               <  1:548566/140‑1 (MQ=255)                                                                                               ggCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGAt                                                   <  2:699375/131‑1 (MQ=255)                                                                                               ggCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGAt                                                   >  1:699375/1‑131 (MQ=255)                                                                                                       gACTGGCAAGATTACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCgg                                  <  1:96677/140‑1 (MQ=255)                                                                                                                   tACGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCg                      >  1:650313/1‑140 (MQ=255)                                                                                                                     cGCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGaa                    >  1:1051288/1‑140 (MQ=255)                                                                                                                      gCACTTCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAAc                   >  2:146545/1‑140 (MQ=255)                                                                                                                          ttCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATc               >  2:216196/1‑140 (MQ=255)                                                                                                                           tCACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCg              >  2:844276/1‑140 (MQ=255)                                                                                                                            cACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCgg             <  2:108183/140‑1 (MQ=255)                                                                                                                            cACCCGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCgg             <  1:861794/140‑1 (MQ=255)                                                                                                                               ccGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGa          >  1:593989/1‑140 (MQ=255)                                                                                                                                cGCCACCACGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGAt         <  1:146545/140‑1 (MQ=255)                                                                                                                                    accacGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCa     >  2:502442/1‑140 (MQ=255)                                                                                                                                       acGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGc  <  1:216196/140‑1 (MQ=255)                                                                                                                                       acGGC‑AAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGc  <  2:650313/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                        TCATTGTTTCCCCGGCGCTTTCGACCAGCACCGAACCGGTATCAACGCGTGAGACGGAGTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATTTTCTTCGAACTATCGGCGATCTCATGC  >  NC_000913/1971371‑1971645

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