Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,494,950	T→C	Y85H (TAT→CAT)	ypdI →	putative lipoprotein involved in colanic acid biosynthesis

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,494,950	0	T	C	100.0%	181.8 / NA	51	Y85H (TAT→CAT)	ypdI	putative lipoprotein involved in colanic acid biosynthesis
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/25);  total (26/25)

AAAAACATCAGCGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCC  >  NC_000913/2494814‑2495068                                                                                                                                         |                                                                                                                       aaaaaCATCAGCGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCAt                                                                                                                      <  2:1022967/139‑1 (MQ=255) aaaaaCATCAGCGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCAt                                                                                                                      <  1:467807/139‑1 (MQ=255)                    acCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCAt                                                                                                                      >  1:51657/1‑120 (MQ=255)                     ccACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATa                                                                                                  <  2:890/139‑1 (MQ=255)                      cacacaCTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATaa                                                                                                 >  2:790308/1‑138 (MQ=255)                      cacacaCTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATaa                                                                                                 <  2:852537/139‑2 (MQ=255)                             tctACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCttt                                                                                          >  1:894000/1‑139 (MQ=255)                               tACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCAt                                                                                                                      <  2:554130/109‑1 (MQ=255)                               tACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCAt                                                                                                                      >  1:554130/1‑109 (MQ=255)                               tACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTca                                                                                        >  1:276724/1‑139 (MQ=255)                               tACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTca                                                                                        >  2:608721/1‑139 (MQ=255)                                 cAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCAt                                                                                                                      >  2:585452/1‑107 (MQ=255)                                  aaaaCTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACAt                                                                                     <  1:608721/139‑1 (MQ=255)                                  aaaaCTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACAt                                                                                     <  1:790308/139‑1 (MQ=255)                                  aaaaCTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACAt                                                                                     <  2:16286/139‑1 (MQ=255)                                     aCTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGc                                                                                  <  2:276724/139‑1 (MQ=255)                                       ttCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCAt                                                                                                                      >  2:603362/1‑101 (MQ=255)                                       ttCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                >  1:560274/1‑139 (MQ=255)                                       ttCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                >  1:657013/1‑139 (MQ=255)                                       ttCAGCCATCACAGGCGATTTTATTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                >  1:262120/1‑139 (MQ=255)                                         cAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGt                                                                              >  2:624309/1‑139 (MQ=255)                                                 acaGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGGAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATc                                                                      <  2:787483/139‑1 (MQ=255)                                                   aGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCAt                                                                    <  1:542486/139‑1 (MQ=255)                                                         ttttGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGc                                                                                  >  1:33762/1‑119 (MQ=255)                                                         ttttGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGc                                                                                  <  2:33762/119‑1 (MQ=255)                                                             atttAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgct                                                          <  2:262120/138‑1 (MQ=255)                                                                 aaGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTAtt                                                                                                                <  2:614319/81‑1 (MQ=255)                                                                 aaGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTAtt                                                                                                                >  1:614319/1‑81 (MQ=255)                                                                      gagagaAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGa                                                                                                    <  2:27283/88‑1 (MQ=255)                                                                      gagagaAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGa                                                                                                    >  1:27283/1‑88 (MQ=255)                                                                        gagaAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACAt                                                                                     >  2:241268/1‑101 (MQ=255)                                                                        gagaAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACAt                                                                                     <  1:241268/101‑1 (MQ=255)                                                                             gaTGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACAt                                                                                     >  1:232558/1‑96 (MQ=255)                                                                             gaTGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACAt                                                                                     <  2:232558/96‑1 (MQ=255)                                                                                gcgcCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCaa                                       <  2:884993/139‑1 (MQ=255)                                                                                              ccTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGa                         >  2:549589/1‑139 (MQ=255)                                                                                                 cacaGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAac                                                 >  2:97537/1‑112 (MQ=255)                                                                                                 cacaGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAac                                                 <  1:97537/112‑1 (MQ=255)                                                                                                        gaaATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCg               <  1:1010502/139‑1 (MQ=255)                                                                                                        gaaATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCg               >  2:184352/1‑139 (MQ=255)                                                                                                           aTGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCt            >  2:827873/1‑139 (MQ=255)                                                                                                            tGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTg           >  1:72068/1‑139 (MQ=255)                                                                                                                gccCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATc                                                                      <  1:523294/76‑1 (MQ=255)                                                                                                                gccCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATc                                                                      >  2:523294/1‑76 (MQ=255)                                                                                                                gccCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCAt       >  2:382610/1‑139 (MQ=255)                                                                                                                gccCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCAt       >  2:569195/1‑139 (MQ=255)                                                                                                                   cccTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGc                                 >  1:91592/1‑110 (MQ=255)                                                                                                                   cccTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGc                                 <  2:91592/110‑1 (MQ=255)                                                                                                                    ccTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATAccc   <  1:549589/139‑1 (MQ=255)                                                                                                                     cTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATAcccc  <  2:657013/139‑1 (MQ=255)                                                                                                                                      aaTCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                <  2:789721/44‑1 (MQ=38)                                                                                                                                         |                                                                                                                       AAAAACATCAGCGATGACGACCACACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCC  >  NC_000913/2494814‑2495068

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