Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	913,656	A→G	R53R (CGT→CGC)	hcp ←	protein S‑nitrosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	913,656	0	A	G	100.0%	78.2 / NA	23	R53R (CGT→CGC)	hcp	protein S‑nitrosylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/9);  total (14/9)

TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCACGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACCACACATCCCCTGCGCGTATGAGCAGCCGTTTCCTGCCGGAGTACGGATA  >  NC_000913/913519‑913788                                                                                                                                          |                                                                                                                                     ttCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTcgcg                                                                                                                                    >  2:475491/1‑140 (MQ=255) ttCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTcgcg                                                                                                                                    >  2:696504/1‑140 (MQ=255)                       tCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGaa                                                                                                              >  2:615319/1‑140 (MQ=255)                       tCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGaa                                                                                                              <  2:646957/140‑1 (MQ=255)                                   cGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAgc                                                                                                  <  2:45817/140‑1 (MQ=255)                                    gTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAgcg                                                                                                 <  2:735661/140‑1 (MQ=255)                                            cAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATg                                                                                         >  2:178785/1‑140 (MQ=255)                                                                        ggTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTccg                                                             >  1:50968/1‑140 (MQ=255)                                                                         gTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTccgc                                                            <  1:611186/140‑1 (MQ=255)                                                                                 tGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTAc                                                    <  1:615319/140‑1 (MQ=255)                                                                                   agaAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACca                                                  <  1:475491/140‑1 (MQ=255)                                                                                   agaAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACca                                                  >  2:570644/1‑140 (MQ=255)                                                                                     aaaaaGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACcaca                                                <  1:696504/140‑1 (MQ=255)                                                                                        aaGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACCACACAt                                             >  1:83309/1‑140 (MQ=255)                                                                                             aCGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGTGTAAATCCAGAAGGTCCGAAGTTTCAGCCGTTTTACCACACATCCCCt                                        >  2:307684/1‑140 (MQ=255)                                                                                             aCGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACCACACATCCCCt                                        >  2:133103/1‑140 (MQ=255)                                                                                                                  aTCGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACCACACATCCCCTGCGCGTATGAGCAGCCGTTTc                   >  1:701690/1‑140 (MQ=255)                                                                                                                    cGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACCACACATCCCCTGCGCGTATGAGCAGCCGTTTCCt                 >  1:708330/1‑140 (MQ=255)                                                                                                                    cGTGGTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACCACACATCCCCTGCGCGTATGAGCAGCCGTTTCCt                 >  2:651214/1‑140 (MQ=255)                                                                                                                       ggTTGATGATGCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACCACACAACCCCTGCGCGTATGAGCAGCCGTTTCCTGcc              <  2:708330/140‑1 (MQ=255)                                                                                                                              gatgCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACCACACATCCCCTGCGCGTATGAGCAGCCGTTTCCTGCCGGAGTAc       >  2:194231/1‑140 (MQ=255)                                                                                                                              gatgCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACCACACATCCACTGGGCGTATGAACAGCCGTTTCCTGCCGGAGTAc       >  2:328407/1‑140 (MQ=255)                                                                                                                                tgCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACCACACATCCCCTGCGCGTATGAGCAGCCGTTTCCTGCCGGAGTACgg     <  1:178785/140‑1 (MQ=255)                                                                                                                                 gCCGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACCACACATCCCCTGCGCGTATGAGCAGCCGTTTCCTGCCGGAGTACGGa    >  2:163600/1‑140 (MQ=255)                                                                                                                                  ccGTATTCGCGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACAACACATCCACTGGGCGGATGAGCAGCCGTTACCTGCCGGAGTACGGAt   >  2:230927/1‑140 (MQ=255)                                                                                                                                   cGTATTCGCGCGCTTATACCGCCCACGCACAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACCACACATCCCCTGCGCGTATGAGCAGCCGTTTCCTGCCGGAGTACGGATa  <  1:230927/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                     TTCAGCGCCTCGCGCAGGGCAATCGCTTCACGAGCGTAGCCGACAATACGCGGAGAATCGAAGTTAACGTTGGTCAGGGTTGAGAAAAAGGCACGTGGCGCGAAGCTGTCAACATCGTGGTTGATGATGCCGTATTCACGCGCTTTTACCGCCCAGGCAGAAAGCCCTTGCAGCGCCGCGATGAGTAAATCCTGAAGGTCAGAAGTTTCCGCCGTTTTACCACACATCCCCTGCGCGTATGAGCAGCCGTTTCCTGCCGGAGTACGGATA  >  NC_000913/913519‑913788

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