Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	477,636	T→C	T328A (ACG→GCG)	pdeB ←	c‑di‑GMP phosphodiesterase PdeB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	477,636	0	T	C	100.0%	107.7 / NA	30	T328A (ACG→GCG)	pdeB	c‑di‑GMP phosphodiesterase PdeB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/14);  total (16/14)

GGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATGCTGCTGTGGATGCTGACGCAGCCAGTCGCCCATATCTTCAAAGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGTTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCA  >  NC_000913/477500‑477761                                                                                                                                         |                                                                                                                              gggATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATGCTGCTGTGGATGCTGACGCAGCCAGTCGCCCATATCTTCAAAGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTg                                                                                                                            <  2:143446/140‑1 (MQ=255)    aTTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATGCTGCTGTGGATGCTGACGCAGCCAGTCGCCCATATCTTCAAAGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCtg                                                                                                                         >  1:256309/1‑140 (MQ=255)               ggTGAGCACGGGGGATTCAAGATTGATCGAAATATGCTGCTGTGGATGCTGACGCAGCCAGTCGCCCATATCTTCAAAGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGaa                                                                                                              <  2:256309/140‑1 (MQ=255)                 tGAGCACGGGGGATTCAAGATTGATCGAAATATGCTGCTGTGGATGCTGACGCAGCCAGTCGCCCATATCTTCAAAGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATa                                                                                                            <  1:231659/140‑1 (MQ=255)                         ggggATTCAAGATTGATCGAAATATGCTGCTGTGGATGCTGACGCAGCCAGTCGCCCATATCTTCAAAGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATc                                                                                                    <  2:72972/140‑1 (MQ=255)                               tCAAGATTGATCGAAATATGCTGCTGTGGATGCTGACGCAGCCAGTCGCCCATATCTTCAAAGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGa                                                                                              <  2:221012/140‑1 (MQ=255)                                        aTCGAAATATGCTGCTGTGGATGCTGACGCAGCCAGTCGCCCATATCTTCAAAGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACt                                                                                     >  2:217329/1‑140 (MQ=255)                                                      ctgTGGATGCTGACGCAGCCAGTCGCCCATATCTTCAAAGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCg                                                                       >  1:90807/1‑140 (MQ=255)                                                        gtgGATGCTGACGCAGCCAGTCGCCCATATCTTCAAAGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGc                                                                     <  2:90807/140‑1 (MQ=255)                                                                tGACGCAGCCAGTCGCCCATATCTTCAAAGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAgcgcgc                                                             <  1:217329/140‑1 (MQ=255)                                                                 gACGCAGCCAGTCGCCCATATCTTCAAAGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTTTCTGCGGCCAGCGCGcc                                                            <  2:229983/140‑1 (MQ=255)                                                                              gCCCATATCTTCAAAGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATATGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCAc                                               >  1:325587/1‑140 (MQ=255)                                                                                        tCAAAGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAAtttt                                     >  1:377705/1‑140 (MQ=255)                                                                                           aaGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGcc                                  <  1:72805/140‑1 (MQ=255)                                                                                               cGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACaa                                         >  2:24648/1‑129 (MQ=255)                                                                                               cGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACaa                                         <  1:24648/129‑1 (MQ=255)                                                                                               cGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTg                              >  1:260962/1‑140 (MQ=255)                                                                                               cGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTg                              >  2:350969/1‑140 (MQ=255)                                                                                                  ttcttATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCt                           >  2:56639/1‑140 (MQ=255)                                                                                                            cAGTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACaa                 <  2:102669/140‑1 (MQ=255)                                                                                                              gTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATc               >  2:147262/1‑140 (MQ=255)                                                                                                              gTAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATc               <  1:350969/140‑1 (MQ=255)                                                                                                               tAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATGAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGga                      <  2:354869/132‑1 (MQ=255)                                                                                                               tAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATCAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTCGCTAAGga                      >  1:354869/1‑132 (MQ=255)                                                                                                               tAGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCg              >  1:136820/1‑140 (MQ=255)                                                                                                                aGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACaa                 <  1:330615/136‑1 (MQ=255)                                                                                                                aGCGTCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACaa                 >  2:330615/1‑136 (MQ=255)                                                                                                                    tCAATGGCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGa         >  2:75044/1‑140 (MQ=255)                                                                                                                         ggCTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTg    >  2:80465/1‑140 (MQ=255)                                                                                                                           cTCAGAAAGGCCCGCTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCa  >  2:17906/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                              GGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATGCTGCTGTGGATGCTGACGCAGCCAGTCGCCCATATCTTCAAAGACGCTTCTTATAATCAGTAGCGTCAATGGCTCAGAAAGGCCCGTTTGCTGTGCCAGCGGAATAAAACTATCTGGTGACAACCAACTACCGTCTGTCTGCGGCCAGCGCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCA  >  NC_000913/477500‑477761

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