Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,104,024	G→T	I162I (ATC→ATA)	wbbK ←	putative glycosyltransferase WbbK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,104,024	0	G	T	100.0%	97.0 / NA	28	I162I (ATC→ATA)	wbbK	putative glycosyltransferase WbbK
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/13);  total (15/13)

GCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTGATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTT  >  NC_000913/2103782‑2104267                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    gCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAAttt                                                                                                                                                                                                                                             >  1:87762/1‑251 (MQ=255)              cTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCa                                                                                                                                                                                                                                 >  2:622847/1‑250 (MQ=255)                 tACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAAt                                                                                                                                                                                                                             >  2:166660/1‑251 (MQ=255)                   ccACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAAttg                                                                                                                                                                                                                           <  1:604963/251‑1 (MQ=255)                         aTAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGaa                                                                                                                                                                                                                      >  1:502326/1‑250 (MQ=255)                                      aTTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATGATAAGCTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGt                                                                                                                                                                                                         >  2:775282/1‑250 (MQ=255)                                             atTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCattactgtttt                                                                                                                                                                                                 <  2:43419/251‑11 (MQ=255)                                                         tttCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAACTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGtttattt                                                                                                                                                                                     >  1:615869/1‑251 (MQ=255)                                                          ttCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGtttatttt                                                                                                                                                                                    >  1:476539/1‑251 (MQ=255)                                                              aTTTCCTTGCTGCACTAATAATAAGCTTGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATTATAAGTAAGTTGGCTTTTATAAGATAATTTAATTTCTTGCCGACTGACAATGATGTTATTTATAGAAGATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCCGTATTTTTTTTAATGTTTATTTtatat                                                                                                                                                                                <  2:37399/251‑1 (MQ=255)                                                              aTTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTtatat                                                                                                                                                                                >  1:97854/1‑251 (MQ=255)                                                                 tCCTTGCTGAAATAATAAAAAGCTCGTAAATTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAAGTTATAATCATCAACAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTGTAATGTTTATTTTatatat                                                                                                                                                                              <  2:575369/250‑1 (MQ=255)                                                                                 aataaGCTCGTAATTTTTAAATACTCGTTGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGATATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACCCTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCAtt                                                                                                                                                              <  1:168658/250‑1 (MQ=255)                                                                                      gCTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATTTTAAAt                                                                                                                                                        <  1:560075/251‑1 (MQ=255)                                                                                       cTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATTTTAAAt                                                                                                                                                        >  2:25315/1‑250 (MQ=255)                                                                                                   aaaTACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATTTTAAATAAGAAAAAGCTa                                                                                                                                            <  1:775282/250‑1 (MQ=255)                                                                                                             ggAACAGCAGTGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATaa                                                                                                                                  >  1:53689/1‑250 (MQ=255)                                                                                                                                                     tCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAgg                                                                                          <  1:111927/250‑1 (MQ=255)                                                                                                                                                                               aGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTT‑AATGTTTATTTTATATATCAGCCCGTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTa                                                              <  1:125221/251‑1 (MQ=255)                                                                                                                                                                                       tttATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTgac                                                         <  1:587932/249‑1 (MQ=255)                                                                                                                                                                                        ttATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTgacgac                                                      >  1:198042/1‑251 (MQ=255)                                                                                                                                                                                             aGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCA‑AATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGc                                                >  2:570019/1‑251 (MQ=255)                                                                                                                                                                                              gATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGc                                                <  2:476539/251‑1 (MQ=255)                                                                                                                                                                                               aTAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATTTTTAATTTATATATCAGCCCGTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGcc                                               <  2:53689/251‑1 (MQ=255)                                                                                                                                                                                                                               ttatttatAGAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGc               >  1:781737/1‑251 (MQ=255)                                                                                                                                                                                                                                        gAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGctc      <  1:780255/251‑1 (MQ=255)                                                                                                                                                                                                                                        gAATATTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAAAACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATTTTATATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCATTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGctc      >  2:734157/1‑251 (MQ=255)                                                                                                                                                                                                                                            atTTCTTTATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCtttt  >  1:693310/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    GCATACGCATTTTCTGTACCACTGATAGTAAGCAGAAATTTAATATTGGATTGTTCTTTCAATTTCCTTGCTGCACTAATAATAAGCTCGTAATTTTTAAATACTCGTGGAACAGCAGGGTAAAATATTGTCAACTCAGAAGGGTTATTCTTAAATTGAGAATCATCATCAGTAAGTTGGCTTTTATCAGATAATTTAATTTCTGGCCGACTGACAATGATGTTATTTATAGAATATTTCTTGATAAATTTTTCTTTCATCCAGAATTGTTGAACAAACACTGCAGTATTTTTTTTAATGTTTATTTTATATATCAGCCCGTATAGCATTTTAAATAAGAAAAAGCTAGGCTCCATAAGAATTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTT  >  NC_000913/2103782‑2104267

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