Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,877,653	T→C	intergenic (‑121/‑62)	yoaF ← / → dgcP	DUF333 domain‑containing lipoprotein YoaF/diguanylate cyclase DgcP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,877,653	0	T	C	100.0%	80.6 / NA	23	intergenic (‑121/‑62)	yoaF/dgcP	DUF333 domain‑containing lipoprotein YoaF/diguanylate cyclase DgcP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/6);  total (17/6)

GATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACA  >  NC_000913/1877521‑1877788                                                                                                                                     |                                                                                                                                        gATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCa                                                                                                                                  <  1:700972/140‑1 (MQ=255)        ttCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGc                                                                                                                           <  1:1025586/140‑1 (MQ=255)        ttCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGc                                                                                                                           <  2:877227/140‑1 (MQ=255)        ttCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGc                                                                                                                           <  2:60281/140‑1 (MQ=255)         tCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCa                                                                                                                          <  1:338185/140‑1 (MQ=255)                      tcACATATTTATCCTAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCAAAGATTAaataat                                                                                                             <  1:802095/140‑1 (MQ=255)                              ttATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCGCTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCAAAGGTTAAATAATCAGCAttt                                                                                                     >  2:961196/1‑140 (MQ=255)                              ttATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCAAAGATTAAATAATCAGCAttt                                                                                                     >  2:216006/1‑140 (MQ=255)                              ttATCCAAGATTAGAGTATCGCGGTATCGTTTTGTCTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCAAAGATTAAATAATCAGCAttt                                                                                                     >  2:891030/1‑140 (MQ=255)                                aTCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTAc                                                                                                   >  2:675714/1‑140 (MQ=255)                                      gATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACccgccg                                                                                             >  1:1066096/1‑140 (MQ=255)                                                tCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAg                                                                                   >  1:341029/1‑140 (MQ=255)                                                tCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAg                                                                                   >  1:31270/1‑140 (MQ=255)                                                tCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAg                                                                                   >  1:898977/1‑140 (MQ=255)                                                                   tgCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGtcagatcag                                                                >  1:1045807/1‑140 (MQ=255)                                                                                       aTTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCa                                            >  2:719213/1‑140 (MQ=255)                                                                                                     tatTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAg                              >  1:52562/1‑140 (MQ=255)                                                                                                      atTACGTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCAAAGATTAAATAAACAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGa                             >  2:745867/1‑140 (MQ=255)                                                                                                           gTCTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGt                        >  2:538617/1‑140 (MQ=255)                                                                                                             cTTGTTTCATCTTTGTTGATGATGCTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCa                      >  1:492513/1‑140 (MQ=255)                                                                                                                 ttttttttttGTTGATTTTTTTTTATCACGCCTAAAAAAATAAAATAATCAGCTTTCCCCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTgg                  <  2:52562/133‑1 (MQ=255)                                                                                                                      tCTTTGTTGATGATGCTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGt             >  1:701323/1‑140 (MQ=255)                                                                                                                           gttgATGATGCTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGt        >  1:546920/1‑140 (MQ=255)                                                                                                                                 gatgCTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACa  >  1:21922/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                        GATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACA  >  NC_000913/1877521‑1877788

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