Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,288,561	(C)7→6	pseudogene (354/2525 nt)	yejO ←	pseudogene, autotransporter outer membrane homology,putative transport, Not classified, putative ATP‑binding component of a transport system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,288,555	0	C	.	100.0%	276.4 / NA	68	pseudogene (360/2525 nt)	yejO	pseudogene, autotransporter outer membrane homology,putative transport, Not classified, putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (36/32);  total (36/32)

GGTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCG  >  NC_000913/2288423‑2288688                                                                                                                                     |                                                                                                                                      ggTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGcccccca                                                                                                                                 <  1:501419/139‑2 (MQ=255)           gAATACCTCCCTGGCCCGCGTTGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCCAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCaa                                                                                                                      <  1:364720/139‑1 (MQ=255)                      tGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           >  1:1475999/1‑139 (MQ=255)                           cgcgTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgcccgcc                                                                                                      >  2:838144/1‑139 (MQ=255)                              gTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTACAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCttt                                                                                                   >  2:989183/1‑139 (MQ=255)                                           tAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCa                                                                                      >  2:475142/1‑139 (MQ=255)                                              cACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTc                                                                                   >  1:1666658/1‑139 (MQ=255)                                               aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCtt                                                                                                    >  1:783800/1‑121 (MQ=255)                                               aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCtt                                                                                                    <  2:783800/121‑1 (MQ=255)                                                ccGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           <  1:115969/113‑1 (MQ=255)                                                ccGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           >  2:115969/1‑113 (MQ=255)                                                 cGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCt                                                                                >  1:668699/1‑139 (MQ=255)                                                  gCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  1:475142/139‑1 (MQ=255)                                                   ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTgc                                                                              <  1:1339246/139‑1 (MQ=255)                                                   ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTgc                                                                              <  1:838144/139‑1 (MQ=255)                                                    cTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCtt                                                                                                                            >  1:1286633/1‑92 (MQ=255)                                                    cTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCtt                                                                                                                            <  2:1286633/92‑1 (MQ=255)                                                         ttCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgccc                                                                                                         <  1:655381/106‑1 (MQ=255)                                                         ttCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgccc                                                                                                         >  2:655381/1‑106 (MQ=255)                                                                 tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATGGCGACAAAGACGG‑CCCCCTATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAc                                                                     <  2:1181339/134‑1 (MQ=255)                                                                 tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  2:1233257/124‑1 (MQ=255)                                                                 tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  1:1233257/1‑124 (MQ=255)                                                                 tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATaacc                                                                     >  1:1181339/1‑132 (MQ=255)                                                                 tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTa                                                                >  2:182058/1‑139 (MQ=255)                                                                      caaatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcaca                                                           <  2:816628/139‑1 (MQ=255)                                                                       aaatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCacac                                                          <  1:182058/139‑1 (MQ=255)                                                                        aatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACAcc                                                         >  1:1226513/1‑139 (MQ=255)                                                                         atgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           <  1:1111935/88‑1 (MQ=255)                                                                         atgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAgcgc                                                                                                           >  2:1111935/1‑88 (MQ=255)                                                                          tgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTaa                                                                                                 >  1:1336650/1‑97 (MQ=255)                                                                          tgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTaa                                                                                                 <  2:1336650/97‑1 (MQ=255)                                                                              aatCCCCCCTCTTTAACCTCCAGTAGTCATCCCTGTTCGATAGCGACAAGGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCattatt                                                   >  2:1095298/1‑139 (MQ=255)                                                                                   ccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCg                                              <  1:989183/139‑1 (MQ=255)                                                                                    cccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGt                                             >  1:1635638/1‑139 (MQ=255)                                                                                           tAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGc                                               <  1:1500231/130‑1 (MQ=255)                                                                                           tAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGc                                               >  2:1500231/1‑130 (MQ=255)                                                                                           tAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTgg                                      >  1:1438655/1‑139 (MQ=255)                                                                                                 ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGcc                                >  2:1473227/1‑139 (MQ=255)                                                                                                  cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCa                                                                          >  2:437702/1‑96 (MQ=255)                                                                                                  cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCa                                                                          <  1:437702/96‑1 (MQ=255)                                                                                                   agtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCttt                                                                                                                           >  2:10325/1‑46 (MQ=38)                                                                                                   agtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCttt                                                                                                                           <  1:10325/46‑1 (MQ=38)                                                                                                   agtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGa                                                 >  1:1198927/1‑120 (MQ=255)                                                                                                   agtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCcgc                              >  1:330159/1‑139 (MQ=255)                                                                                                   agtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGCGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGa                                                 <  2:1198927/120‑1 (MQ=255)                                                                                                    gtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCcgcg                             <  2:668699/139‑1 (MQ=255)                                                                                                        tcctccCTTTTCGATAGCGACAAAGACGG‑CTCCCCATCTTTTTCa                                                                                                                       <  1:1422387/45‑1 (MQ=37)                                                                                                        tcctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCa                                                                                                                       >  2:1422387/1‑45 (MQ=38)                                                                                                        tcctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATc                                                                                               >  2:1460316/1‑69 (MQ=255)                                                                                                        tcctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATc                                                                                               <  1:1460316/69‑1 (MQ=255)                                                                                                        tcctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATa                                                                        <  1:537759/92‑1 (MQ=255)                                                                                                        tcctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATa                                                                        >  2:537759/1‑92 (MQ=255)                                                                                                        tcctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                             >  1:300731/1‑103 (MQ=255)                                                                                                        tcctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                                                             <  2:300731/103‑1 (MQ=255)                                                                                                         cctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcaca                                                           <  1:352803/104‑1 (MQ=255)                                                                                                         cctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcaca                                                           >  2:352803/1‑104 (MQ=255)                                                                                                          ctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTg                                  <  2:1587711/128‑1 (MQ=255)                                                                                                          ctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTg                                  >  1:1587711/1‑128 (MQ=255)                                                                                                          ctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCAccc                       >  1:1070938/1‑139 (MQ=255)                                                                                                          ctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCAccc                       >  2:1510984/1‑139 (MQ=255)                                                                                                            cccTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgcc                                                                                                          <  2:1177095/54‑1 (MQ=39)                                                                                                            cccTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgcc                                                                                                          >  1:1177095/1‑54 (MQ=255)                                                                                                            cccTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTg                                                                  >  2:571036/1‑94 (MQ=255)                                                                                                            cccTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTg                                                                  <  1:571036/94‑1 (MQ=255)                                                                                                               ttttCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               <  2:653080/94‑1 (MQ=255)                                                                                                               ttttCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               >  1:653080/1‑94 (MQ=255)                                                                                                                     aTAGCGACAAATACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTAGTTTGCCGCGCCACCCATTTGAACAAt            <  1:1095298/139‑1 (MQ=255)                                                                                                                          gACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGAt       <  2:1635638/139‑1 (MQ=255)                                                                                                                               aGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCg  <  2:1438655/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                      GGTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCG  >  NC_000913/2288423‑2288688

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