Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,969,993	T→G	15.9%	T84P (ACT→CCT)	fliC ←	flagellin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,969,993	0	T	G	15.9%	104.8 / 14.2	44	T84P (ACT→CCT)	fliC	flagellin
Reads supporting (aligned to +/- strand):  ref base T (20/17);  new base G (3/4);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CCGAAAGTGGTGGTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTT  >  NC_002947/4969845‑4970134                                                                                                                                                     |                                                                                                                                              ccGAAAGTGGTGGTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAACCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGt                                                                                                                                               >  3:209918/1‑149 (MQ=255)     aaGTGGTGGTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGcc                                                                                                                                           <  2:245644/149‑1 (MQ=255)          gttgtgtTAGCTATACGTGTGAGCTGGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTAGAGCCGTTCGCCGATTGCAGGGCCAGTTCTCGCATACGCTGCATGATGTTGCCTGAAGGCGCGATgg                                                                                                                                       <  7:5009/145‑1 (MQ=255)          gtggtgTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGGCGCCATGgc                                                                                                                                      >  7:133459/1‑149 (MQ=255)             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                                                                                                                              >  5:185583/1‑134 (MQ=255)                                 cTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGGCGCCATGGCGACTTcag                                                                                                                              <  6:185583/134‑1 (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:35903/1‑149 (MQ=255)                                                       tCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCtt                                                                                         >  7:151798/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                                                        >  2:155618/1‑149 (MQ=255)                                                                                        ttCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGt                                                        <  1:147059/149‑1 (MQ=255)                                                                                            aGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGa                                                    >  8:199841/1‑149 (MQ=255)                                                                                            aGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGa                                                    >  3:109434/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)                                                                                                                                                     |                                                                                                                                              CCGAAAGTGGTGGTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCAGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTT  >  NC_002947/4969845‑4970134

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: ‑