Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,900,687	C→T	E512K (GAA→AAA)	bglH ←	carbohydrate‑specific outer membrane porin, cryptic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,900,687	0	C	T	100.0%	113.9 / NA	34	E512K (GAA→AAA)	bglH	carbohydrate‑specific outer membrane porin, cryptic
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/17);  total (17/17)

ACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTATATTGGTCCCAAATATAGGCCGGGCGAACAACCGCACGAATAGATTCGAAATCAGAGTG  >  NC_000913/3900450‑3900921                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           aCCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGat                                                                                                                                                                                                                                >  1:478127/1‑250 (MQ=255)          acaTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCg                                                                                                                                                                                                                      <  1:2133/251‑1 (MQ=255)                       ccGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATAAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAAttt                                                                                                                                                                                                         <  1:544557/251‑1 (MQ=255)                                  ggTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAg                                                                                                                                                                                              <  2:198635/251‑1 (MQ=255)                                       cGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCa                                                                                                                                                                                          >  1:59556/1‑250 (MQ=255)                                             ggctggcATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATa                                                                                                                                                                                     <  1:272173/249‑1 (MQ=255)                                                                    cccACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAAGTTCCGGACGCTAGGTCAACATTCTGGTATTGACTTTAAAGGTATGGa                                                                                                                                                            <  2:525336/251‑1 (MQ=255)                                                                          gAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGc                                                                                                                                                       <  1:576830/250‑1 (MQ=255)                                                                              aGAAATACAGCTGGCCATAACCAGCGTTAAAGCAGCAATTATTATATTCATCGACATATTCATTTCATTAAGTTCAATTACCACCAGAGTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAACGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATCCTGGTATTGACTTTAACGGTAAGGaatatcgttgt                                                                                                                                                  >  2:292382/1‑242 (MQ=255)                                                                                   tACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTACTATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTtat                                                                                                                                             >  1:604224/1‑251 (MQ=255)                                                                                                 gCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGTGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTa                                                                                                                               <  2:59556/251‑1 (MQ=255)                                                                                                         tAAAGCAGCAATTTTTATATTCATCGTAATATTCATTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGGGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATttatt                                                                                                                       >  1:205497/1‑251 (MQ=255)                                                                                                              cagcaATTTTTATATTCATCGTAATAGTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTTCTAt                                                                                                                  <  2:184039/251‑1 (MQ=255)                                                                                                               agcaATTTTTATAATCATCGTTATATTCCTTTCAGTAATTTCCATTACCGCCAGATTTCAGCCTGGGCACCGATAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGTCATCCAGGTCGTTTTTCAGGGCTTTGATATGCGTGGCGTAGAAACGAATTTCCGGGCGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGATTCATTAAATTTATTACTAtt                                                                                                                 <  2:356418/251‑1 (MQ=255)                                                                                                                                  gTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTgg                                                                                              >  1:324517/1‑251 (MQ=255)                                                                                                                                                   aaggtCCAGTACCACCAGATTTCAGCCTGTTCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCAGCCAGTGCGTGTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGAGATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATGAAAATTATTACGATTCGCATCTTTTTATTGCAGGGTGAATTAGCCCAggtc                                                                             <  2:388988/247‑4 (MQ=255)                                                                                                                                                      ttCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAac                                                                          >  2:503860/1‑251 (MQ=255)                                                                                                                                                                 ccaGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGt                                                               >  2:241604/1‑251 (MQ=255)                                                                                                                                                                     aTTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTtata                                                           <  2:201552/251‑1 (MQ=255)                                                                                                                                                                       ttCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTatat                                                          <  1:324091/250‑1 (MQ=255)                                                                                                                                                                             cTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTATATTGGTccc                                                   >  2:511942/1‑251 (MQ=255)                                                                                                                                                                              tGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTATATTGGTccc                                                   >  2:424966/1‑250 (MQ=255)                                                                                                                                                                                ggCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTATATTGGTCCCaaa                                                <  2:21638/251‑1 (MQ=255)                                                                                                                                                                                ggCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTATATTGGTCCCaa                                                 >  1:362746/1‑250 (MQ=255)                                                                                                                                                                                 gCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTATATTGGTCCCAAat                                               <  1:298055/251‑1 (MQ=255)                                                                                                                                                                                  cACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTATATTGGTCCCAAata                                              <  1:581024/251‑1 (MQ=255)                                                                                                                                                                                                 ggTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCCGTTCGTTTTTCAGGGCTTTAATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTATATTGGTCCCAAATATAGGCCGGGCGaac                                >  1:260003/1‑250 (MQ=255)                                                                                                                                                                                                    cgtcTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTATATTGGTCCCAAATATAGGCCGGGCGAACAAcc                            >  1:428048/1‑251 (MQ=255)                                                                                                                                                                                                             tGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCACCACCTGTCTGGTTATATTGGTCCCAAATATAGGCAGGGCGAACAACCGCACGAAt                    >  2:25462/1‑250 (MQ=255)                                                                                                                                                                                                                     aaGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATATGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTATATTGGTCCCa                                                  >  1:396586/1‑212 (MQ=255)                                                                                                                                                                                                                     aaGGTGAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTATATTGGTCCCa                                                  <  2:396586/212‑1 (MQ=255)                                                                                                                                                                                                                          gAAGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATTCTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTATATTGGTCCCAAATATAGGCCGGGCGAACAACCGCACGAATAGATTCGAAATCa       <  2:486735/250‑1 (MQ=255)                                                                                                                                                                                                                           aaGCCATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTATATTGGTACCAAATATAGGCCGGGCGAACAACCGCACGAATAGATTCGAAATc        >  2:190679/1‑248 (MQ=255)                                                                                                                                                                                                                              ccATCCAGTTCGTTTTTCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTATATTGGTCCCAAATATAGGCCGGGCGAACAACCGCACGAATAGATTCGAAATCAGAGTg  >  2:408182/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           ACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGTCGTCTTTATTGTCTTCAAAGGTGAAGCCATCCAGTTCGTTTTCCAGGGCTTTGATATACGTGGCGTAGAAACGAATTTCCGGACGCGAGGTCAACATACTGGTATTGACTTTAAAGGTATGGAAGAGCGTGGTTTTATAACCAGACTCATTAAATTTATTACTATTCGCATCTTTGTTTTGCTGGGTGAAATAGCCCAGTTCAACACCTGTCTGGTTATATTGGTCCCAAATATAGGCCGGGCGAACAACCGCACGAATAGATTCGAAATCAGAGTG  >  NC_000913/3900450‑3900921

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