Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,610,848	C→T	intergenic (‑168/+59)	uxaB ← / ← dgcF	tagaturonate reductase/putative diguanylate cyclase DgcF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,610,848	0	C	T	100.0%	89.4 / NA	25	intergenic (‑168/+59)	uxaB/dgcF	tagaturonate reductase/putative diguanylate cyclase DgcF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/10);  total (15/10)

CTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCA  >  NC_000913/1610775‑1610922                                                                          |                                                                           cTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGa                                                                       <  1:301982/79‑1 (MQ=255)   gatgatTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCa                                                                    >  2:302007/1‑80 (MQ=255)       atTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATaaa                                                                <  2:301983/80‑1 (MQ=255)         tAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTg                                                             >  2:302010/1‑81 (MQ=255)         tAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTg                                                             >  1:302016/1‑81 (MQ=255)              aGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAg                                                         <  1:301984/80‑1 (MQ=255)              aGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGa                                                        <  1:301985/81‑1 (MQ=255)                      tATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACAtt                                                 <  2:301986/80‑1 (MQ=255)                        tGAGAAAAGTGTGGCGCGGGTCATGGTTTAATCGAGGAAAAACGCCTTTTTCTAGATCATAGAGTGGTAGAACACATTGCa                                              <  2:301987/81‑1 (MQ=255)                         gagaAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCAt                                             >  1:302024/1‑81 (MQ=255)                                 tgtgGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATcg                                     <  2:301988/81‑1 (MQ=255)                                  gtgGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCTTAAAGTGGTAGAACACATTGCATTCAAATcgc                                    >  1:302013/1‑81 (MQ=255)                                     gcgcGGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATcgcgc                                  <  1:301989/80‑1 (MQ=255)                                       gcgGATCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGt                                >  1:302012/1‑80 (MQ=255)                                            tCATGGTTTAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGaa                          >  1:302029/1‑81 (MQ=255)                                                    tAATCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATg                  >  2:302020/1‑81 (MQ=255)                                                     aaTCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGt                 >  1:302032/1‑81 (MQ=255)                                                      aTCGAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGt                 >  2:302014/1‑80 (MQ=255)                                                         gAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGa             >  2:302027/1‑81 (MQ=255)                                                         gAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGa             >  2:302025/1‑81 (MQ=255)                                                         gAGGAAAAACGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGa             >  2:302021/1‑81 (MQ=255)                                                               aaaCGCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACtt       >  1:302030/1‑81 (MQ=255)                                                                   gCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCt    >  1:302017/1‑80 (MQ=255)                                                                   gCCTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTc   <  2:301990/81‑1 (MQ=255)                                                                    ccTTTTTCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCa  <  2:301991/81‑1 (MQ=255)                                                                          |                                                                           CTGATGATTAATGAGCAGTTTTATGAGAAAAGTGTGGCGCGGATCATGGTTTAATCGAGGAAAAACGCCTTTTCCTGGATCATAAAGTGGTAGAACACATTGCATTCAAATCGCGCGTAATGAATAAAGATGTCAGACAACTTCCTCA  >  NC_000913/1610775‑1610922

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