Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,407,462	T→C	T112A (ACG→GCG)	ydaM ←	diguanylate cyclase, csgD regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,407,462	0	T	C	100.0%	174.1 / NA	47	T112A (ACG→GCG)	ydaM	diguanylate cyclase, csgD regulator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/27);  total (20/27)

ACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGT  >  NC_000913/1407335‑1407592                                                                                                                                |                                                                                                                                   acgaTTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCtt                                                                                                                         >  2:343784/1‑139 (MQ=255)     tttGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGAt                                                                                                                     >  1:875741/1‑139 (MQ=255)             tctcGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGa                                                                                                             <  2:880304/139‑1 (MQ=255)                         tcaatcaaCAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCg                                                                                                 >  1:159788/1‑139 (MQ=255)                          caatcgaCAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCGTAAAAGCCTTGTTTTTTGCGCTGATAATTTGCCGGGCGACTGTCTTTTAATCCCAGCGCTGCCGGCGCTTCGATACCTTGGAAAATAATCACCgg                                                                                                <  2:481594/139‑1 (MQ=255)                             tcaaCAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGttc                                                                                             >  1:608754/1‑139 (MQ=255)                                       tGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTc                                                                                   >  2:167546/1‑139 (MQ=255)                                         gTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAg                                                                                 >  1:40533/1‑139 (MQ=255)                                          tGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGc                                                                                <  2:875741/139‑1 (MQ=255)                                                    gTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCaaaa                                                                      >  2:971390/1‑139 (MQ=255)                                                           aaaaGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCg                                                               >  2:590349/1‑139 (MQ=255)                                                            aaaGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGa                                                              <  1:488616/139‑1 (MQ=255)                                                                 gcgcATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCt                                                         <  2:59836/139‑1 (MQ=255)                                                                 gcgcATAAAAACCTTGTTTTTAGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGTAAAACAAGCCGACAGCt                                                         <  2:700954/139‑1 (MQ=255)                                                                   gcATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCa                                                       >  2:389529/1‑139 (MQ=255)                                                                   gcAAAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAATACAAGCCGACAGCTCa                                                       >  2:761088/1‑139 (MQ=255)                                                                        aaaaCCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCt                                                  <  2:159788/139‑1 (MQ=255)                                                                        aaaaCCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCt                                                  <  1:971390/139‑1 (MQ=255)                                                                           aCCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGtt                                               <  1:607118/139‑1 (MQ=255)                                                                               tGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCct                                           <  1:617525/139‑1 (MQ=255)                                                                               tGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCct                                           >  2:897423/1‑139 (MQ=255)                                                                                    tttGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCtcttct                                      <  2:559080/139‑1 (MQ=255)                                                                                    tttGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCtcttct                                      <  1:905000/139‑1 (MQ=255)                                                                                           tGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCgg                                                                                                                              <  2:875869/44‑1 (MQ=255)                                                                                           tGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCgg                                                                                                                              >  1:875869/1‑44 (MQ=255)                                                                                              tAATTTGCCGAGCGACGGGCTCTTAATCCCAGCGCTGCTGGCGCTTCGATCGCTTCGTACATAATCACCGGTTCTGGTTCTCTCAGCTTTCGGAAAACAAGCCGACAGCTCAATGCTGTTGGCTGTTCTTAACGCTGaa                            <  1:279841/139‑1 (MQ=255)                                                                                                   tGCCGAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTg                       <  2:847969/139‑1 (MQ=255)                                                                                                       gAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGa                                                                                                                      >  2:524202/1‑40 (MQ=255)                                                                                                       gAGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGa                                                                                                                      <  1:524202/40‑1 (MQ=255)                                                                                                        aGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCtcttct                                      >  1:90297/1‑119 (MQ=255)                                                                                                        aGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCtcttct                                      <  2:90297/119‑1 (MQ=255)                                                                                                        aGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGAt                  <  1:291542/139‑1 (MQ=255)                                                                                                        aGCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGAt                  <  2:40533/139‑1 (MQ=255)                                                                                                         gCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCtt                                                                                                                >  2:75091/1‑44 (MQ=255)                                                                                                         gCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCtt                                                                                                                <  1:75091/44‑1 (MQ=255)                                                                                                         gCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAaata                                                                                                        <  1:902438/52‑1 (MQ=255)                                                                                                         gCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAaata                                                                                                        >  2:902438/1‑52 (MQ=255)                                                                                                         gCGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGAtt                 >  1:575797/1‑139 (MQ=255)                                                                                                          cGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACGTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGAttt                <  1:705633/139‑1 (MQ=255)                                                                                                          cGACTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGAttt                <  1:167546/139‑1 (MQ=255)                                                                                                            aCTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGAt                                                                                                                     >  2:181945/1‑36 (MQ=255)                                                                                                            aCTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGAt                                                                                                                     <  1:181945/36‑1 (MQ=255)                                                                                                            aCTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATACTCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCt                                                  <  2:273961/103‑1 (MQ=255)                                                                                                            aCTGGCTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCt                                                  >  1:273961/1‑103 (MQ=255)                                                                                                                 cTTTTAAACCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGataat         <  1:343784/139‑1 (MQ=255)                                                                                                                        aCCCAGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCgtgt  <  1:885501/139‑1 (MQ=255)                                                                                                                            aGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACa                                                            <  2:1061012/77‑1 (MQ=255)                                                                                                                            aGCGCTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACa                                                            >  1:1061012/1‑77 (MQ=255)                                                                                                                                |                                                                                                                                   ACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGACTGGCTTTTAAACCCAGCGTTGCCGGCGCTTCGATACCTTCGAAAATAATCACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGT  >  NC_000913/1407335‑1407592

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