Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,564,818	A→G	S306P (TCA→CCA)	dosP ←	oxygen‑sensing c‑di‑GMP phosphodiesterase DosP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,564,818	0	A	G	100.0%	69.5 / NA	20	S306P (TCA→CCA)	dosP	oxygen‑sensing c‑di‑GMP phosphodiesterase DosP
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/11);  total (9/11)

CGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGA  >  NC_000913/1564683‑1564949                                                                                                                                        |                                                                                                                                    cgcgTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGcc                                                                                                                                 >  2:417421/1‑140 (MQ=255)  gcgTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCGTCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGccc                                                                                                                                >  1:21680/1‑140 (MQ=255)     ttCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAgt                                                                                                                             >  2:240039/1‑140 (MQ=255)                  cTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATccc                                                                                                                >  1:370245/1‑140 (MQ=255)                       ttCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATCCCGTTgc                                                                                                           <  1:417421/140‑1 (MQ=255)                         cTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATCCCGTTgcgc                                                                                                         <  2:245449/140‑1 (MQ=255)                                          ttAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAAc                                                                                        <  2:381041/140‑1 (MQ=255)                                           tAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCg                                                                                            <  1:408534/135‑1 (MQ=255)                                           tAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCg                                                                                            >  2:408534/1‑135 (MQ=255)                                                      aTCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGtt                                                                            <  2:155573/140‑1 (MQ=255)                                                                gcgcgcCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGAt                                                                  >  1:67834/1‑140 (MQ=255)                                                                     ccATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCg                                                                              <  1:178795/123‑1 (MQ=255)                                                                     ccATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCg                                                                              >  2:178795/1‑123 (MQ=255)                                                                               gACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCGATGTTACGACaa                                                   <  2:67834/140‑1 (MQ=255)                                                                                aCGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCGATGTTACGACaaa                                                  <  2:21680/140‑1 (MQ=255)                                                                                      ggTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGAttt                                            >  2:42631/1‑140 (MQ=255)                                                                                            tGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCa                                      <  1:144648/140‑1 (MQ=255)                                                                                                      tGCGCATTTTGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCTATGTTACGACAAATGATTTCCCCCATTTCATGaaa                            <  2:376134/140‑1 (MQ=255)                                                                                                             tttGAATTTCTGCACCGTGGGAAGATGGCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCggg                     <  2:332659/140‑1 (MQ=255)                                                                                                                                ggAAGATGGCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCga  >  2:226218/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                    CGCGTTCGATAAAGGCGCTGGTTTCTGCTCCTGACGAGGTTTTAATTTGCAGGATCCCCGCAGGCGCGCCATCACGCTGACGAATGGTCGCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGA  >  NC_000913/1564683‑1564949

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