Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,021,501	G→A	*138* (TAG→TAA)	fliN →	flagellar motor switching and energizing component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,021,501	0	G	A	100.0%	186.2 / NA	57	*138* (TAG→TAA)	fliN	flagellar motor switching and energizing component
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/27);  total (30/27)

CGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAGTGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGAT  >  NC_000913/2021369‑2021633                                                                                                                                     |                                                                                                                                     cgAACTACTGGATATTCTGATCATCGGTTAGTTAATCGCCTAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATAACCGATATCATTACTCCGTCTGAGCGAATGCGACGCCTGAGCCGTTAAtgatga                                                                                                                                >  2:1600303/1‑139 (MQ=255) cgAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAAtgatga                                                                                                                                <  1:139687/139‑1 (MQ=255)                     tCAACGGTTATTTAATCGCCCAGGGCGAAGTTGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAAt                                                                                                                              <  2:149716/121‑1 (MQ=255)                     tCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAAt                                                                                                                              >  1:149716/1‑121 (MQ=255)                         cGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCtt                                                                                                        >  2:1153422/1‑139 (MQ=255)                         cGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCtt                                                                                                        >  1:114461/1‑139 (MQ=255)                         cGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCtt                                                                                                        >  2:141616/1‑139 (MQ=255)                         cGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCtt                                                                                                        >  2:606046/1‑139 (MQ=255)                                      gCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTc                                                                                           >  2:805613/1‑139 (MQ=255)                                      gCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTc                                                                                           >  1:365565/1‑139 (MQ=255)                                               gAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGcc                                                                                  <  1:1600303/139‑1 (MQ=255)                                                 aGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCAc                                                                                <  1:141616/139‑1 (MQ=255)                                                      tggtCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCa                                                                                 <  2:27657/133‑1 (MQ=255)                                                      tggtCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCa                                                                                 >  1:27657/1‑133 (MQ=255)                                                      tggtCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCActgctg                                                                           >  1:1364364/1‑139 (MQ=255)                                                                   tAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAgcggcg                                                              <  1:1153422/139‑1 (MQ=255)                                                                            cGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCg                                                     >  1:1293694/1‑139 (MQ=255)                                                                             gTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGc                                                    <  2:1293694/139‑1 (MQ=255)                                                                             gTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCCGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGc                                                    <  1:805613/139‑1 (MQ=255)                                                                              tGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACg                                                                                                                       >  1:214403/1‑71 (MQ=255)                                                                              tGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACg                                                                                                                       <  2:214403/71‑1 (MQ=255)                                                                                cGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCa                                                                                 >  1:296395/1‑107 (MQ=255)                                                                                cGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCa                                                                                 <  2:296395/107‑1 (MQ=255)                                                                                  gATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCatta                                               <  1:606046/139‑1 (MQ=255)                                                                                  gATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCatta                                               <  2:114461/139‑1 (MQ=255)                                                                                             tCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTa                                                                                                                    <  2:842009/59‑1 (MQ=255)                                                                                             tCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTa                                                                                                                    >  1:842009/1‑59 (MQ=255)                                                                                                  aCTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGcc                                                                                  >  1:296186/1‑88 (MQ=255)                                                                                                  aCTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGcc                                                                                  <  2:296186/88‑1 (MQ=255)                                                                                                      cGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAAt                                                                                                                              <  2:915770/40‑1 (MQ=255)                                                                                                      cGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAAt                                                                                                                              >  1:915770/1‑40 (MQ=255)                                                                                                      cGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCt                                  <  1:600432/132‑1 (MQ=255)                                                                                                      cGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCt                                  >  2:600432/1‑132 (MQ=255)                                                                                                      cGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCgct                              <  2:369255/136‑1 (MQ=255)                                                                                                      cGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCgct                              >  1:369255/1‑136 (MQ=255)                                                                                                        tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTctgc                                                                                        >  1:1285123/1‑76 (MQ=255)                                                                                                        tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTctgc                                                                                        <  2:1285123/76‑1 (MQ=255)                                                                                                        tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAgcgg                                                                >  1:1222090/1‑100 (MQ=255)                                                                                                        tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAgcgg                                                                <  2:1222090/100‑1 (MQ=255)                                                                                                        tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGc                                                    >  2:496777/1‑112 (MQ=255)                                                                                                        tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGc                                                    <  1:496777/112‑1 (MQ=255)                                                                                                        tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCtgctgc                         >  2:1446246/1‑139 (MQ=255)                                                                                                        tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCtgctgc                         >  2:1464746/1‑139 (MQ=255)                                                                                                        tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCtgctgc                         >  1:746436/1‑139 (MQ=255)                                                                                                        tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCtgctgc                         >  2:771547/1‑139 (MQ=255)                                                                                                        tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCtgctgc                         >  2:663567/1‑139 (MQ=255)                                                                                                                aaTGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGt                                                                      <  2:629960/86‑1 (MQ=255)                                                                                                                aaTGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGt                                                                      >  1:629960/1‑86 (MQ=255)                                                                                                                 aTGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCcgcg                                                                                                   <  1:155636/56‑1 (MQ=255)                                                                                                                 aTGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCcgcg                                                                                                   >  2:155636/1‑56 (MQ=255)                                                                                                                  tGCGCCGCCTGAGCCGTTAATGATGAATAACCACGc                                                                                                                      <  1:1204350/36‑1 (MQ=255)                                                                                                                  tGCGCCGCCTGAGCCGTTAATGATGAATAACCACGc                                                                                                                      >  2:1204350/1‑36 (MQ=255)                                                                                                                  tGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAgcggcg                                                              <  1:167194/92‑1 (MQ=255)                                                                                                                  tGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAgcggcg                                                              >  2:167194/1‑92 (MQ=255)                                                                                                                     gccgccTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTaaaa            <  1:663567/139‑1 (MQ=255)                                                                                                                          cTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGtt       <  1:1464746/139‑1 (MQ=255)                                                                                                                               ccGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGAt  <  2:1364364/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                     CGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAGTGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGAT  >  NC_000913/2021369‑2021633

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