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%	179.3 / NA	56	*138* (TAG→TAA)	fliN	flagellar motor switching and energizing component
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/28);  total (28/28)

CGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAGTGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTA  >  NC_000913/2021365‑2021620                                                                                                                                         |                                                                                                                        tgggcgAACCACTGGTTATTCTGATCAACAGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAAtg                                                                                                                       <  2:186092/138‑1 (MQ=255) cgggcgAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAAtg                                                                                                                       >  1:1556872/1‑139 (MQ=255)   ggcgAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAAtgat                                                                                                                     <  2:918364/139‑1 (MQ=255)         ccACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATaa                                                                                                               <  1:1251741/139‑1 (MQ=255)                 tatTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTa                                                                                                       <  2:1556872/139‑1 (MQ=255)                   ttCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACt                                                                                                     >  2:614089/1‑139 (MQ=255)                        aTCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCa                                                                                                <  1:614089/139‑1 (MQ=255)                         tCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCG‑CTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAAt                                                                                              >  1:1911024/1‑139 (MQ=255)                             cGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCt                                                                                                        >  1:1708181/1‑126 (MQ=255)                             cGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAAAATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCt                                                                                                        <  2:1708181/126‑1 (MQ=255)                                               gggCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTctgctg                                                                         >  2:1003541/1‑139 (MQ=255)                                                   gAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGcc                                                                     <  1:705109/139‑1 (MQ=255)                                                    aaGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCG‑CTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCAc                                                                   <  2:1911024/139‑1 (MQ=255)                                                     aGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCa                                                                                                            >  2:62681/1‑98 (MQ=255)                                                     aGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCa                                                                                                            <  1:62681/98‑1 (MQ=255)                                                     aGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCAc                                                                   >  2:1034064/1‑139 (MQ=255)                                                                               gCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTctg                                                                            <  1:287217/104‑1 (MQ=255)                                                                               gCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTctg                                                                            >  2:287217/1‑104 (MQ=255)                                                                                 gTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGc                                       <  1:1393053/139‑1 (MQ=255)                                                                                  tGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGTGCCGCCTGAGCCGTTAAtg                                                                                                                       <  1:672856/58‑1 (MQ=255)                                                                                   gCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAAt                                                                                                                 <  1:2032588/63‑1 (MQ=255)                                                                                   gCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAAt                                                                                                                 >  2:2032588/1‑63 (MQ=255)                                                                                   gCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCa                                     >  1:1034649/1‑139 (MQ=255)                                                                                    cGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCcgcg                                                                                      >  1:158575/1‑89 (MQ=255)                                                                                    cGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCcgcg                                                                                      <  2:158575/89‑1 (MQ=255)                                                                                     ggATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCActgc                                                                <  1:559124/110‑1 (MQ=255)                                                                                     ggATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCActgc                                                                >  2:559124/1‑110 (MQ=255)                                                                                                 tCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCAc                                                                   >  1:2033619/1‑95 (MQ=255)                                                                                                 tCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCAc                                                                   <  2:2033619/95‑1 (MQ=255)                                                                                                    ttACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGcc                             <  2:726074/130‑1 (MQ=255)                                                                                                    ttACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGcc                             >  1:726074/1‑130 (MQ=255)                                                                                                      aCTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAgcggcg                                                 <  2:1474362/108‑1 (MQ=255)                                                                                                      aCTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTATCTGCTGCGCCACTGCTGCAGGTGAgcggcg                                                 >  1:1474362/1‑108 (MQ=255)                                                                                                       cTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCattatt                                >  2:1990182/1‑124 (MQ=255)                                                                                                       cTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCattatt                                <  1:1990182/124‑1 (MQ=255)                                                                                                          cGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCatta                                  <  1:440080/119‑1 (MQ=255)                                                                                                          cGTCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCatta                                  >  2:440080/1‑119 (MQ=255)                                                                                                            tCTGAGCGAATGCGCCGCCTGAGCCGTTAAtgatga                                                                                                                   >  1:1019490/1‑36 (MQ=255)                                                                                                            tCTGAGCGAATGCGCCGCCTGAGCCGTTAAtgatga                                                                                                                   <  2:1019490/36‑1 (MQ=255)                                                                                                            tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTa                                                                                                       >  1:2031937/1‑48 (MQ=255)                                                                                                            tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTa                                                                                                       <  2:2031937/48‑1 (MQ=255)                                                                                                            tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCcgcg                                                                                      >  2:1748434/1‑65 (MQ=255)                                                                                                            tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCcgcg                                                                                      <  1:1748434/65‑1 (MQ=255)                                                                                                            tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCtgctgc            <  1:1003541/139‑1 (MQ=255)                                                                                                            tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCtgctgc            >  1:627736/1‑139 (MQ=255)                                                                                                            tCTGAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCtgctgc            >  2:1992917/1‑139 (MQ=255)                                                                                                               gAGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCctg         >  2:1419667/1‑139 (MQ=255)                                                                                                                aGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAgcggc                                                  >  1:507585/1‑97 (MQ=255)                                                                                                                aGCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAgcggc                                                  <  2:507585/97‑1 (MQ=255)                                                                                                                 gCGAATGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCctggc       <  1:1992917/139‑1 (MQ=255)                                                                                                                    aaTGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTa                                                                                                       >  1:905698/1‑40 (MQ=255)                                                                                                                    aaTGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTa                                                                                                       <  2:905698/40‑1 (MQ=255)                                                                                                                      tGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACt                                                                                                     <  1:1722161/40‑1 (MQ=255)                                                                                                                      tGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACt                                                                                                     >  2:1722161/1‑40 (MQ=255)                                                                                                                      tGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTa  >  1:307416/1‑139 (MQ=255)                                                                                                                      tGCGCCGCCTGAGCCGTTAATGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCCGCCTGGCTGGTa  >  1:931607/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                        CGGGCGAACCACTGGATATTCTGATCAACGGTTATTTAATCGCCCAGGGCGAAGTGGTGGTCGTTGCCGATAAATATGGCGTGCGGATCACCGATATCATTACTCCGTCTGAGCGAATGCGCCGCCTGAGCCGTTAGTGATGAATAACCACGCTACTGTGCAATCTTCCGCGCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTA  >  NC_000913/2021365‑2021620

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