Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,465,668	G→C	pseudogene (277/2518 nt)	ydbA →	putative outer membrane protein, N‑terminal fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,465,668	0	G	C	100.0%	206.4 / NA	53	pseudogene (277/2518 nt)	ydbA	putative outer membrane protein, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (26/27);  total (26/27)

CAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTC  >  NC_000913/1465535‑1465800                                                                                                                                      |                                                                                                                                     cAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGatcatcatc                                                                                                                                <  1:563179/140‑1 (MQ=255)        gAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCg                                                                                                                         <  1:247544/140‑1 (MQ=255)             gCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGcaaca                                                                                                                    <  2:587622/140‑1 (MQ=255)              cTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGcaacaa                                                                                                                   <  2:864337/140‑1 (MQ=255)              cTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGcaacaa                                                                                                                   >  2:131556/1‑140 (MQ=255)                     gCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCa                                                                                                            <  1:642327/140‑1 (MQ=255)                     gCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCa                                                                                                            <  1:518131/140‑1 (MQ=255)                       ttACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACt                                                                                                          >  2:654452/1‑140 (MQ=255)                        tACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTc                                                                                                         >  2:873678/1‑140 (MQ=255)                        tACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTc                                                                                                         >  1:577589/1‑140 (MQ=255)                               tGAATGCTCTGCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGa                                                                                                  >  2:601887/1‑140 (MQ=255)                               tGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGa                                                                                                  >  1:647669/1‑140 (MQ=255)                               tGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGa                                                                                                  >  1:524708/1‑140 (MQ=255)                                gAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGac                                                                                                 <  1:607911/140‑1 (MQ=255)                                 aaTGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGaca                                                                                                <  1:491414/140‑1 (MQ=255)                                 aaTGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGaca                                                                                                <  1:607534/140‑1 (MQ=255)                                 aaTGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGaca                                                                                                <  1:131556/140‑1 (MQ=255)                                    gCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGaca                                                                                                <  1:783460/137‑1 (MQ=255)                                    gCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGaca                                                                                                >  2:783460/1‑137 (MQ=255)                                     ctctCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACAcccc                                                                                            <  2:684889/140‑1 (MQ=255)                                     ctctCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACAcccc                                                                                            >  1:684889/1‑140 (MQ=255)                                     ctctCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACAcccc                                                                                            <  1:17120/140‑1 (MQ=255)                                     ctctCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACAcccc                                                                                            <  1:530471/140‑1 (MQ=255)                                      tctcCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACAccccc                                                                                           >  2:745352/1‑140 (MQ=255)                                         ccTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtga                                                                                        >  2:850878/1‑140 (MQ=255)                                              aGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCtgatgatg                                                                                   <  1:873678/140‑1 (MQ=255)                                                tCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGaa                                                                                 >  2:372911/1‑140 (MQ=255)                                                    aCAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAg                                                                             >  2:83395/1‑140 (MQ=255)                                                      aaaaCCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAgac                                                                           <  2:577589/140‑1 (MQ=255)                                                        aaCCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAgacga                                                                         <  1:601887/140‑1 (MQ=255)                                                                 ccATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGt                                                                <  1:850878/140‑1 (MQ=255)                                                                    tGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAga                                                                            <  2:66189/125‑1 (MQ=255)                                                                    tGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAga                                                                            >  1:66189/1‑125 (MQ=255)                                                                    tGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTcc                                                             >  1:219876/1‑140 (MQ=255)                                                                    tGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTcc                                                             >  2:450819/1‑140 (MQ=255)                                                                                ggCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTgg                                                 <  2:768442/140‑1 (MQ=255)                                                                                ggCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTgg                                                 >  2:531455/1‑140 (MQ=255)                                                                                    ctgctgCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGga                                             >  1:666677/1‑140 (MQ=255)                                                                                    ctgctgCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGga                                             >  1:751977/1‑140 (MQ=255)                                                                                      gctgcGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACTCTCCAGTTCCCCCAACTCCTGGCGgaga                                           <  1:401782/140‑1 (MQ=255)                                                                                               aTCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGataat                                  >  2:490591/1‑140 (MQ=255)                                                                                                    ggCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCgg                             >  2:859430/1‑140 (MQ=255)                                                                                                    ggCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCgg                             >  1:196110/1‑140 (MQ=255)                                                                                                     gCCTTAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGa                            <  1:450819/140‑1 (MQ=255)                                                                                                         tAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGAt                        <  2:751977/140‑1 (MQ=255)                                                                                                         tAGCCGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGAt                        <  1:531455/140‑1 (MQ=255)                                                                                                            ccGTAGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCg                     <  2:647669/140‑1 (MQ=255)                                                                                                                aGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGgatg                 >  2:796417/1‑140 (MQ=255)                                                                                                                aGTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGgatg                 >  1:468826/1‑140 (MQ=255)                                                                                                                 gTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGgatga                <  1:372911/140‑1 (MQ=255)                                                                                                                 gTGGAACTAAACGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGgatga                <  2:524708/140‑1 (MQ=255)                                                                                                                          aaCGACGATGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTac       <  1:859430/140‑1 (MQ=255)                                                                                                                            cgacgaTGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTacac     >  1:466082/1‑140 (MQ=255)                                                                                                                               cgaTGATCATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTc  >  2:755597/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                     CAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCACGGCCTTAGCCGTAGTGGAACTAAACGACGATGATGATCATCATCATCGCAACAATTCTCCACTCCCACCGACACCCCCTGATGATGAATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGATCCGGATGATACGCCTACACCTC  >  NC_000913/1465535‑1465800

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