Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,129,796	G→A	S12F (TCC→TTC)	flgN ←	flagellar biosynthesis protein FlgN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,129,796	0	G	A	96.2%	71.4 / NA	26	S12F (TCC→TTC)	flgN	flagellar biosynthesis protein FlgN
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (14/11);  minor base C (1/0);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGTTATGTTGATTCATTTGGCGTAGTTGCTGCGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGGACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCC  >  NC_000913/1129549‑1130039                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    gccggtattgtGATTCATTTGGCGTAGTTGCTGCGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATAGCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACa                                                                                                                                                                                                                                                  <  2:333052/241‑1 (MQ=255)                    tGGCGTAGTTGCTGCGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGc                                                                                                                                                                                                                               >  1:319162/1‑251 (MQ=255)                                        cacaGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCa                                                                                                                                                                                                            <  1:97784/250‑1 (MQ=255)                                             tAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTAct                                                                                                                                                                                                      >  1:62797/1‑251 (MQ=255)                                               aTTTCCTGCCAGCGTTGACTAATATCATCGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGATGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTActct                                                                                                                                                                                                    <  1:312014/251‑1 (MQ=255)                                                      gctagcgttgtcgaatttcgattcTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTTGTCTAGTTAATCCCGCGTTTGCAGCAGTGAGCTGTTTTGTTCTGTAATCTATTGCAACTGGCTGCTGTTTAGCTGCCTTTTATAGAGATGTTGCTGATCTTGGTCCATTTCCGGTTTGGGAGCGTTAAGCACACCGAACATCGTTTCTAGGATCTCTGCAAGACGTGTTATACGATCAGTTACTCTTCAAGTc                                                                                                                                                                                             <  1:156571/228‑1 (MQ=255)                                                                        cAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGCTGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGAt                                                                                                                                                                           >  1:49792/1‑251 (MQ=255)                                                                                              ttCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCg                                                                                                                                                     >  1:197193/1‑251 (MQ=255)                                                                                                                     ggTAATCCAGCGTCGCCATCAGTGAGCTTTTTTGTTCTGTAATCCATTGCTACTGGCTGCCGATGATCTGCCCCATAGAGAGATGTTTCTGCTTGTGATCCATTACCGTTTTGAGAGCGTTAATCACAGCGAACATCGGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGTTGCGTTGATCAGCGCATCGGCAATTTTGACGGTGTCCATTTTTAGTTAACCggt                                                                                                                              <  2:300461/251‑3 (MQ=255)                                                                                                                             aGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCg                                                                                                                      <  2:319162/251‑1 (MQ=255)                                                                                                                                gTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGGTGATCTGCACCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTAAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGGCTTGCTGCGCTTCGGTGATCAACGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAGTCGccc                                                                                                                   >  1:228839/1‑250 (MQ=255)                                                                                                                                         agTGGGCTTTTATGTTCTGTAATCCATTGCAACTGGCTGCCGATGATCTGCCCCATAGAGAGATGTAGCTTCTCTTGATCCATTACCGTTTTTAGATCGTTAAGCACAGCGCTCATCTGGTCGAGGATCACTGCAAGACGTGTCATACGATCAGTTACTCAGCAAGTCTTGCTGCGATTCGATGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACAGTTCCGAAACGCAAGTTTAAACg                                                                                                          >  2:233961/1‑251 (MQ=255)                                                                                                                                                   tttGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTGCACCGTTACGAATCGCCAGTTTTAACGCTTCGacac                                                                                                 <  2:307826/250‑1 (MQ=255)                                                                                                                                                                 cATTGCAACTGGCTGCCGTTGATTTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGatat                                                                                  <  1:370937/251‑1 (MQ=255)                                                                                                                                                                     gCAACTGGCTGCCGTGGAGCTGCCCCATAGAGAGATGTTGCGGCGCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGTTGTCCATTTTTAGTTCACCGTTTCGAAACGCCAGTTTTAACGCTTCTACACGTTCAAGATTGATATCAct                                                                              <  1:234743/251‑1 (MQ=255)                                                                                                                                                                                  gTTGATCTGCCCCATAGAGAGATGTTTCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTAAAGATTGCTATCACTGCTGCCGGGTTGc                                                                 >  1:294135/1‑251 (MQ=255)                                                                                                                                                                                    tGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCAt                                                               >  1:284571/1‑251 (MQ=255)                                                                                                                                                                                                           tgctgcTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTATGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCCTCGGCAATTTTGCCGGTGTACATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCtt                                         >  1:255189/1‑250 (MQ=255)                                                                                                                                                                                                                     gATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACAct                              >  1:305470/1‑251 (MQ=255)                                                                                                                                                                                                                     gATCCATTACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACAct                              <  1:115473/251‑1 (MQ=255)                                                                                                                                                                                                                           ttACCGTTTTGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTAGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACActggtgct                        >  1:260150/1‑251 (MQ=255)                                                                                                                                                                                                                                   ttGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAg                 <  1:135951/250‑1 (MQ=255)                                                                                                                                                                                                                                    tGAGATCGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGc               >  1:148450/1‑251 (MQ=255)                                                                                                                                                                                                                                       gaTCGTTAAGCACGGAGAGCATCTGGTCGAGGATCTCTGCAAGACGGGACATACGATCAGTTACTCTGTAAGTCTTACTGCGCTTCGGTGATCAGCGCATCGGCGATTTTGCCGATGACCATTTTTAGTTCACCGTTACGAAGCGCCAGTTTTAACGCTTCGACACGTTCAAGCTTGATAGCCCTGATGCCGAGTTTCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACActgatgcggttggtggcgtg            >  2:139858/1‑238 (MQ=255)                                                                                                                                                                                                                                          cGTTAAGCACAGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGgttgt         >  1:315715/1‑251 (MQ=255)                                                                                                                                                                                                                                                 cacaGCGAACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGttttttcc  >  1:105684/1‑248 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    GCCGTTATGTTGATTCATTTGGCGTAGTTGCTGCGTTTTCACAGTAATTTCCTGCCAGCGTTGACTAATATCAACGCTATTTGCTGTATTGGGTTCTTTCCTGCGTAACTGTTCGAGGTAATCCAGCGTCGCCAGCAGTGAGCTTTTTTGTTCTGTAATCCATTGCAACTGGCTGCCGTTGATCTGCCCCATAGAGAGATGTTGCTGCTCTTGATCCATTACCGTTTTGAGATCGTTAAGCACAGCGGACATCTGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCC  >  NC_000913/1129549‑1130039

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