Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,976,988	C→T	G51D (GGC→GAC)	motA ←	motility protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,976,988	0	C	T	100.0%	161.5 / NA	48	G51D (GGC→GAC)	motA	motility protein A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/22);  total (26/22)

TCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTA  >  NC_000913/1976866‑1977121                                                                                                                           |                                                                                                                                      tCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATa                                                                                                                      >  1:628657/1‑140 (MQ=255)  cTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATaa                                                                                                                     <  2:951973/140‑1 (MQ=255)    gCCGCGATTCCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAAc                                                                                                                   >  1:622359/1‑140 (MQ=255)       gcgATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGAc                                                                                                                <  2:302401/140‑1 (MQ=255)         gATTCCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGAccc                                                                                                              <  2:622359/140‑1 (MQ=255)               gCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATccc                                                                                                        >  2:160745/1‑140 (MQ=255)                    cAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGAccc                                                                                                              <  1:342025/129‑1 (MQ=255)                    cAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGAccc                                                                                                              >  2:342025/1‑129 (MQ=255)                                         cagATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAg                                                                              <  2:948695/140‑1 (MQ=255)                                           gATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCg                                                                            <  1:653823/140‑1 (MQ=255)                                                      aTTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATagag                                                                 >  2:62563/1‑140 (MQ=255)                                                      aTTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATagag                                                                 >  1:274955/1‑140 (MQ=255)                                                      aTTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATagag                                                                 <  2:274955/140‑1 (MQ=255)                                                      aTTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATagag                                                                 >  1:474198/1‑140 (MQ=255)                                                        tGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATagag                                                                 >  2:41890/1‑138 (MQ=255)                                                        tGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATagag                                                                 <  1:41890/138‑1 (MQ=255)                                                        tGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTg                                                               >  1:152010/1‑140 (MQ=255)                                                        tGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTg                                                               >  2:608791/1‑140 (MQ=255)                                                        tGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTg                                                               >  2:94709/1‑140 (MQ=255)                                                        tGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTg                                                               >  1:939790/1‑140 (MQ=255)                                                        tGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTg                                                               >  1:152361/1‑140 (MQ=255)                                                        tGCTTTGGTGGATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTg                                                               >  2:8292/1‑140 (MQ=255)                                                             tGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCa                                                          >  2:602258/1‑140 (MQ=255)                                                               gtgtATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATaga                                                                  >  1:126507/1‑130 (MQ=255)                                                               gtgtATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATaga                                                                  <  2:126507/130‑1 (MQ=255)                                                               gtgtATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAg                                                        <  1:94709/140‑1 (MQ=255)                                                               gtgtATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAg                                                        >  2:623861/1‑140 (MQ=255)                                                               gtgtATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAg                                                        >  1:429028/1‑140 (MQ=255)                                                               gtgtATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAg                                                        >  2:304293/1‑140 (MQ=255)                                                                tgtATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTcc                                                           <  1:859766/136‑1 (MQ=255)                                                                tgtATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTcc                                                           >  2:859766/1‑136 (MQ=255)                                                                  tATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCt                                                     <  1:718582/140‑1 (MQ=255)                                                                  tATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCt                                                     >  2:718582/1‑140 (MQ=255)                                                                    tttGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTc                                                   >  1:969756/1‑140 (MQ=255)                                                                       ggAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCAc                                                <  2:257117/140‑1 (MQ=255)                                                                        gAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCAcc                                               <  2:593981/140‑1 (MQ=255)                                                                          gCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCt                                                             >  2:440890/1‑124 (MQ=255)                                                                          gCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCt                                                             <  1:440890/124‑1 (MQ=255)                                                                           cgacgaAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGt                                            <  1:602258/140‑1 (MQ=255)                                                                             acgaAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGtca                                          >  1:385416/1‑140 (MQ=255)                                                                                                     cAGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTAc                  <  1:62563/140‑1 (MQ=255)                                                                                                     cAGCGTCCCTTAGATCGTTTTTTTATTATTGCCCCCGATAAACGACCAACTCCCTCCACCGCCAATCATCACCAGTGTGGGGGTTTGATATATCGCTCAAAGGCTCCCCCGGGTCATCAAATAACCGCCGAAAACTGTAc                  <  2:208476/140‑1 (MQ=255)                                                                                                      aGCGTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTAcc                 >  2:575152/1‑140 (MQ=255)                                                                                                         gTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCgag              >  2:635781/1‑140 (MQ=255)                                                                                                         gTGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCAAAGGCTTCCACCGGTCATAAAATAACCGCCGAAAACTGTACCgag              <  1:160745/140‑1 (MQ=255)                                                                                                          tGCCTTTAATCGCTTTGTCATTATTGCCGACGATAAACGCCCCAATCCCTGCCCCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCgaga             <  1:200059/140‑1 (MQ=255)                                                                                                             cTTTAATCGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGaaca          <  2:152010/140‑1 (MQ=255)                                                                                                                    cGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCATGGCTTCCACCGGTCATAAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGt   <  2:474198/140‑1 (MQ=255)                                                                                                                    cGCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGt   <  2:152361/140‑1 (MQ=255)                                                                                                                     gCTTTGTCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTa  <  1:575152/140‑1 (MQ=255)                                                                                                                           |                                                                                                                                      TCTGCCGCGATTTCGCCATCAACCGATAAAGCAGAGCCAGCAGATCCATATACATTGCTTTGGTGTATTTGGAGCGACGAAACAGCAACGGCAGCGCCTTCAGCGTGCCTTTAATCGCTTTGCCATTATTGCCGACGATAAACGACCCAATCCCTGCACCGGCAATAATCACCAGTTCAGCGGGTTGATAGAGTGCTCCAAGGCTTCCACCGGTCATCAAATAACCGCCGAAAACTGTACCGAGAACAACCAGGTA  >  NC_000913/1976866‑1977121

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