Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,187,588	G→A	pseudogene (162/2439 nt)	yqiG →	pseudogene, fimbrial export usher family,putative membrane, Not classified, putative membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,187,588	0	G	A	100.0%	214.1 / NA	66	pseudogene (162/2439 nt)	yqiG	pseudogene, fimbrial export usher family,putative membrane, Not classified, putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/32);  total (34/32)

ATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGATGGTTGT  >  NC_000913/3187455‑3187723                                                                                                                                      |                                                                                                                                        aTGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAAt                                                                                                                                    <  2:1072277/139‑1 (MQ=255) aTGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAAt                                                                                                                                    >  2:1540400/1‑139 (MQ=255)  tGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATc                                                                                                                                   >  1:1312782/1‑139 (MQ=255)       ttGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAGTTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCt                                                                                                                              >  2:1188842/1‑139 (MQ=255)        tGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCtt                                                                                                                             >  2:327033/1‑139 (MQ=255)                acaaAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCa                                                                                                                     <  1:1254706/139‑1 (MQ=255)                   aaGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATg                                                                                                                  <  2:1606553/139‑1 (MQ=255)                      aTTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                                    <  2:247151/134‑1 (MQ=255)                      aTTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                                    >  1:247151/1‑134 (MQ=255)                         tAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTAtt                                                                                                            <  2:536218/139‑1 (MQ=255)                                aGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTACGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGAt                                                                                                     >  2:1188419/1‑139 (MQ=255)                                aGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGAt                                                                                                     >  1:1548390/1‑139 (MQ=255)                                    gAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAAcc                                                                                                 <  1:327033/139‑1 (MQ=255)                                            tGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAAt                                                                                                                   <  2:1245672/113‑1 (MQ=255)                                            tGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAAt                                                                                                                   >  1:1100498/1‑113 (MQ=255)                                            tGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAAt                                                                                                                   <  2:1100498/113‑1 (MQ=255)                                            tGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAAt                                                                                                                   >  1:1245672/1‑113 (MQ=255)                                                  cGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGc                                                                                   >  2:1139500/1‑139 (MQ=255)                                                   gTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCAGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCg                                                                                  <  2:1344249/139‑1 (MQ=255)                                                    ttAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGt                                                                                 <  1:1592093/139‑1 (MQ=255)                                                      aaCCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                                    >  1:352048/1‑102 (MQ=255)                                                      aaCCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                                    <  2:352048/102‑1 (MQ=255)                                                      aaCCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGa                                                                               <  1:1188419/139‑1 (MQ=255)                                                         cTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGacac                                                                            <  2:703059/139‑1 (MQ=255)                                                         cTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGacac                                                                            <  1:1139500/139‑1 (MQ=255)                                                         cTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGacac                                                                            <  1:388985/139‑1 (MQ=255)                                                            tCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCag                                                                         >  1:391416/1‑139 (MQ=255)                                                             cccAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGa                                                                                                                           >  2:380337/1‑88 (MQ=255)                                                             cccAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGa                                                                                                                           <  1:380337/88‑1 (MQ=255)                                                             cccAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCaga                                                                        >  1:1231126/1‑139 (MQ=255)                                                             cccAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCaga                                                                        <  2:1382357/139‑1 (MQ=255)                                                             cccAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCag                                                                         >  1:1433221/1‑138 (MQ=255)                                                             cccAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCag                                                                         <  2:1433221/138‑1 (MQ=255)                                                              ccAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCagag                                                                       <  2:1520846/139‑1 (MQ=255)                                                                  tttGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGtt                                                                                      <  2:216859/120‑1 (MQ=255)                                                                  tttGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGtt                                                                                      >  1:216859/1‑120 (MQ=255)                                                                     gAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATg                                                                                                                  >  1:1513042/1‑89 (MQ=255)                                                                     gAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATg                                                                                                                  <  2:1513042/89‑1 (MQ=255)                                                                       aaCTGATAGCAAATTACCCGTCGGCAAATATTCACTAAGCTCTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACa                                                              <  2:430451/139‑1 (MQ=255)                                                                          tGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACatta                                                           >  1:176248/1‑139 (MQ=255)                                                                            aTGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACattatt                                                         >  1:823979/1‑139 (MQ=255)                                                                                  aaTTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGAtt                                                   <  2:1312782/139‑1 (MQ=255)                                                                                             ggCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAAt                                        >  2:255765/1‑139 (MQ=255)                                                                                                  atatTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGa                                                                               >  2:1321834/1‑95 (MQ=255)                                                                                                  atatTCACTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGa                                                                               <  1:1321834/95‑1 (MQ=255)                                                                                                        aCTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCa                                                                     >  1:178884/1‑99 (MQ=255)                                                                                                        aCTAAGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCa                                                                     <  2:178884/99‑1 (MQ=255)                                                                                                            aGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                                    >  1:790625/1‑48 (MQ=255)                                                                                                            aGCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCaa                                                                                                                    <  2:790625/48‑1 (MQ=255)                                                                                                             gCACTCTGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCa                        >  1:1554069/1‑139 (MQ=255)                                                                                                                 tctGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGacac                                                                            >  2:446616/1‑83 (MQ=255)                                                                                                                 tctGATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGacac                                                                            <  1:446616/83‑1 (MQ=255)                                                                                                                    gATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACattatt                                                         <  1:976702/99‑1 (MQ=255)                                                                                                                    gATTAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACattatt                                                         >  2:976702/1‑99 (MQ=255)                                                                                                                      ttAATAATAAGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCt               >  2:1593805/1‑139 (MQ=255)                                                                                                                        aataataaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGa             >  1:48129/1‑139 (MQ=255)                                                                                                                           aataaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCg          >  2:464605/1‑139 (MQ=255)                                                                                                                           aataaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCg          >  1:375269/1‑139 (MQ=255)                                                                                                                            ataaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGa         <  2:1548390/139‑1 (MQ=255)                                                                                                                             taaGAGGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGAt        <  1:1593805/139‑1 (MQ=255)                                                                                                                               agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAAc                                                               >  1:1554995/1‑82 (MQ=255)                                                                                                                               agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAAc                                                               <  2:1554995/82‑1 (MQ=255)                                                                                                                               agagGACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGATgg      >  2:495667/1‑139 (MQ=255)                                                                                                                                   gACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGtt                                                                                      <  2:700062/55‑1 (MQ=255)                                                                                                                                   gACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGtt                                                                                      >  1:700062/1‑55 (MQ=255)                                                                                                                                   gACACCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGATGgttgt  >  1:53646/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                        ATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGGCAAATATTCACTAAGCACTCTGATTAATAATAAGAGGACGCCAATCCACCTTGACCTCCAATGGGTATTAATTGATAACCAAACTGCAGTTTGCGTGACACCAGAGCAATTAACATTATTAGGATTTACTGATGAATTTATTGAAAAAACTCAGCAAAACCTGATCGATGGTTGT  >  NC_000913/3187455‑3187723

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: ‑