Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,970,004	C→A	12.2%	G80V (GGC→GTC)	fliC ←	flagellin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,970,004	0	C	A	12.2%	112.0 / 6.7	41	G80V (GGC→GTC)	fliC	flagellin
Reads supporting (aligned to +/- strand):  ref base C (20/16);  new base A (2/3);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

GTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAGGCCA  >  NC_002947/4969857‑4970141                                                                                                                                                    |                                                                                                                                          gtgtTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGcc                                                                                                                                          >  5:88166/1‑149 (MQ=255)      aGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGGCGCCATGGCGACTTca                                                                                                                                      <  8:133459/148‑1 (MQ=255)       gCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGGCGCCATGGCGACTTcag                                                                                                                                     <  1:206950/148‑1 (MQ=255)        cAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTcagca                                                                                                                                   <  8:29666/149‑1 (MQ=255)             cGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGttt                                                                                                                               <  6:88166/148‑1 (MQ=255)                     cTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGGCGCCATGGCGACTTcag                                                                                                                                     <  6:185583/134‑1 (MQ=255)                     cTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGGCGCCATGGCGACTTcag                                                                                                                                     >  5:185583/1‑134 (MQ=255)                        gCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATa                                                                                                                       >  2:128791/1‑145 (MQ=255)                        gCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAACCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATaga                                                                                                                     <  4:209918/147‑1 (MQ=255)                                          cTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCt                                                                                                 >  7:206237/1‑149 (MQ=255)                                           tCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCtt                                                                                                >  7:151798/1‑149 (MQ=255)                                           tCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCtt                                                                                                >  7:35903/1‑149 (MQ=255)                                             ctgctgCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCtt                                                                                                <  5:136640/147‑1 (MQ=255)                                                   cATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGc                                                                                          <  8:151798/147‑1 (MQ=255)                                                     tGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGTGAAAGTTCACGCATACGCTGCATTATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTTCGCAATAGAGATACAATCGTTGGCGTTATTTACCGCAAc                                                                                       >  1:37657/1‑148 (MQ=255)                                                       gAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCt                                                                                    <  8:206237/149‑1 (MQ=255)                                                                tcgtcATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCGTCAGCAGTTTGCGCAATAGAGATACCATCGTTTGCGTTCTGGACCGCCACGCTCCGACCtat                                                                           >  4:28776/1‑147 (MQ=255)                                                                   tcATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCtttga                                                                         <  7:235379/148‑1 (MQ=255)                                                                   tcATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCtttga                                                                         >  8:235379/1‑148 (MQ=255)                                                                    cATCGCTGTTCGAGCCGTTCGACTATTGCAGTGCCAGTTCACGCATACGCTTCATGATGTTGCCTGAAGTCTCCATGGCTCCTTAAGCAGTTTGCTCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGAGCtttgat                                                                        <  2:37657/148‑1 (MQ=255)                                                                        gCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCtttgatttga                                                                    <  6:179363/148‑1 (MQ=255)                                                                        gCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAg                                                                   >  1:119246/1‑149 (MQ=255)                                                                          ggTTCGAGCCGTTCGCCGATTGCAGGGACAGTTCAAGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTg                                                                 <  3:28776/148‑1 (MQ=255)                                                                            ttCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTgg                                                                >  1:79603/1‑148 (MQ=255)                                                                            ttCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGt                                                               <  1:147059/149‑1 (MQ=255)                                                                            ttCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGt                                                               >  2:155618/1‑149 (MQ=255)                                                                                aGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGa                                                           >  3:109434/1‑149 (MQ=255)                                                                                aGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGa                                                           >  8:199841/1‑149 (MQ=255)                                                                                  ccGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACg                                                         >  7:200629/1‑149 (MQ=255)                                                                                    gTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGt                                                       <  1:128791/149‑1 (MQ=255)                                                                                     ttCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGtt                                                      <  2:79603/149‑1 (MQ=255)                                                                                           gATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACAATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTTATTTGAGTGGTCACAGGGTTGgaga                                                 >  8:140145/1‑148 (MQ=255)                                                                                                   ggCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGc                                        <  8:35903/149‑1 (MQ=255)                                                                                                         ttCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTgcgg                                  >  7:66036/1‑149 (MQ=255)                                                                                                            aCGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGcgt                               >  7:149036/1‑149 (MQ=255)                                                                                                             cGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGcgtc                              <  7:199841/149‑1 (MQ=255)                                                                                                             cGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGcgtc                              <  1:155618/149‑1 (MQ=255)                                                                                                                  aCGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCAGCAGGCCTGCGGCGTCGTCtt                         >  7:225394/1‑149 (MQ=255)                                                                                                                     cTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCttttg                      >  5:72724/1‑147 (MQ=255)                                                                                                                     cTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTgg                      >  5:76480/1‑149 (MQ=255)                                                                                                                                   cTGAAGGCGCCATGGCGACTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGAtttt         >  1:184009/1‑148 (MQ=255)                                                                                                                                           gCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAGGCCa  >  1:246514/1‑147 (MQ=255)                                                                                                                                                    |                                                                                                                                          GTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAGGCCA  >  NC_002947/4969857‑4970141

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