Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	478,533	C→T	A29T (GCC→ACC)	ylaB ←	putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	478,533	0	C	T	100.0%	197.9 / NA	61	A29T (GCC→ACC)	ylaB	putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/34);  total (27/34)

AGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGCCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGT  >  NC_000913/478402‑478651                                                                                                                                    |                                                                                                                       aGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGtcagccag                                                                                                                 >  1:625017/1‑139 (MQ=255) aGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGtcagccag                                                                                                                 >  1:442987/1‑139 (MQ=255) aGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGtcagccag                                                                                                                 >  2:332215/1‑139 (MQ=255)  gcgcCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGa                                                                                                                >  2:30432/1‑139 (MQ=255)                ttCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACa                                                                                                  <  1:766645/139‑1 (MQ=255)                        ccTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCaa                                                                                          <  2:568796/139‑1 (MQ=255)                         cTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGTCAGGCAGCAAt                                                                                         >  2:593263/1‑139 (MQ=255)                          tGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATa                                                                                        >  2:479130/1‑139 (MQ=255)                               cgcATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGaa                                                                                   <  1:4496/139‑1 (MQ=255)                                                tCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                             >  2:780425/1‑112 (MQ=255)                                                tCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                             <  1:780425/112‑1 (MQ=255)                                                                tCGTATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCATATGCTTAAGCCGACAGGTAGCAATACTGAAAGAATCAGTACTcc                                                                     >  1:517089/1‑120 (MQ=255)                                                                tCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACa                                                                                                  <  1:719227/91‑1 (MQ=255)                                                                tCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACa                                                                                                  >  2:719227/1‑91 (MQ=255)                                                                tCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTcc                                                                     <  2:517089/120‑1 (MQ=255)                                                                  gAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATg                                                >  1:319156/1‑139 (MQ=255)                                                                  gAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATg                                                <  1:227465/139‑1 (MQ=255)                                                                   aaTAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGaa                                                                                   <  2:93686/103‑1 (MQ=255)                                                                   aaTAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGaa                                                                                   >  1:93686/1‑103 (MQ=255)                                                                   aaTAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGt                                               >  2:428297/1‑139 (MQ=255)                                                                        cGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCgtgt                                          >  2:645489/1‑139 (MQ=255)                                                                            gcgGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTtctc                                      <  1:428297/139‑1 (MQ=255)                                                                             cgGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAgg                                                                                                >  2:145024/1‑80 (MQ=255)                                                                             cgGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAgg                                                                                                <  1:145024/80‑1 (MQ=255)                                                                             cgGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCa                                     >  2:205720/1‑139 (MQ=255)                                                                                       aGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATc                           <  1:30432/139‑1 (MQ=255)                                                                                       aGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATc                           <  2:243531/139‑1 (MQ=255)                                                                                       aGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATc                           <  2:25548/139‑1 (MQ=255)                                                                                         gTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGtt                                         <  2:585431/123‑1 (MQ=255)                                                                                         gTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGtt                                         >  1:585431/1‑123 (MQ=255)                                                                                             ccAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTa                                                                         <  2:119472/87‑1 (MQ=255)                                                                                             ccAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTa                                                                         >  1:119472/1‑87 (MQ=255)                                                                                             ccAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCa                               <  1:735020/129‑1 (MQ=255)                                                                                             ccAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCa                               >  2:735020/1‑129 (MQ=255)                                                                                              cAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCttatta                    >  2:608368/1‑139 (MQ=255)                                                                                                gCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAac                  >  2:37768/1‑139 (MQ=255)                                                                                                 ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                 <  2:552460/139‑1 (MQ=255)                                                                                                 ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                 <  2:608266/139‑1 (MQ=255)                                                                                                 ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                 <  1:410355/139‑1 (MQ=255)                                                                                                 ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                 <  1:479130/139‑1 (MQ=255)                                                                                                 ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                 <  1:538930/139‑1 (MQ=255)                                                                                                 ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                 <  2:737844/139‑1 (MQ=255)                                                                                                 ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                 <  2:92025/139‑1 (MQ=255)                                                                                                 ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                 >  1:666621/1‑139 (MQ=255)                                                                                                  tctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAacac                <  1:622107/139‑1 (MQ=255)                                                                                                    ttCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAg              <  1:205720/139‑1 (MQ=255)                                                                                                      cAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTtctc                                      <  1:357249/113‑1 (MQ=255)                                                                                                      cAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTtctc                                      >  2:357249/1‑113 (MQ=255)                                                                                                          aaaCGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAgg        <  2:666621/139‑1 (MQ=255)                                                                                                            aCGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGAt                                                 >  2:190555/1‑96 (MQ=255)                                                                                                            aCGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGAt                                                 <  1:190555/96‑1 (MQ=255)                                                                                                            aCGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGa      <  2:442987/139‑1 (MQ=255)                                                                                                            aCGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGa      <  1:593263/139‑1 (MQ=255)                                                                                                              gATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGaaa    >  2:54650/1‑139 (MQ=255)                                                                                                                tGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCata                             >  1:8032/1‑112 (MQ=255)                                                                                                                tGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCata                             <  2:8032/112‑1 (MQ=255)                                                                                                                tGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGt  >  1:216097/1‑139 (MQ=255)                                                                                                                tGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGt  >  1:36955/1‑139 (MQ=255)                                                                                                                       aCCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                 <  2:753410/117‑1 (MQ=255)                                                                                                                       aCCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                 >  1:753410/1‑117 (MQ=255)                                                                                                                           gctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAgg        <  1:388355/122‑1 (MQ=255)                                                                                                                           gctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAgg        >  2:388355/1‑122 (MQ=255)                                                                                                                                tGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAg                                                                           <  2:15914/50‑1 (MQ=255)                                                                                                                                tGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAg                                                                           >  1:15914/1‑50 (MQ=255)                                                                                                                                    |                                                                                                                       AGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGCCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGT  >  NC_000913/478402‑478651

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