Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,933,154	G→T	Q179K (CAA→AAA)	fucA ←	L‑fuculose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,933,154	0	G	T	100.0%	70.5 / NA	21	Q179K (CAA→AAA)	fucA	L‑fuculose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (14/7);  total (14/7)

GTTGCTTTACGAAATTACTCTTCAATTCGTAACCCATAGGTTTTGAATTTCTCCAGCACTACGGCAATCTCTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCTTACGATTTTTGAGAGCCAGCGCAACATGTTCAGA  >  NC_000913/2933027‑2933280                                                                                                                                |                                                                                                                               gttgCTTTACGAAATTACTCTTCAATTCGTAACCCATAGGTTTTGAATTTCTCCAGCACTACGGCAATCTCTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCacttc                                                                                                                    >  2:198161/1‑140 (MQ=255)          cGAAATTACTCTTCAATTCGTAACCCATAGGTTTTGAATTTCTCCAGCACTACGGCAATCTCTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATg                                                                                                           >  1:353420/1‑140 (MQ=255)                       cAATTCGTAACCCATAGGTTTTGAATTTCTCCAGCACTACGGCAATCTCTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAAc                                                                                              >  1:298894/1‑140 (MQ=255)                          ttCGTAACCCATAGGTTTTGAATTTCTCCAGCACTACGGCAATCTCTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCt                                                                                           >  2:264240/1‑140 (MQ=255)                                            tGAATTTCTCTAGCACTACGGCAATCTCTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTca                                                                         >  1:80262/1‑140 (MQ=255)                                            tGAATTTCTCCAGCACTACGGCAATCTCTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTca                                                                         >  1:252364/1‑140 (MQ=255)                                            tGAATTTCTCCAGCACTACGGCAATCTCTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTca                                                                         >  2:301651/1‑140 (MQ=255)                                            tGAATTTCTCCAGCACTACGGCAATCTCTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTca                                                                         <  2:306011/140‑1 (MQ=255)                                                    tcCAGCACTACGGCAATCTCTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGAt                                                                 <  2:208721/140‑1 (MQ=255)                                                     ccAGCACTACGGCAATCTCTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATa                                                                <  2:101241/140‑1 (MQ=255)                                                     ccAGCACTACGGCAATCTCTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATa                                                                <  2:106691/140‑1 (MQ=255)                                                        gCACTACGGCAATCTCTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGc                                                             <  2:320428/140‑1 (MQ=255)                                                                   aTCTCTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGAtgttg                                                  >  1:131659/1‑140 (MQ=255)                                                                    tctcTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGAtgttgt                                                 >  2:61590/1‑140 (MQ=255)                                                                          cATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACaaa                                           >  1:118306/1‑140 (MQ=255)                                                                                              gggTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCTTACGATTTTTgag                       >  2:91158/1‑140 (MQ=255)                                                                                                     tAATCGCCAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCTTACGATTTTTGAGAGCCAgc                <  2:142254/140‑1 (MQ=255)                                                                                                           ccAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCTTACGATTTTTGAGAGc                    >  2:348960/1‑130 (MQ=255)                                                                                                           ccAGGGTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCTTACGATTTTTGAGAGc                    <  1:348960/130‑1 (MQ=255)                                                                                                              gggTCGTCAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCTTACGATTTTTGAGAGCCAGCGCAACATGt       >  1:109334/1‑140 (MQ=255)                                                                                                                   gtcAGGTAAAGTTTCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCTTACGATTTTTGAGAGCCAGCGCAACATGTTCAGa  >  2:71741/1‑140 (MQ=255)                                                                                                                                |                                                                                                                               GTTGCTTTACGAAATTACTCTTCAATTCGTAACCCATAGGTTTTGAATTTCTCCAGCACTACGGCAATCTCTTCATCGCTCAGCACTGGCACCGGGTCCGTAATCGCCAGGGTCGTCAGGTAAAGTTGCGCCAGCACTTCAACTTCATGCGCCAGCCATAACGCTTTTTCCAGATTCACCTCACAAGCGATAAGCCCATGATGTTGTAACAAAGTTGCCTTACGATTTTTGAGAGCCAGCGCAACATGTTCAGA  >  NC_000913/2933027‑2933280

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: ‑