Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,453,265	C→T	pseudogene (1737/2646 nt)	yfcU ←	pseudogene, FimD fimbrial export usher family,putative membrane, Not classified, putative outer membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,453,265	0	C	T	100.0%	239.1 / NA	73	pseudogene (1737/2646 nt)	yfcU	pseudogene, FimD fimbrial export usher family,putative membrane, Not classified, putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/38);  total (35/38)

CCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAT  >  NC_000913/2453131‑2453395                                                                                                                                       |                                                                                                                                   ccGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACa                                                                                                                                >  2:433987/1‑139 (MQ=255)         cAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTa                                                                                                                        >  1:1027697/1‑139 (MQ=255)               aTAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGcc                                                                                                                  <  1:359570/139‑1 (MQ=255)                 aGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCTATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                 <  2:973522/122‑1 (MQ=255)                 aGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCAAAATTTAc                                                                                                                                 >  1:973522/1‑122 (MQ=255)                     agacgcgTCATAGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCAGTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                               >  2:873389/4‑136 (MQ=255)                     aGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                               <  1:873389/136‑1 (MQ=255)                       tCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTc                                                                                                          >  1:2458/1‑139 (MQ=255)                         gcgTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACa                                                                                                                                <  1:1417877/115‑1 (MQ=255)                         gcgTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACa                                                                                                                                >  2:1417877/1‑115 (MQ=255)                         gcgTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCAt                                                                                                        >  1:1036729/1‑139 (MQ=255)                          cgTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATa                                                                                                       >  2:467136/1‑139 (MQ=255)                            tCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACt                                                                                                     <  2:2458/139‑1 (MQ=255)                              aTCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCg                                                                                                   >  1:308784/1‑139 (MQ=255)                                        cTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                 >  1:131954/1‑99 (MQ=255)                                        cTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                 <  2:131954/99‑1 (MQ=255)                                                aCCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGt                                                                                     >  1:1104453/1‑135 (MQ=255)                                                aCCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGt                                                                                     <  2:1104453/135‑1 (MQ=255)                                                 ccGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAg                                                                                                >  2:45987/1‑123 (MQ=255)                                                 ccGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAg                                                                                                <  1:45987/123‑1 (MQ=255)                                                  cGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                 <  2:656346/89‑1 (MQ=255)                                                  cGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                 >  1:656346/1‑89 (MQ=255)                                                        ggCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                         <  2:1027697/139‑1 (MQ=255)                                                        ggCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACCGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                         <  2:1036729/139‑1 (MQ=255)                                                                   gtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGaa                                                              >  1:712952/1‑139 (MQ=255)                                                                     cccACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGGATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACCCTCATATTACGAAtg                                                            <  2:1006085/139‑1 (MQ=255)                                                                          tGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                               <  2:53504/83‑1 (MQ=255)                                                                          tGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCgg                                                                                                               >  1:53504/1‑83 (MQ=255)                                                                          tGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                       <  2:308784/139‑1 (MQ=255)                                                                               tAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCatat                                                  <  1:467136/139‑1 (MQ=255)                                                                                      cGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTc                                                                        >  2:964068/1‑110 (MQ=255)                                                                                      cGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTc                                                                        <  1:964068/110‑1 (MQ=255)                                                                                            aGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGga                                     <  1:834282/139‑1 (MQ=255)                                                                                                aaCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCg                                                                                                                >  1:715119/1‑60 (MQ=255)                                                                                                aaCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCg                                                                                                                <  2:715119/60‑1 (MQ=255)                                                                                                aaCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTa                                                                                                 <  2:160502/75‑1 (MQ=255)                                                                                                aaCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTa                                                                                                 >  1:160502/1‑75 (MQ=255)                                                                                                aaCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGaaa                                                                             >  2:394978/1‑95 (MQ=255)                                                                                                aaCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGaaa                                                                             <  1:394978/95‑1 (MQ=255)                                                                                                  cGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCat                               <  1:433987/139‑1 (MQ=255)                                                                                                     tACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACt                                                                                                     <  1:529805/66‑1 (MQ=255)                                                                                                     tACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACt                                                                                                     >  2:529805/1‑66 (MQ=255)                                                                                                       cTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAAct                                                                           >  2:1327728/1‑89 (MQ=255)                                                                                                       cTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACg                                                                           <  1:1327728/90‑1 (MQ=255)                                                                                                            gtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGc                                                                                               >  2:936200/1‑65 (MQ=255)                                                                                                            gtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGc                                                                                               <  1:936200/65‑1 (MQ=255)                                                                                                            gtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCatg                              <  2:814714/130‑1 (MQ=255)                                                                                                            gtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCatg                              >  1:814714/1‑130 (MQ=255)                                                                                                             tCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAg                                                                                                >  2:768553/1‑63 (MQ=255)                                                                                                             tCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAg                                                                                                <  1:768553/63‑1 (MQ=255)                                                                                                               gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCg                                                                                                                    <  1:490969/41‑1 (MQ=255)                                                                                                               gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCg                                                                                                                    >  2:490969/1‑41 (MQ=255)                                                                                                                    aCGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtgc                                                           >  2:683764/1‑93 (MQ=255)                                                                                                                    aCGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtgc                                                           <  1:683764/93‑1 (MQ=255)                                                                                                                    aCGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTg             >  1:1400722/1‑139 (MQ=255)                                                                                                                     cGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGc            <  2:90216/139‑1 (MQ=255)                                                                                                                           aCTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGt                                                                                     >  2:1238951/1‑60 (MQ=255)                                                                                                                           aCTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGt                                                                                     <  1:1238951/60‑1 (MQ=255)                                                                                                                              gAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                          >  1:1277113/1‑52 (MQ=255)                                                                                                                              gAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAg                                                                                          <  2:1277113/52‑1 (MQ=255)                                                                                                                              gAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGa                                                                               >  2:283808/1‑63 (MQ=255)                                                                                                                              gAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGa                                                                               <  1:283808/63‑1 (MQ=255)                                                                                                                               aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGaaa                                                                             >  2:533996/1‑64 (MQ=255)                                                                                                                               aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGaaa                                                                             <  1:533996/64‑1 (MQ=255)                                                                                                                               aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtg                                                            >  2:1104605/1‑81 (MQ=255)                                                                                                                               aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtg                                                            <  1:1104605/81‑1 (MQ=255)                                                                                                                               aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTg             >  1:353265/1‑128 (MQ=255)                                                                                                                               aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTg             <  2:353265/128‑1 (MQ=255)                                                                                                                               aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAt  <  2:712952/139‑1 (MQ=255)                                                                                                                                gCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCatg                              <  2:443979/110‑1 (MQ=255)                                                                                                                                gCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCatg                              >  1:443979/1‑110 (MQ=255)                                                                                                                                gCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAt  >  2:410499/1‑138 (MQ=255)                                                                                                                                gCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAt  <  1:410499/138‑1 (MQ=255)                                                                                                                                       |                                                                                                                                   CCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAT  >  NC_000913/2453131‑2453395

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