Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,331,613	A→G	pseudogene (2790/4605 nt)	yfaS ←	putative macroglobulin family protein, N‑terminal fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,331,613	0	A	G	100.0%	206.1 / NA	56	pseudogene (2790/4605 nt)	yfaS	putative macroglobulin family protein, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/25);  total (31/25)

AGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACGTGAAATACGCT  >  NC_000913/2331485‑2331748                                                                                                                                 |                                                                                                                                        agcgtcagcgtcTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCAt                                                                                                                              >  2:117750/1‑140 (MQ=255)     tcagcgtcTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTag                                                                                                                          >  2:953033/1‑140 (MQ=255)     tcagcgtcTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTag                                                                                                                          >  1:121713/1‑140 (MQ=255)        gcgtcTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAgaga                                                                                                                       >  1:20378/1‑140 (MQ=255)                 gcATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTcc                                                                                                              <  2:749296/140‑1 (MQ=255)                    tctcAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGaa                                                                                                           <  1:163560/140‑1 (MQ=255)                     ctcAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAAc                                                                                                          >  1:706019/1‑140 (MQ=255)                             cTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGCAGAGATTTTTTTCCGAACGCAGATaa                                                                                                  >  2:1053404/1‑140 (MQ=255)                                    tttAGTCACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCACGTAGAGATTTTTTTCCGAACGCAGATAAGCACGcc                                                                                           <  1:83905/140‑1 (MQ=255)                                        gTCACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTg                                                                                       <  1:783126/140‑1 (MQ=255)                                         tCACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGc                                                                                      >  2:118485/1‑140 (MQ=255)                                          cACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTg                                                                                       <  2:217749/138‑1 (MQ=255)                                          cACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTg                                                                                       >  1:217749/1‑138 (MQ=255)                                                 gcCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGCAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAg                                                                              <  1:1053404/140‑1 (MQ=255)                                                   cACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGcc                                                                            <  1:117750/140‑1 (MQ=255)                                                   cACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGcc                                                                            <  2:20378/140‑1 (MQ=255)                                                   cACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGcc                                                                            >  1:324696/1‑140 (MQ=255)                                                      cGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACTCCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGt                                                                         <  1:556981/140‑1 (MQ=255)                                                        gTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCg                                                                       <  2:706019/140‑1 (MQ=255)                                                          tCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGA‑TTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCg                                                                    <  1:953033/140‑1 (MQ=255)                                                                 ccTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATATACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACACAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATc                                                              >  1:763411/1‑140 (MQ=255)                                                                   tgctgACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATccc                                                            >  1:1086783/1‑140 (MQ=255)                                                                   tgctgACTGAAGAT‑AACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCa                                                           >  1:1077002/1‑140 (MQ=255)                                                                     ctgactgaAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCAc                                                          >  2:533648/1‑140 (MQ=255)                                                                      tgactgaAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCAcg                                                         >  1:112327/1‑140 (MQ=255)                                                                       gactgaAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCAcgc                                                        >  1:103987/1‑140 (MQ=255)                                                                        actgaAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATATACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACACAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCAcgcg                                                       <  2:763411/140‑1 (MQ=255)                                                                        actgaAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCAcgcg                                                       <  2:643193/140‑1 (MQ=255)                                                                         ctgaAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCAcgcgc                                                      <  1:118485/140‑1 (MQ=255)                                                                                     aGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCa                                          <  2:112327/140‑1 (MQ=255)                                                                                      gTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGAt                                                >  1:403975/1‑133 (MQ=255)                                                                                      gTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGAt                                                <  2:403975/133‑1 (MQ=255)                                                                                       tCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGc                                        >  2:922043/1‑140 (MQ=255)                                                                                       tCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGc                                        >  2:369464/1‑140 (MQ=255)                                                                                       tCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGc                                        >  2:1157833/1‑140 (MQ=255)                                                                                                 ggTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGc                                        <  2:458496/130‑1 (MQ=255)                                                                                                 ggTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGc                                        >  1:458496/1‑130 (MQ=255)                                                                                                 ggTTTGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGa                              >  2:16072/1‑140 (MQ=255)                                                                                                   tttGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAAt                            <  1:772632/140‑1 (MQ=255)                                                                                                   tttGTCGCCCACGCGATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAAt                            <  1:533648/140‑1 (MQ=255)                                                                                                              cgcgATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAgg                 >  1:12703/1‑140 (MQ=255)                                                                                                                 gATACACCGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAAc              <  2:103987/140‑1 (MQ=255)                                                                                                                        cGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACGTGAAAt       >  2:131091/1‑140 (MQ=255)                                                                                                                        cGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACGTGAAAt       <  1:207409/140‑1 (MQ=255)                                                                                                                        cGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACGTGAAAt       >  2:125149/1‑140 (MQ=255)                                                                                                                        cGTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACGTGAAAt       >  1:1062087/1‑140 (MQ=255)                                                                                                                         gTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACGTGAAATa      >  2:829871/1‑140 (MQ=255)                                                                                                                         gTTGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACGTGAAATa      >  1:695294/1‑140 (MQ=255)                                                                                                                           tGGCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACGTGAAATACg    <  1:1157833/140‑1 (MQ=255)                                                                                                                            ggCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATc                    <  2:417730/123‑1 (MQ=255)                                                                                                                            ggCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATc                    >  1:417730/1‑123 (MQ=255)                                                                                                                             gCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACg             <  2:126852/129‑1 (MQ=255)                                                                                                                             gCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACg             >  1:126852/1‑129 (MQ=255)                                                                                                                             gCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACGTGAAATACGCt  <  1:16072/140‑1 (MQ=255)                                                                                                                             gCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACGTGAAATACGCt  >  1:835026/1‑140 (MQ=255)                                                                                                                             gCATGCTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACGTGAAATACGCt  >  1:621594/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                                        AGCGTCAGCGTCTGGCGCATCTCAGCGCCTGCAAATTTAGTCACCAGCGCCACCGGTTCGTTATCCTGCTGACTGAAGATAAACAGTCCTGCCGCCGGTTTGTCGCCCACGCGATACACCGTTGGCATACTCCACTTCATGTAGAGATTTTTTTCCGAACGCAGATAAGCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACGTGAAATACGCT  >  NC_000913/2331485‑2331748

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