Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,565,918	G→A	A408V (GCC→GTC)	dosC ←	diguanylate cyclase DosC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,565,918	0	G	A	95.7%	57.6 / ‑6.2	23	A408V (GCC→GTC)	dosC	diguanylate cyclase DosC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (14/8);  minor base C (0/1);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

TCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGGCTTTCAGTTTGGTTTTCTCCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAAT  >  NC_000913/1565678‑1566150                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         tCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCCCCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGACTTTCAGttt                                                                                                                                                                                                                                >  2:79301/1‑251 (MQ=255)      gtcGGGGAAAAATATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCAGCCTGGCGCATATAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGTATGTGAGAGGGCAATATCTTCGCCGTTTGCGCCTTTCAGTTTGGttt                                                                                                                                                                                                                           <  2:58435/249‑1 (MQ=255)             aaaaaaaTGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCAc                                                                                                                                                                                                                    >  1:336801/1‑251 (MQ=255)                          cggcggCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGATTCTACCTCGTCTTTTGGCGCTATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCACCCGACTGCGAAt                                                                                                                                                                                                        >  1:276458/1‑250 (MQ=255)                                         cATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCACCCGACTGCGAATACGTTCTGCGGTACGt                                                                                                                                                                                        >  1:269131/1‑251 (MQ=255)                                            cGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCACCCGACTGCGAATACGTTCTGCGGTACGTaa                                                                                                                                                                                      >  1:179655/1‑250 (MQ=255)                                                 aGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCACCCGACTGCGAATACGTTCTGCGGTACGTAAcgtttc                                                                                                                                                                                <  1:300844/251‑1 (MQ=255)                                                 aGCTTCATGATTACCTCTGCAGCCTGGCGCATTTAAAGACTGGTTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTAATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCTATGTAAATTGAGAGGGCAATATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCACCCGACTGCGAATACGTTCTGCGGTACGTAAcgttt                                                                                                                                                                                 <  2:382944/250‑1 (MQ=255)                                                                          ggCGCATCTAAAGACTGGCTTTCCAGAGTTCAACGCGGTTTCTACCTCGTCTTTTGGCGATATACAGATCTTCATCGGCTATTTGAATGAGGCTCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGTAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCTCCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACa                                                                                                                                                       <  2:13880/251‑1 (MQ=255)                                                                                               tCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGccccc                                                                                                                                   <  1:79301/250‑1 (MQ=255)                                                                                                       tCAACACGGCTTCTAACTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACAATTAAACATGGCGGCACCGTTGGAAAGTGAGAGGGCAATATCTTCGCCATTTGCGACTTTAAGTTTGGTTTTCTCAACCCGACTGCGACTACGTTCTGCGGTACGTAACGTTTAGCTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGa                                                                                                                           >  2:56420/1‑250 (MQ=255)                                                                                                                                     atatACAGAGCTTCATCGGCTATTTGAATGATGCGCTCATAGTCAGGATGATCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTTAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTAATGCGGACGTTGttat                                                                                             >  2:391783/1‑247 (MQ=255)                                                                                                                                     atatACAGAGCTTACTCGGCAATTTGAATGACGCGATCATAGTCAGGATGCTCATTAAACATGGAGGCACCGGTGGAAAGGCAAAGGGCAATATCCTAGCCCTCTGCGACTTTCTGTTTGGTTTTCTCCACTCGACTACGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCACAAGCTTCAGTCACAACAATGATAAATTCATCGCCCCCGTACCGGAAAACATAATCACTACTGCGGACGTTGTNAt                                                                                             >  2:355335/1‑250 (MQ=255)                                                                                                                                           agagCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGc                                                                                      >  2:344592/1‑251 (MQ=255)                                                                                                                                                              gAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGa                                                                   <  1:216924/251‑1 (MQ=255)                                                                                                                                                                          cATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACTTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCa                                                       <  2:68406/251‑1 (MQ=255)                                                                                                                                                                                          aTTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACg                                       >  1:9274/1‑251 (MQ=255)                                                                                                                                                                                                   ggcgtgacggATGGAGAGTGAGAGGGCGATGTCTTCGCCCTGTGCGACTTTCAGGTTGGTTGTCTCCTGCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTAATGCGGACGTTGGCATAAAAGGCCTGGGAGTCTGTACGCAGAAGTTCATCACCAGGGTTAGGGTCCCACGTATCGTTGa                              <  2:276458/242‑1 (MQ=255)                                                                                                                                                                                                      ggcACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGAtct                           >  2:282009/1‑251 (MQ=255)                                                                                                                                                                                                                  aaGTGAGAGGGCAATATCTTGGAGGTTTGCCACTTTCAGTTTGGTTTTCTAAACGAGACTGAGAATAAGTGAGGAGGTACGTAACGTTTAGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGAAAAAGTAGCGGAAAAAATAATAACTAATGAGGAAGTTGTCATAAAAGGCATGAGAGACGTGAAGCAGAAGGTGAGAACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTAt               <  1:124475/251‑1 (MQ=255)                                                                                                                                                                                                                         agGGCAATATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTc        >  1:197085/1‑251 (MQ=255)                                                                                                                                                                                                                              aaTATCTTCGCCGTTTGCGACTTTCAGTTTGGTTTTCTCCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCaata    >  1:230969/1‑250 (MQ=255)                                                                                                                                                                                                                               atatCTTCGCAGCTTGCGACTTTCCGATTGGTTTTCTCCACCCGACTGCCAGTACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTCCAGTCAAAACAATGATAAATTCATCGCCACCGTAGCGGACAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCAACGTATCGTTGATCTCTTTGAATTTATCAACGTCaataat  >  2:167837/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         TCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGAGAGGGCAATATCTTCGCCGTTTGCGGCTTTCAGTTTGGTTTTCTCCACCCGACTGCGAATACGTTCTGCGGTACGTAACGTTTCGTTTTCAGAAGCTTCAGTCAAAACAATGATAAATTCATCGCCCCCGTAGCGGAAAACATAATCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAAT  >  NC_000913/1565678‑1566150

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