Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,094,491	T→C	intergenic (‑219/+232)	rhaD ← / ← rhaA	rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,094,491	0	T	C	100.0%	88.0 / NA	32	intergenic (‑219/+232)	rhaD/rhaA	rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/16);  total (16/16)

GACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATC  >  NC_000913/4094366‑4094621                                                                                                                              |                                                                                                                                   gACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAgg                                                                                                                       >  2:709232/1‑139 (MQ=21) gACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAgg                                                                                                                       <  1:709232/139‑1 (MQ=25) gACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAgg                                                                                                                       <  1:884497/139‑1 (MQ=21)                             aCTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCgg                                                                                           <  1:538851/139‑1 (MQ=255)                                aTTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaa                                                                                        >  2:466139/1‑139 (MQ=37)                                aTTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGc                                                                                          >  1:567581/1‑137 (MQ=37)                                aTTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGc                                                                                          <  2:567581/137‑1 (MQ=25)                                 ttGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaat                                                                                       <  2:446129/139‑1 (MQ=25)                                   gCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatca                                                                                     <  2:367343/139‑1 (MQ=21)                                        ccTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGcc                                                                                >  2:780506/1‑139 (MQ=37)                                         cTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCt                                                                               >  1:446129/1‑139 (MQ=21)                                                 cTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGAc                                                                       <  1:34971/139‑1 (MQ=18)                                                  tGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACg                                                                      >  2:360075/1‑139 (MQ=18)                                                           aCTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGt                                                             <  1:360075/139‑1 (MQ=18)                                                               gtgCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTa                                                         >  2:155539/1‑139 (MQ=18)                                                                  ccGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCAATTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCa                                                      >  2:669776/1‑139 (MQ=17)                                                                    gCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAgg                                                    <  1:466139/139‑1 (MQ=18)                                                                     cATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                                                   >  2:1013982/1‑139 (MQ=18)                                                                     cATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGc                                                   >  1:919502/1‑139 (MQ=18)                                                                       tCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAg                                                                                                              <  1:697726/78‑1 (MQ=21)                                                                       tCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAg                                                                                                              >  2:697726/1‑78 (MQ=21)                                                                                   tGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGcc                                     <  1:623945/139‑1 (MQ=37)                                                                                         aTGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTg                               >  2:623945/1‑139 (MQ=37)                                                                                          tGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGt                              >  2:1027680/1‑139 (MQ=255)                                                                                              aCGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAc                                            <  1:835828/121‑1 (MQ=17)                                                                                              aCGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAc                                            >  2:835828/1‑121 (MQ=18)                                                                                              aCGCTGTCGCGTCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGc                          <  1:780506/139‑1 (MQ=255)                                                                                                         tCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCt                                                 <  2:356545/105‑1 (MQ=18)                                                                                                         tCTTATCAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCt                                                 >  1:356545/1‑105 (MQ=18)                                                                                                               cAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                                                  >  2:395123/1‑82 (MQ=25)                                                                                                               cAGGCCTACAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                                                  <  1:395123/82‑1 (MQ=18)                                                                                                                      aCAACTGTCGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATc  <  1:155539/139‑1 (MQ=38)                                                                                                                              |                                                                                                                                   GACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATC  >  NC_000913/4094366‑4094621

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