Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,678,937	T→A	Q84L (CAG→CTG)	yhjH ←	cyclic‑di‑GMP phosphodiesterase, FlhDC‑regulated

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,678,937	0	T	A	100.0%	255.9 / NA	73	Q84L (CAG→CTG)	yhjH	cyclic‑di‑GMP phosphodiesterase, FlhDC‑regulated
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (32/41);  total (32/41)

CTGACGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACC  >  NC_000913/3678840‑3679032                                                                                                  |                                                                                                cTGACGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGc                                                                                               >  2:413065/1‑100 (MQ=255) cTGACGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGc                                                                                               <  1:400389/100‑1 (MQ=255)   gACGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGctc                                                                                             >  2:1641711/1‑100 (MQ=255)     cGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCtt                                                                                           >  1:1243871/1‑100 (MQ=255)     cGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCtt                                                                                           >  1:768946/1‑100 (MQ=255)        aGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTca                                                                                        <  1:345490/100‑1 (MQ=255)         gggCGATGAGCGTAGGGTCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTcac                                                                                       >  1:1528515/1‑100 (MQ=255)         gggCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTcac                                                                                       >  1:1239405/1‑100 (MQ=255)          ggCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTcaca                                                                                      >  1:1887221/1‑100 (MQ=255)          ggCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTcaca                                                                                      <  1:163810/100‑1 (MQ=255)          ggCGATGAGCGTAGGGCCATAAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCttt                                                                                          <  2:1609566/96‑1 (MQ=255)           gCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACaa                                                                                     <  2:883393/100‑1 (MQ=255)             gATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAAcc                                                                                   >  1:2033769/1‑100 (MQ=255)               tGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTc                                                                                 <  2:1298163/100‑1 (MQ=255)                 aGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCa                                                                               <  2:1424503/100‑1 (MQ=255)                     tAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCg                                                                           >  1:1475719/1‑100 (MQ=255)                     tAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCg                                                                           >  2:1539188/1‑100 (MQ=255)                       gggCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAATTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGAt                                                                         <  2:1001439/100‑1 (MQ=255)                       gggCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGAt                                                                         >  2:1203526/1‑100 (MQ=255)                        ggCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACGACCTCCATCCGATg                                                                        >  2:1391099/1‑100 (MQ=255)                          ccATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGTAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGc                                                                      <  2:1434297/100‑1 (MQ=255)                          ccATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCCCCAGCAAATCAATCAGCTCTTTCACAACTTCCATCCGATgac                                                                      <  1:1382947/100‑3 (MQ=255)                             tCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGa                                                                   >  1:21173/1‑100 (MQ=255)                              cAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGAc                                                                  >  2:372959/1‑100 (MQ=255)                                atatTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACgg                                                                >  1:562270/1‑100 (MQ=255)                                    tGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGAt                                                            <  2:1008259/100‑1 (MQ=255)                                     gACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGAtt                                                           <  2:351816/100‑1 (MQ=255)                                     gACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGAtt                                                           <  1:372239/100‑1 (MQ=255)                                     gACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGAtt                                                           <  1:1391108/100‑1 (MQ=255)                                     gACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGAtt                                                           >  1:168274/1‑100 (MQ=255)                                       ccGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTc                                                         >  1:1989223/1‑100 (MQ=255)                                         gATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAg                                                       <  2:653252/100‑1 (MQ=255)                                         gATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAg                                                       <  2:562268/100‑1 (MQ=255)                                           taccAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTa                                                     <  2:1239399/98‑1 (MQ=255)                                            gCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTaa                                                    <  1:413066/100‑1 (MQ=255)                                             ccAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTaaa                                                   >  1:1664432/1‑100 (MQ=255)                                             ccAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCACTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGaaaa                                                   >  2:1892572/1‑100 (MQ=255)                                                 cagGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAg                                               <  1:117044/100‑1 (MQ=255)                                                  agGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGc                                              >  1:711521/1‑100 (MQ=255)                                                  agGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGc                                              <  2:690405/100‑1 (MQ=255)                                                  agGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGc                                              <  2:833093/100‑1 (MQ=255)                                                   ggCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGc                                              >  2:431072/1‑99 (MQ=255)                                                       gTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATcc                                         <  2:1664419/100‑1 (MQ=255)                                                        tGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCcg                                        >  2:431321/1‑100 (MQ=255)                                                             tataAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCggcggc                                   <  2:1562374/100‑1 (MQ=255)                                                             tataAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCggcggc                                   <  1:431322/100‑1 (MQ=255)                                                               taAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAg                                 >  1:670751/1‑100 (MQ=255)                                                               taAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAg                                 <  1:1657708/100‑1 (MQ=255)                                                                naagaagTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAgg                                >  2:951360/2‑100 (MQ=255)                                                                  agaagTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCCGGCg                              >  2:727127/1‑100 (MQ=255)                                                                    aagTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGtt                            <  2:96727/100‑1 (MQ=255)                                                                      gTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGc                          >  2:1038798/1‑100 (MQ=255)                                                                       tCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCg                         >  1:78188/1‑100 (MQ=255)                                                                       tCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCg                         <  2:1528506/100‑1 (MQ=255)                                                                       tCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCg                         <  2:1059138/100‑1 (MQ=255)                                                                        cGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGa                        >  2:1658573/1‑100 (MQ=255)                                                                         ggCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGaa                       >  1:874475/1‑100 (MQ=255)                                                                         ggCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGaa                       >  2:659396/1‑100 (MQ=255)                                                                          gCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAg                      <  2:670749/100‑1 (MQ=255)                                                                           cTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAgg                     <  2:853250/100‑1 (MQ=255)                                                                             ttttGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGt                   <  2:564319/100‑1 (MQ=255)                                                                                tGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCa                <  2:78187/100‑1 (MQ=255)                                                                                  cgcCAGCAAATCACTCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCaan              <  1:951364/100‑2 (MQ=255)                                                                                  cgcCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAg              <  1:1539197/100‑1 (MQ=255)                                                                                   gcCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAgg             <  2:1767513/100‑1 (MQ=255)                                                                                       gCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCTGCAGGCGTTGCGAAGGGTTCAAGGGATg         <  1:1707105/100‑1 (MQ=255)                                                                                        cAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGc        >  1:916308/1‑100 (MQ=255)                                                                                        cAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGc        <  1:1203532/100‑1 (MQ=255)                                                                                            tcaatcAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGa    <  2:1243865/100‑1 (MQ=255)                                                                                             caatcAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGac   >  2:499546/1‑100 (MQ=255)                                                                                             caatcAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGac   <  1:68361/100‑1 (MQ=255)                                                                                              aatcAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCCGGCGTTGCGAAGGGTTCAAGGGATGCGTGacc  <  1:1933571/100‑1 (MQ=255)                                                                                              aatcAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGacc  <  1:1658586/100‑1 (MQ=255)                                                                                                  |                                                                                                CTGACGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACC  >  NC_000913/3678840‑3679032

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