Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,270,736	G→A	R4H (CGC→CAC)	pdeN →	putative c‑di‑GMP phosphodiesterase PdeN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,270,736	0	G	A	100.0%	51.1 / NA	19	R4H (CGC→CAC)	pdeN	putative c‑di‑GMP phosphodiesterase PdeN
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/7);  total (12/7)

GAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGCCTTCAGTTTTTAGTGGCCTGGCATAAGCACGAAGTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGAC  >  NC_000913/2270489‑2270976                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 gaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACacgc                                                                                                                                                                                                                                               >  1:388366/1‑251 (MQ=255)     cGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACACGCTccc                                                                                                                                                                                                                                           <  1:339427/251‑1 (MQ=255)                     ttACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACACGCTCCCAATTTTGGACGTAAGc                                                                                                                                                                                                                           >  2:9889/1‑251 (MQ=255)                                                                     gATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACACGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTg                                                                                                                                                                            >  2:278154/1‑250 (MQ=255)                                                                          aaTCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATAAGAGAAAATGATGTTTCCGCGTCTTGAACCAGGCTATAGTCCGGTCATTGTTATATTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACACGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCCGGCGTAATGCTTGCGATACTGGTGAGTTGc                                                                                                                                                                      >  1:240002/1‑251 (MQ=255)                                                                                                                    aGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACACGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGCCTTCAGTTTTTAGTGGCCTGGCATAAGCACGAAGTCAAATa                                                                                                                             <  1:184826/250‑1 (MQ=255)                                                                                                                                        catgcatcaGAGAAAATTATGTTTCCGCGTCTTGATCCAGTCTATAGTACGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACACGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGCCTTCAGTTTTTAGTGGCCTGGCATAAGCACGAAGTCAAATACGACACACTGATTACCGACGt                                                                                                        <  1:299765/247‑1 (MQ=255)                                                                                                                                                                             cAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACACGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGCCTTCAGTTTTTAGTGGCCTGGCATAAGCACGAAGTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGaaa                                                                   <  1:278154/251‑1 (MQ=255)                                                                                                                                                                                 cTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACACGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTACTGATTGCGATACTGGTGAGTTGCCTTCAGTTTTTAGTGGCCTGGCATAAGCACGAAGTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATcc                                                                >  1:414095/1‑250 (MQ=255)                                                                                                                                                                                                  ttATCTTTTGAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACACGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATCCTGGTGAGTTGCCTTCAGTTTTTAGTGGCCTGGCATAAGCACGACGTCAAATACGACACACTGATTACCGCCGTCCAAAAGAATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTcc                                                >  1:226510/1‑249 (MQ=255)                                                                                                                                                                                                      cTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACACGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGCCTTCAGTTTTTAGTGGCCTGGCATAAGCACGAAGTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGc                                          >  1:415461/1‑251 (MQ=255)                                                                                                                                                                                                                  gtCGTAATTTCAGGAAATTAACGGAATCATGTTCATACACGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGCCTTCAGTTTTTAGTGGCCTGGCATAAGCACGAAGTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGAt                               <  1:147875/250‑1 (MQ=255)                                                                                                                                                                                                                        atgtCAGGAAATTAACGGAATCATGTTCATACACTCACACAATTTTGGACGTAAGCTCCTGCTTCCCTGAATGGTTGCAGGCGTTATGATTGCGATATTGGTGAGTTGCCTTCATTTTTTAGTGGCCTGGCATAAGCACGAAGTCAAATCCGCCACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGTCGACTTTAAATCCACTACTGACCGGCTCCAGCCGCTGCCCTTAGATACCtccc                        <  2:401623/248‑4 (MQ=255)                                                                                                                                                                                                                             aGGAAATTAACGGAATCATGTTCATACACGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGCCTTCAGTTTTTAGTGGCCTGGCATAAGCACGAAGTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCa                    >  1:44665/1‑250 (MQ=255)                                                                                                                                                                                                                                    tAACGGAATCATGTTCATACACGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGCCTTCAGTTTTTAGTGGCCTGGCATAAGCACGAAGTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAAcc            >  1:352591/1‑251 (MQ=255)                                                                                                                                                                                                                                      aCGGAATCATGTTCATACACGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGCCTTCAGTTTTTAGTGGCCTGGCATAAGCACGAAGTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAAcccc          <  1:182146/251‑1 (MQ=255)                                                                                                                                                                                                                                       cGGAATCATGTTCATACACGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGCCTTCAGTTTTTAGTGGCCTGGCATAAGCACGAAGTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCg         >  2:122477/1‑251 (MQ=255)                                                                                                                                                                                                                                        ggAATCATGTTCATACACGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGCCTTCAGTTTTTAGTGGCCTGGCATAAGCACGAAGTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGCCCTTAGATACCTGCCAGCAAGCTAACCCCg         >  2:246837/1‑250 (MQ=255)                                                                                                                                                                                                                                              aTGTTCATACACGCTCCCCATTTTGCACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGCTTGAGATACTGGTGAGTTGCCTTCAGATTTTAGTGGCCTGGCATAAGCACGAAGTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGAc  >  2:181027/1‑251 (MQ=255)                                                                                                                                                                                                                                              aTGTTCATACACGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGCCTTCAGTTTTTAGTGGCCTGGCATAAGCACGAAGTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGCTACCTGCCAGCAAGCTAACCCCGCACTGAc  >  2:284452/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 GAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGACGTAAGCTCCTGCTTACCTGCATTGTTGCAGGCGTAATGATTGCGATACTGGTGAGTTGCCTTCAGTTTTTAGTGGCCTGGCATAAGCACGAAGTCAAATACGACACACTGATTACCGACGTACAAAAGTATCTCGATACCTATTTTGCCGACCTGAAATCCACTACTGACCGGCTCCAGCCGCTGACCTTAGATACCTGCCAGCAAGCTAACCCCGAACTGAC  >  NC_000913/2270489‑2270976

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