Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,034,536	G→A	pseudogene (486/663 nt)	yedS →	pseudogene, outer membrane protein homology, putative outer membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,034,536	0	G	A	100.0%	220.4 / NA	67	pseudogene (486/663 nt)	yedS	pseudogene, outer membrane protein homology, putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/33);  total (34/33)

GTCGTATATGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAGGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTATGTCAA  >  NC_000913/2034402‑2034669                                                                                                                                       |                                                                                                                                      gtCGTATATGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                                   <  2:35010/139‑1 (MQ=255) gtCGTATATGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                                   <  1:169174/139‑1 (MQ=255)   cGTATATGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGAt                                                                                                                                 <  2:523179/139‑1 (MQ=255)   cGTATATGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGAt                                                                                                                                 <  1:36251/139‑1 (MQ=255)     tataTGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATgg                                                                                                                               >  2:1207976/1‑139 (MQ=255)      atatGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGc                                                                                                                              <  1:1602461/139‑1 (MQ=255)            cGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGc                                                                                                                              <  2:1530574/133‑1 (MQ=255)            cGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGc                                                                                                                              >  1:1530574/1‑133 (MQ=255)            cGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAAc                                                                                                                        >  1:707159/1‑139 (MQ=255)            cGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAAc                                                                                                                        >  2:1424467/1‑139 (MQ=255)            cGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAAc                                                                                                                        >  2:113527/1‑139 (MQ=255)               aGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActtt                                                                                                                     >  2:127236/1‑139 (MQ=255)                  aGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgct                                                                                                                  <  1:1207976/139‑1 (MQ=255)                      tGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGaa                                                                                                                         >  1:11816/1‑128 (MQ=255)                      tGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGaa                                                                                                                         <  2:11816/128‑1 (MQ=255)                      tGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACtttgctttgc                                                                                                              <  1:1424467/139‑1 (MQ=255)                      tGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACtttgctttgc                                                                                                              <  2:1546765/139‑1 (MQ=255)                      tGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACtttgctttgc                                                                                                              >  2:384164/1‑139 (MQ=255)                      tGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACtttgctttgc                                                                                                              >  2:688774/1‑139 (MQ=255)                       ggACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCa                                                                                                             <  1:113527/139‑1 (MQ=255)                                               gtggCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctttg                                                                                                               <  2:226621/113‑1 (MQ=255)                                               gtggCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctttg                                                                                                               >  1:226621/1‑113 (MQ=255)                                                            tGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                                   <  2:748371/80‑1 (MQ=255)                                                            tGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                                   >  1:748371/1‑80 (MQ=255)                                                            tGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCggg                                                                        <  1:688774/139‑1 (MQ=255)                                                            tGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCggg                                                                        >  1:872400/1‑139 (MQ=255)                                                             gCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGc                                                                       >  1:1418406/1‑139 (MQ=255)                                                                ggCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAgg                                                                    <  2:707159/139‑1 (MQ=255)                                                                 gcgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                                   <  1:65499/75‑1 (MQ=255)                                                                 gcgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                                   >  2:65499/1‑75 (MQ=255)                                                                 gcgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGGActttgctttg                                                                                                               >  1:899051/1‑95 (MQ=255)                                                                 gcgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctttg                                                                                                               <  2:899051/95‑1 (MQ=255)                                                                   gcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAg                                                                                                      >  1:1564265/1‑102 (MQ=255)                                                                   gcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAg                                                                                                      <  2:1564265/102‑1 (MQ=255)                                                                      gACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGa                                                                                           >  1:788575/1‑110 (MQ=255)                                                                      gACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGa                                                                                           <  2:788575/110‑1 (MQ=255)                                                                      gACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGtt                                                              >  1:1286935/1‑139 (MQ=255)                                                                       aCAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                                   <  1:622880/69‑1 (MQ=255)                                                                       aCAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                                   >  2:622880/1‑69 (MQ=255)                                                                       aCAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTc                                                             <  2:872400/139‑1 (MQ=255)                                                                                   aCGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTCCGGTAGTGGTGa                                                 <  1:127236/139‑1 (MQ=255)                                                                                       tCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGc                                                                                                                              >  1:757812/1‑58 (MQ=255)                                                                                       tCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGc                                                                                                                              <  2:757812/58‑1 (MQ=255)                                                                                         gTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCg                                           >  2:1625080/1‑139 (MQ=255)                                                                                          tGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGt                                                               <  1:1101951/118‑1 (MQ=255)                                                                                          tGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGt                                                               >  2:1101951/1‑118 (MQ=255)                                                                                                 ggTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCttt                                   >  1:30710/1‑139 (MQ=255)                                                                                                 ggTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCttt                                   >  1:673996/1‑139 (MQ=255)                                                                                                 ggTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCCTTGCAATATCAGGGTAATAATg                                                                                            >  1:295500/1‑82 (MQ=255)                                                                                                 ggTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCCTTGCAATATCAGGGTAATAATg                                                                                            <  2:295500/82‑1 (MQ=255)                                                                                                     ttGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                                   <  2:1370900/39‑1 (MQ=255)                                                                                                     ttGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                                   >  1:1370900/1‑39 (MQ=255)                                                                                                           ccTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTcca                                                           >  1:1084382/1‑104 (MQ=255)                                                                                                           ccTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTCCg                                                           <  2:1084382/105‑1 (MQ=255)                                                                                                               ccGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctt                                                                                                                 <  2:526691/47‑1 (MQ=255)                                                                                                               ccGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctt                                                                                                                 >  1:526691/1‑47 (MQ=255)                                                                                                                cGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTg                    <  1:921807/139‑1 (MQ=255)                                                                                                                    aCAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATgg                >  1:967374/1‑139 (MQ=255)                                                                                                                     cAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGt               >  2:1070951/1‑139 (MQ=255)                                                                                                                       aTGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGttt             >  2:826388/1‑139 (MQ=255)                                                                                                                         ggCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCg           >  1:1611000/1‑139 (MQ=255)                                                                                                                         ggCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCg           >  1:1225619/1‑139 (MQ=255)                                                                                                                          gCTTCTTTGGTTAAGTCGATGGCTTGAACtttgctttgc                                                                                                              >  1:814878/1‑39 (MQ=255)                                                                                                                          gCTTCTTTGGTTAAGTCGATGGCTTGAACtttgctttgc                                                                                                              <  2:814878/39‑1 (MQ=255)                                                                                                                           cttcttTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGt         <  2:1225619/139‑1 (MQ=255)                                                                                                                            ttcttTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTa        <  2:1418406/139‑1 (MQ=255)                                                                                                                                  ggTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGATCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTATGTCaa  <  1:1070951/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      GTCGTATATGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAGGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTATGTCAA  >  NC_000913/2034402‑2034669

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