Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,878,732	Δ2 bp	coding (1018‑1019/1026 nt)	yeaP →	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,878,731	0	A	.	100.0%	146.3 / NA	31	coding (1017/1026 nt)	yeaP	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (10/21);  total (10/21)
*	NC_000913	1,878,732	0	C	.	100.0%	146.3 / NA	31	coding (1018/1026 nt)	yeaP	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (10/21);  total (10/21)

ATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGAT  >  NC_000913/1878583‑1878875                                                                                                                                                     ||                                                                                                                                                aTTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTatt                                                                                                                                                >  1:146367/1‑149 (MQ=255)     gTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTattcct                                                                                                                                             <  1:404215/150‑6 (MQ=255)                  gCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGTCGTATTcac                                                                                                                              <  2:611672/150‑1 (MQ=255)                           ggCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGa                                                                                                                    <  2:213998/151‑1 (MQ=255)                           ggCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGa                                                                                                                    <  2:199581/151‑1 (MQ=255)                                agtagAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTgg                                                                                                                >  2:621412/1‑150 (MQ=255)                                agtagAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTgg                                                                                                                <  1:244681/150‑1 (MQ=255)                                agtagAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGt                                                                                                               >  2:177208/1‑151 (MQ=255)                                       gTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACa                                                                                                        <  2:606782/151‑1 (MQ=255)                                       gTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACa                                                                                                        <  2:309907/151‑1 (MQ=255)                                         cGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACAt                                                                                                       <  1:152038/150‑1 (MQ=255)                                            tCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCg                                                                                                    <  2:493415/150‑1 (MQ=255)                                                      acCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGtatta                                                                                         <  1:60167/151‑1 (MQ=255)                                                            gCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTAttttt                                                                                    <  1:177213/150‑1 (MQ=255)                                                             cAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTcc                                                                                  >  2:344671/1‑151 (MQ=255)                                                              aGACAGTTCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACTCCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGTTATTAttttt                                                                                    <  2:285432/148‑1 (MQ=255)                                                                                    gATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCgg                                                           >  1:141438/1‑151 (MQ=255)                                                                                       atTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTc                                                         <  1:288976/150‑1 (MQ=255)                                                                                                 aTCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGc                                              <  1:104230/151‑1 (MQ=255)                                                                                                  tCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGc                                              >  1:312854/1‑150 (MQ=255)                                                                                                   cAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCAt                                            <  1:344678/151‑1 (MQ=255)                                                                                                        gaAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGttt                                       >  1:143759/1‑151 (MQ=255)                                                                                                                acaAACAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTgg                               <  1:236467/151‑1 (MQ=255)                                                                                                                    aCAGAAAACACCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTggcggc                           <  2:297762/151‑1 (MQ=255)                                                                                                                     cAGAAAACCCCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCt                          <  1:71777/151‑1 (MQ=255)                                                                                                                             acCTTTTGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGcc                  >  1:89918/1‑151 (MQ=255)                                                                                                                                 tttGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCgg              >  1:124395/1‑151 (MQ=255)                                                                                                                                   tGTCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCggg             >  1:331303/1‑150 (MQ=255)                                                                                                                                     tCGCGCATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGc          >  2:562789/1‑151 (MQ=255)                                                                                                                                          cATCCAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATaaa     <  1:437680/151‑1 (MQ=255)                                                                                                                                             ccAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGAt  <  1:268136/151‑1 (MQ=255)                                                                                                                                             ccAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGAt  <  1:159865/151‑1 (MQ=255)                                                                                                                                             ccAGCGCT‑‑ATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGAt  <  2:70548/151‑1 (MQ=255)                                                                                                                                                     ||                                                                                                                                                ATTTGTATTATCCCGGTGCCAGTCTTGGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGAT  >  NC_000913/1878583‑1878875

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