Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,237,235	C→T	Q173Q (CAG→CAA)	mglC ←	D‑galactose/methyl‑galactoside ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,237,235	0	C	T	100.0%	101.3 / NA	29	Q173Q (CAG→CAA)	mglC	D‑galactose/methyl‑galactoside ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/13);  total (16/13)

CCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGG  >  NC_000913/2236993‑2237476                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  ccGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTGGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCaa                                                                                                                                                                                                                                            >  1:163853/1‑250 (MQ=255)      aCGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGt                                                                                                                                                                                                                                        >  2:275625/1‑249 (MQ=255)                                  gCATAGAACACGCCAGACAACGCGTAGACCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTgg                                                                                                                                                                                                          <  1:230988/251‑1 (MQ=255)                                          cacGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGccc                                                                                                                                                                                                  <  1:275625/251‑1 (MQ=255)                                                 gACAACGCGTAGGTCATAAGCAGGTTCAGGCAGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATTGCAAAAATGTTCTAACCGACGCGGGTTTTGTTCCACAACACCCAGACGAACGCAACCGCAATAAACGCGTAGAAGGAGATGTAAGAGAGCCGGAAACTCCCCAGCGCCACAAAGCCTTGCGCACAGGACGAGAACCCACCGTCAAACCCACAAATTGGCGACGCCCCGACac                                                                                                                                                                                            >  2:227979/1‑249 (MQ=255)                                                                atcaGCAGGTTAAGGCCGCCGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAGGGTAGAGAAGCCAATGACAAAACCAGAAATTGGCGACGCCACGACAAAGACGAAATAGAGCg                                                                                                                                                                             >  2:180830/1‑250 (MQ=255)                                                                                          acCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCg                                                                                                                                                  >  1:110433/1‑251 (MQ=255)                                                                                               attcTTTTGCCGCTCCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCg                                                                                                                                             <  2:71982/248‑1 (MQ=255)                                                                                                    tttGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCg                                                                                                                                         <  1:82616/250‑1 (MQ=255)                                                                                                             cccGGGTTCCCGCCAATGGCAAAAATGTGCTTACCGAAGCGGTGTTTGTTCCACAACACCCAGACGATCGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGTGACGCAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTATAGAAGCCTCTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGGTGAGGCCATTGACGATGATCATCGTGCCGAGCGTGTTAAGGaa                                                                                                                               <  2:269481/250‑1 (MQ=255)                                                                                                                                   aaaTGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGAGATGTCAGAGAGACGGACACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGCTGCTCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCg                                                                                                         >  2:200843/1‑251 (MQ=255)                                                                                                                                         tCTTACCGAAGCGGGTTTTGTTCCACAACAACCAGACGACCGCCACCACCATCAACGCGTAGAAGGTGATGTAAGAGAGACAGAAACCCCCCAGCGCGACAAAGCCTTGAGCAAACGTAGAGAAACCACAGTCACAACCAGAAATTGGCGACGCCCCGACAACGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATAATCGTGCCGCGCGTGGTAATGAACGGCGTCACGATGAGATACGCGATAATc                                                                                                   >  1:286882/1‑251 (MQ=255)                                                                                                                                                         tttGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGACGTAATAGAGCGAGTTGATGCCATAGACGATGATAATCGTGCCGAGCGTGGTAATGCACGGCGTCACGTTGAGAAAAgagagaagaaaaaggaggagaaaa                                                                                    >  2:183680/1‑227 (MQ=255)                                                                                                                                                              tCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAAt                                                                                                                                  <  1:202879/199‑1 (MQ=255)                                                                                                                                                              tCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGAAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAAt                                                                                                                                  >  2:202879/1‑199 (MQ=255)                                                                                                                                                                cacaACACACAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCCAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGAGCAAACCGATCAc                                                                            >  2:251266/1‑251 (MQ=255)                                                                                                                                                                     cacCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCAc                                                                       >  1:277161/1‑251 (MQ=255)                                                                                                                                                                          aGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGTCGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCTCCAATg                                                                  <  1:251266/251‑1 (MQ=255)                                                                                                                                                                              gAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGcac                                                              >  1:283089/1‑251 (MQ=255)                                                                                                                                                                                               aCGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGa                                                           <  1:176110/237‑1 (MQ=255)                                                                                                                                                                                               aCGCGAAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCCCAGa                                                           >  2:176110/1‑237 (MQ=255)                                                                                                                                                                                                               gTAAGAGAGACGTAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGTTAGAGAAGCCGCTGTCAAAACCAGAAATTGGCGACGCCCCGACAAATGCGTAATAGAGCGAGTTGATGCCATAGACGAGGATCAGCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGAGCACCGCACCAATGTCACAGACAATCAGAATAACCAGCGCAATCGGCAtc                              <  2:163853/250‑1 (MQ=255)                                                                                                                                                                                                                   gagagaCGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAACGGTAGAGAAACCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCACTGGCCCAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGc                         >  1:143928/1‑251 (MQ=255)                                                                                                                                                                                                                     gagaCGGAAACTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCa                       <  1:199759/251‑1 (MQ=255)                                                                                                                                                                                                                         cGGAAACTCCCCAGCGCGACAAATCCTTGAGAAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACACAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATAATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATAAGACCGTTGATCAAACCGCTCACCGCACCAATGGCACAGACCATCAGAATACACAGCGCAATCTGCATCGTCGCCATTTc                   >  2:71165/1‑251 (MQ=255)                                                                                                                                                                                                                            aaaCTCCCCAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCgg                >  1:106530/1‑251 (MQ=255)                                                                                                                                                                                                                                cccccAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAac            <  1:7583/250‑1 (MQ=255)                                                                                                                                                                                                                                  cccAGCGCGACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAacac          <  1:78799/251‑1 (MQ=255)                                                                                                                                                                                                                                          gACAAAGCCTTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTgg  <  1:301229/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  CCGATACGTCCGGCTTCTAACATCCCGCCAAAGGCATAGAACACGCCAGACAACGCGTAGATCATCAGCAGGTTCAGGCCGACGTTGACACCAGATACTTTTGCCGCTTCCGGGTTACCGCCAATGGCAAAAATGTTCTTACCGAAGCGGGTTTTGTTCCACAACACCCAGACGAACGCCACCGCAATCAACGCGTAGAAGGTGATGTAAGAGAGACGGAAACTCCCCAGCGCGACAAAGCCCTGAGCAAAGGTAGAGAAGCCACTGTCAAAACCAGAAATTGGCGACGCCCCGACAAAGTCGTAATAGAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGG  >  NC_000913/2236993‑2237476

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