Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,407,343	T→A	G151G (GGA→GGT)	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,343	0	T	A	100.0%	230.2 / NA	67	G151G (GGA→GGT)	ydaM	diguanylate cyclase, csgD regulator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (27/40);  total (27/40)

CCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCG  >  NC_000913/1407247‑1407437                                                                                                 |                                                                                               cccGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCa                                                                                             >  2:1438671/1‑100 (MQ=255) cccGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCa                                                                                             <  2:86572/100‑1 (MQ=255)      gCATATTTATTTCCCGGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCg                                                                                        <  2:129600/100‑1 (MQ=255)      gCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCg                                                                                        >  1:760143/1‑100 (MQ=255)        atatTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTg                                                                                      >  2:613262/1‑100 (MQ=255)        atatTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTg                                                                                      >  2:569462/1‑100 (MQ=255)        atatTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTg                                                                                      <  2:428152/100‑1 (MQ=255)         tatttattTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGa                                                                                     <  1:524311/100‑1 (MQ=255)         tatttattTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGa                                                                                     >  2:524315/1‑100 (MQ=255)         tatttattTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGa                                                                                     >  2:1402728/1‑100 (MQ=255)             tatttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACggg                                                                                 <  1:1242397/100‑1 (MQ=255)             tatttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACggg                                                                                 <  1:803117/100‑1 (MQ=255)             tatttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACgg                                                                                  >  1:279859/1‑99 (MQ=255)              atttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGt                                                                                >  1:1228457/1‑100 (MQ=255)               tttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGtc                                                                               <  1:1781729/100‑1 (MQ=255)               tttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGtc                                                                               <  1:322831/100‑1 (MQ=255)               tttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGtc                                                                               <  2:1459795/100‑1 (MQ=255)                ttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGTATCGACGATTTGACCATCTCGTGACGGGt                                                                                <  1:886868/98‑1 (MQ=255)                 tCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGtcaa                                                                             >  1:658743/1‑100 (MQ=255)                   ccAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGtcaatc                                                                           <  2:1075628/100‑1 (MQ=255)                    cAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGtcaatca                                                                          <  1:438885/100‑1 (MQ=255)                    cAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGtcaatca                                                                          <  2:1611243/100‑1 (MQ=255)                      ggTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATcaac                                                                        <  2:1589085/100‑1 (MQ=255)                      ggTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATcaac                                                                        >  1:905856/1‑100 (MQ=255)                       gTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATcaaca                                                                       <  2:87429/100‑1 (MQ=255)                         aTGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATcaacaa                                                                      >  2:1745515/1‑99 (MQ=255)                           gTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGAGAGGTTAGCATCGACGATTTGACAATCTCGTGACGGGTCAATCAACAACAt                                                                   >  2:1380224/1‑100 (MQ=255)                               cTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGt                                                               <  1:18223/100‑1 (MQ=255)                                 ggCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGt                                                               >  1:897542/1‑98 (MQ=255)                                    acaTCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTACaga                                                          <  2:1279838/100‑1 (MQ=255)                                      aTCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGt                                                        <  2:461316/100‑1 (MQ=255)                                       tCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGtt                                                       <  2:1801341/100‑1 (MQ=255)                                         gTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAg                                                     <  2:760147/100‑1 (MQ=255)                                           ttCATGGTTATAACCATAGTAATTGAGTTCTGCGAGGTTAGAATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTc                                                   <  2:718822/100‑1 (MQ=255)                                           ttCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTc                                                   <  1:757925/100‑1 (MQ=255)                                            tCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCa                                                  <  2:889457/100‑1 (MQ=255)                                            tCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCa                                                  <  2:918607/100‑1 (MQ=255)                                              aTGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGa                                                >  1:1336573/1‑100 (MQ=255)                                                ggTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGaaa                                              >  2:1692747/1‑100 (MQ=255)                                                   tataACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGaaaaaa                                           >  2:1496045/1‑100 (MQ=255)                                                     taACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAgc                                         <  1:1402717/100‑1 (MQ=255)                                                     taACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAgc                                         >  1:1445643/1‑100 (MQ=255)                                                     taACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAgc                                         >  1:680749/1‑100 (MQ=255)                                                     taACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAg                                          >  1:844801/1‑99 (MQ=255)                                                      aaCCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAgcg                                        <  2:463081/100‑1 (MQ=255)                                                        ccATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAgcgcg                                      <  2:1228468/100‑1 (MQ=255)                                                           tAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCAt                                   >  1:1340611/1‑100 (MQ=255)                                                            aGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATa                                  >  1:457272/1‑100 (MQ=255)                                                               aaTTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATaaaa                               >  2:82385/1‑100 (MQ=255)                                                                 ttGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAAc                             >  1:353073/1‑100 (MQ=255)                                                                        cggcgAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGttt                      <  2:767156/100‑1 (MQ=255)                                                                         ggcgAGGTTAGCATCGGCGATTTGACCATCTAGTGACGTGGCAATCAATAACATTGGTGCAGAGTTAGTCAGAAAAAAGTGCTTATAAAAACCTTGttt                      <  1:222001/99‑1 (MQ=255)                                                                         ggcgAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGAATAAAAACCTTGtttt                     <  1:1019741/100‑1 (MQ=255)                                                                           cgAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGtttttt                   <  2:1482761/100‑1 (MQ=255)                                                                            gAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTg                  <  1:569458/100‑1 (MQ=255)                                                                              ggTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTgcg                >  2:473593/1‑100 (MQ=255)                                                                               gTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTgcgc               <  1:1438660/100‑1 (MQ=255)                                                                                ttAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCt              >  1:923624/1‑100 (MQ=255)                                                                                 tAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTg             <  1:791056/100‑1 (MQ=255)                                                                                  aGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGa            <  1:1186854/100‑1 (MQ=255)                                                                                    cATCGACAATTTGACCAGCTCGTGACGGGTCAATGACCAACATTCGTGCAGAAAAAAAAAAAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCt              <  1:1332766/96‑1 (MQ=255)                                                                                     aTCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATaa         >  2:346294/1‑100 (MQ=255)                                                                                       cgacgaTTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAAtt       <  2:1336583/100‑1 (MQ=255)                                                                                        gacgaTTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAAttt      >  1:1592551/1‑100 (MQ=255)                                                                                         acgaTTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTg     <  1:1820643/100‑1 (MQ=255)                                                                                          cgaTTTGACCATCTCGTGACGGGTCAATCAATAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGc    <  2:88203/100‑1 (MQ=255)                                                                                            aTTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCg  <  2:1126534/100‑1 (MQ=255)                                                                                                 |                                                                                               CCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCG  >  NC_000913/1407247‑1407437

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