Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,022,885	G→A	*90* (TAG→TAA)	fliQ →	flagellar biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,022,885	0	G	A	100.0%	178.5 / NA	55	*90* (TAG→TAA)	fliQ	flagellar biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/26);  total (29/26)

CGAAATGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAGCCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTACCGAAACGGGTAA  >  NC_000913/2022750‑2023019                                                                                                                                        |                                                                                                                                       cGAAATGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCg                                                                                                                                     >  2:210729/1‑139 (MQ=255) cGAAATGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCg                                                                                                                                     <  2:1074693/139‑1 (MQ=255)      tGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTa                                                                                                                                >  1:1685972/1‑139 (MQ=255)      tGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTa                                                                                                                                <  1:210729/139‑1 (MQ=255)         cGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGt                                                                                                                             <  1:201866/139‑1 (MQ=255)           cTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTg                                                                                                                           <  2:1685972/139‑1 (MQ=255)                  ttgttCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGa                                                                                                                    <  1:678026/136‑1 (MQ=255)                  ttATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGT‑GATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGAc                                                                                                                   <  2:1667392/139‑1 (MQ=255)                          aaaaTCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGaa                                                                                                            >  2:775305/1‑139 (MQ=255)                            aaTCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACa                                                                                                          <  2:1261464/139‑1 (MQ=255)                                          ttATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGt                                                                                            <  1:73992/139‑1 (MQ=255)                                               gCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAAc                                                                                       >  2:113631/1‑139 (MQ=255)                                                   ttattattGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTACAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGt                                                                                   >  1:1443198/1‑139 (MQ=255)                                                        attGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTc                                                                              >  1:115292/1‑139 (MQ=255)                                                         ttGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTACAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCt                                                                             <  2:1443198/139‑1 (MQ=255)                                                           gCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTgg                                                                           >  2:1471607/1‑139 (MQ=255)                                                                aCCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTg                                                                                                                           >  2:1758238/1‑86 (MQ=255)                                                                aCCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTg                                                                                                                           <  1:1758238/86‑1 (MQ=255)                                                                aCCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGtt                                                                      >  2:291518/1‑139 (MQ=255)                                                                    tGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                                 >  1:1408270/1‑108 (MQ=255)                                                                    tGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                                 <  2:1408270/108‑1 (MQ=255)                                                                    tGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTg                                                                  >  1:26742/1‑139 (MQ=255)                                                                        tGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGcgcg                                                              <  2:319469/139‑1 (MQ=255)                                                                          cTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAg                                                                                                                >  2:1715581/1‑87 (MQ=255)                                                                          cTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAg                                                                                                                <  1:1715581/87‑1 (MQ=255)                                                                            cAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCt                                                          <  1:781871/139‑1 (MQ=255)                                                                              aTCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGa                                                                                                                    >  1:40333/1‑79 (MQ=255)                                                                              aTCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGa                                                                                                                    <  2:40333/79‑1 (MQ=255)                                                                               tCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTgcg                                                                >  1:1712406/1‑130 (MQ=255)                                                                               tCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTgcg                                                                <  2:1712406/130‑1 (MQ=255)                                                                                cTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGgcg                                                      <  1:775305/139‑1 (MQ=255)                                                                                    tGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGa                                                  >  2:1084627/1‑139 (MQ=255)                                                                                           ttACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCTCGTGCTgg                                                        >  2:1055418/1‑126 (MQ=255)                                                                                           ttACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTgg                                                        <  1:1055418/126‑1 (MQ=255)                                                                                              cGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                                 <  2:1452090/82‑1 (MQ=255)                                                                                              cGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                                 >  1:1452090/1‑82 (MQ=255)                                                                                                    cACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGAt                                  >  1:1721812/1‑139 (MQ=255)                                                                                                      ccTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTc                                <  2:679088/139‑1 (MQ=255)                                                                                                         tGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGa                             >  1:632214/1‑139 (MQ=255)                                                                                                                   ccTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACtt                                                                               <  2:589558/79‑1 (MQ=255)                                                                                                                   ccTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAAAGGCTATCCTGGTTAAACCTGTACtt                                                                               >  1:589558/1‑79 (MQ=255)                                                                                                                       ccGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTgcg                                                                >  2:509917/1‑90 (MQ=255)                                                                                                                       ccGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTgcg                                                                <  1:509917/90‑1 (MQ=255)                                                                                                                            tatCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACg                      <  1:1654587/127‑1 (MQ=255)                                                                                                                            tatCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACg                      >  2:1654587/1‑127 (MQ=255)                                                                                                                                 tcGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                                 <  2:864736/47‑1 (MQ=255)                                                                                                                                 tcGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                                 >  1:864736/1‑47 (MQ=255)                                                                                                                                 tcGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACtt                                                                               >  1:1326769/1‑65 (MQ=255)                                                                                                                                 tcGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACtt                                                                               <  2:1326769/65‑1 (MQ=255)                                                                                                                                 tcGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGcc                                     >  1:638170/1‑107 (MQ=255)                                                                                                                                 tcGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGcc                                     <  2:638170/107‑1 (MQ=255)                                                                                                                                 tcGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTACCGAAACggg     >  2:470842/1‑139 (MQ=255)                                                                                                                                 tcGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTACCGAAACggg     >  2:1203519/1‑139 (MQ=255)                                                                                                                                 tcGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGCTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTACCGAAACggg     >  1:936443/1‑139 (MQ=255)                                                                                                                                    ggTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTACCGAAACGGGTaa  >  2:1259161/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                       CGAAATGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAGCCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTACCGAAACGGGTAA  >  NC_000913/2022750‑2023019

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