Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,632	A→G	I225I (ATT→ATC)	arcA ←	DNA‑binding transcriptional dual regulator ArcA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,639,632	0	A	G	93.8%	39.4 / ‑4.0	17	I225I (ATT→ATC)	arcA	DNA‑binding transcriptional dual regulator ArcA
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (8/7);  total (8/9)
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.95e-01

TTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTC  >  NC_000913/4639387‑4639874                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   gttatgGGAACCAGTGTGATGGTGGTGGGGCTGTCCGGCATGATGTTTGTGACTCGTAAATTCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGAtggt                                                                                                                                                                                                                                               <  2:89709/247‑1 (MQ=255)       ggAACCAGTGTGCTGGAGGTGGCGCCGGCCGGCATGATGTTTGGGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTAAAATGGACGATCAATTACGAATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAACAAGCGCCGTTTTTTTTGTCGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAAACTTCACCGTGGATGGGGGCgat                                                                                                                                                                                                                                         >  1:308714/1‑251 (MQ=255)            cAGTGTGCGTGTGGTGGCGCTGGCCGGCATGATGTTTTTGACTCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGAttt                                                                                                                                                                                                                                    <  1:356645/251‑1 (MQ=255)                                                       cGTAATATCGACTGGGATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAATGGACGAGCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGACGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcggcg                                                                                                                                                                                         <  1:172394/251‑1 (MQ=255)                                                                                                         agaagTTACAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAATCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACATAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGTCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCTCGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCa                                                                                                                                       <  1:369281/251‑1 (MQ=255)                                                                                                                aCAATGGACGATCAATTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTttctt                                                                                                                                >  1:101707/1‑251 (MQ=255)                                                                                                                      gACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTcagcag                                                                                                                          <  1:101945/251‑1 (MQ=255)                                                                                                                                             aTAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCCCGGCCGGTCATTTTCTTCAGCAGTTCAGCACGCGACTGAATTTTGc                                                                                                   >  1:68289/1‑251 (MQ=255)                                                                                                                                               aaGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCt                                                                                                 >  1:301040/1‑251 (MQ=255)                                                                                                                                                 ggTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTgg                                                                                               <  1:146153/251‑1 (MQ=255)                                                                                                                                                                                  gTTTTTTTTGACGGGGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGTTCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTTCCTGGGTTTTCACAGAAGGGAAGAATGGAACGGAACTc                                                              >  1:333898/1‑251 (MQ=255)                                                                                                                                                                                               gtggtAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGGGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGAATGAATTTTGCCTGGGTTTTCACAGTAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGc                                                 >  1:357967/1‑251 (MQ=255)                                                                                                                                                                                                 ggtAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCCCGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGGCTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCtt                                               >  1:224012/1‑251 (MQ=255)                                                                                                                                                                                                    aaaGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTa                                            <  1:98443/251‑1 (MQ=255)                                                                                                                                                                                                                 cTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTAGTGTGGTTTCAGCTCCCGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGGCTGAATTTTGCCTGGGTTTTCACCGACGTGAAGCAAGGCGCGGAACTCGCTGCGCGGCAGCATGTACTGCTCGCCATc                                >  1:57159/1‑250 (MQ=255)                                                                                                                                                                                                                                 gAAGCGATAACCTGTGCCGGGGATGGGGGCTATGATTTCCGGCGTATCCGGCTTAGATTCGAAAAGTTTACGAATCCGGCGTATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGCCTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGTCGCGGAACTCGCGGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCaacgaac               <  2:308999/251‑1 (MQ=255)                                                                                                                                                                                                                                               cACCGTGTATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTc  <  1:339697/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   TTGTTGGGAACCAGTGTGCTGGTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGGTCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCTGGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCATCAGGGCCGATCAACGAACGGCTGTTGATGTC  >  NC_000913/4639387‑4639874

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