Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,328,622	T→C	D496D (GAT→GAC)	rdcA →	regulator of diguanylate cyclase RdcA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,328,622	0	T	C	100.0%	87.7 / NA	33	D496D (GAT→GAC)	rdcA	regulator of diguanylate cyclase RdcA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/15);  total (18/15)

GTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTG  >  NC_000913/4328379‑4328866                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     gTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGGTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCAGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTgg                                                                                                                                                                                                                                               >  2:4001/1‑250 (MQ=255)       tGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGCCGCGCGAATACCTGGATTTTCGTGCGCACGGGCTACCCGTAGCCTGTGCACAATTGCGACAAAATATCCACCAGATCGACGAAAGTTTGCCGCTATTGCAACTCAATCAGGCCCAGGTGAGCGGCGAGATTAAACATGAAATCGCGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCCAGCCGCGCTGACGGTGc                                                                                                                                                                                                                                         >  1:241413/1‑251 (MQ=255)              cGATCTGACGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGAATTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGACGAAAGTTTGCAGCTATTGCAACTAAATCAGGCGCAGGTGAGCGGCGCGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCCACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAgc                                                                                                                                                                                                                                  >  2:168204/1‑251 (MQ=255)                           aTAAACTGTTGAATTACGCGCAGCAGGCGCGGGCATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGCCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTc                                                                                                                                                                                                                       >  2:198242/1‑249 (MQ=255)                                 tGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGgatgat                                                                                                                                                                                                               >  2:196591/1‑251 (MQ=255)                                        ttACGCGCAGAAGGCGCGCGCATACCTGGATTTTCGCGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGCCGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGc                                                                                                                                                                                                        >  2:172501/1‑251 (MQ=255)                                                       gcgGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCCCCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTc                                                                                                                                                                                         <  1:564/251‑1 (MQ=255)                                                        cgGGAATACCGTGATTTTCGTGCGCACGGGTGAAACGGCGCTTGTGAACAATTGCGGCAAAATAGCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGTCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGTGCTGAAGGGGCAGTTATGCGCCTTGTTTCAGGATGTTTCGGCGGCGTTAAGCGAGATTCg                                                                                                                                                                                        <  2:170656/251‑1 (MQ=255)                                                                           gTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGACGTTAAGCGGGATTCGTACCCGCTTTGAGACACTg                                                                                                                                                                     >  2:150670/1‑251 (MQ=255)                                                                                              cGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATGGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAAGCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCGGTGAGACACTGTTACAGACGGCGCAGAaca                                                                                                                                                  <  1:245523/251‑1 (MQ=255)                                                                                                    tgAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATGGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTTGCACTAACCTCCGCCAACCACTTTCTGCGTCAACAACAAGCCGCGCTGAAGGTGCAGTGAGCCGCCTTGTTTCAGGAAGATTCGGAGCCGTTAAGCGAGAGTCGAACCCGCTGTGAAACACTGTTACAGACGGCACAGAACACCATc                                                                                                                                             >  2:80290/1‑250 (MQ=255)                                                                                                                gCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTa                                                                                                                                <  1:51032/251‑1 (MQ=255)                                                                                                                    aaTATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCt                                                                                                                            >  1:223581/1‑251 (MQ=255)                                                                                                                     atatCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACGTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGAGGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTg                                                                                                                           <  1:167412/251‑1 (MQ=255)                                                                                                                             caGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAACCATGAAATCGAGCTGGCC‑TGACCTCCGCCAACCACTTTCTGCGTCAACCGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGTAGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCCGAACCCCCTCAGTCGGGACTTTACGCTGCGTTTTggc                                                                                                                  >  2:176161/1‑249 (MQ=255)                                                                                                                             caGATCCAAAAAAGCTTACAGCTATAGGAACTCAATCAGGCACCGGAGAACGGCGAGGTTACACATGAAATCGAGCTGGCCCTGACATCAGCCAACCACTTTCTGAGACAACAGCCAGACGCGCTGAAGGCGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGTGATTCGTCCCCGCTGAGAGACACTGTTACAGACGGCGCACAACACCGCCAGACGCGGCTTTACGCTGCGTTTTGc                                                                                                                   >  2:113007/1‑251 (MQ=255)                                                                                                                                        aaGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCTCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTGAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAAt                                                                                                        <  1:42337/251‑1 (MQ=255)                                                                                                                                          gttttccgccaTTGCAACTCAACCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAAACGGGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACCGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGAATTTACGCTGCGTTTTGCCGAGCTTGAATcc                                                                                                      <  2:68585/241‑1 (MQ=255)                                                                                                                                            gtgtaGCTATTGCAACTCAATCAGGCGCATGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGCCCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGAAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCAc                                                                                                    >  2:24035/5‑251 (MQ=255)                                                                                                                                                                        aGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCGTTGAGACGCTGTTACAGACTGCTCAGAACACCATCAGTCGCGACTTTACGCTGCGTGTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTAtt                                                                        <  1:265378/251‑1 (MQ=255)                                                                                                                                                                                ggcgAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTTCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCAt                                                                <  1:172501/251‑1 (MQ=255)                                                                                                                                                                                          aaCATGAAATCGAGCTTGCCCTCACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGTTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTGTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGcaacaa                                                      <  1:32267/251‑1 (MQ=255)                                                                                                                                                                                            cATGAAATCGAGCTGGCCCTGCCCTCCGCCAACCACTTGTGGCGTCAACAGCAAGGCGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGt                                                    <  2:95094/251‑1 (MQ=255)                                                                                                                                                                                                        cTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAAtt                                        >  1:2407/1‑251 (MQ=255)                                                                                                                                                                                                         tGGCCCTGACCTCCGCCAAACACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCCGTTCGCCGCCTTGTTTCCGGCTGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACACTTTGCAGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTg                                       >  2:246524/1‑251 (MQ=255)                                                                                                                                                                                                                     ccGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAggg                           >  1:17852/1‑251 (MQ=255)                                                                                                                                                                                                                      cGCCAACCCCTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCCCCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCATCAAGTCAACATGGAATCGAGAGAGTCAGGGt                          <  1:80290/251‑1 (MQ=255)                                                                                                                                                                                                                       gCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGtt                         >  2:155783/1‑251 (MQ=255)                                                                                                                                                                                                                         cAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTc                       >  1:69906/1‑251 (MQ=255)                                                                                                                                                                                                                            ccACTTTCTGCGTCAACAGCAAGACGCGCGGAAGGGGCAGTTATCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTGACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGGTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTc                    <  2:101690/251‑1 (MQ=255)                                                                                                                                                                                                                                      cGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTc           >  2:213054/1‑250 (MQ=255)                                                                                                                                                                                                                                         cAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCAtt        <  1:182226/250‑1 (MQ=255)                                                                                                                                                                                                                                              gCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTg  <  1:213054/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     GTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTG  >  NC_000913/4328379‑4328866

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