Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,978,083	+A	coding (115/351 nt)	flhD ←	DNA‑binding transcriptional dual regulator FlhD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,978,083	1	.	A	94.4%	56.1 / ‑2.4	18	coding (115/351 nt)	flhD	DNA‑binding transcriptional dual regulator FlhD
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (6/11);  total (6/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TT‑TATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAAT  >  NC_000913/1977845‑1978314                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        gATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTGCAGTA‑TATGCCGAGAc                                                                                                                                                                                                                              <  2:358473/251‑1 (MQ=255)                    gcAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCttt                                                                                                                                                                                                           <  2:150504/251‑1 (MQ=255)                                          aaCATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCCAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCtgt                                                                                                                                                                                     <  2:187154/251‑1 (MQ=255)                                                           gtgtTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGa                                                                                                                                                                    <  1:130803/251‑1 (MQ=255)                                                                                    aaTTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCCACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCCTTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTc                                                                                                                                           >  2:154075/1‑251 (MQ=255)                                                                                                          cgGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCatt                                                                                                                      >  2:153396/1‑250 (MQ=255)                                                                                                           gggAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCattat                                                                                                                    <  2:209496/251‑1 (MQ=255)                                                                                                            ggAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCattatt                                                                                                                   <  2:6151/251‑1 (MQ=255)                                                                                                                   tGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCGGCAACTCGGAGGTATGCATTATTcccac                                                                                                              <  2:273973/249‑1 (MQ=255)                                                                                                                      gTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTcccaccc                                                                                                            >  2:310722/1‑248 (MQ=255)                                                                                                                               gTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAAc                                                                                                <  1:311563/251‑1 (MQ=255)                                                                                                                                                                cAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATg                                                                                                                                                          <  1:53609/160‑1 (MQ=255)                                                                                                                                                                cAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATg                                                                                                                                                          >  2:53609/1‑160 (MQ=255)                                                                                                                                                                                                     gaggAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAGGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCAtt                          <  2:110643/249‑1 (MQ=255)                                                                                                                                                                                                            agTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAAt                     >  1:180546/1‑249 (MQ=255)                                                                                                                                                                                                                              tGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAAGAGTTGCGATAAGCTGCAAt  <  1:100941/250‑1 (MQ=255)                                                                                                                                                                                                                                        tcttcA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGACCAATCCAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAg                                                                                                          >  1:138953/1‑136 (MQ=255)                                                                                                                                                                                                                                        tcttcA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAg                                                                                                          <  2:138953/136‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        GATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TT‑TATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAAT  >  NC_000913/1977845‑1978314

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