Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	735,498	C→A	pseudogene (1279/1434 nt)	ybfO →	conserved protein, rhs‑like YbfO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	735,498	0	C	A	100.0%	199.8 / NA	57	pseudogene (1279/1434 nt)	ybfO	conserved protein, rhs‑like YbfO
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (26/31);  total (26/31)

CGACTACACACTGTATTGCAGGAAGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGCCAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAACAGCAGGGATAGATGATTTGGGTGAAATTTTAGCTAAGATGAAACAG  >  NC_000913/735374‑735626                                                                                                                             |                                                                                                                                 cGACTACATACTGTATTGTAGGTAGTGGGTCAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACt                                                                                                                          <  1:246230/133‑1 (MQ=255)        acacTGTATTGCAGGAAGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTac                                                                                                                   >  2:520421/1‑133 (MQ=255)          acTGTATTGCAGGAAGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTacca                                                                                                                 <  2:475070/133‑1 (MQ=255)          acTGTATTCCAGGAAATGGGTAAAACCAAAACTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTacca                                                                                                                 <  1:213451/133‑1 (MQ=255)           cTGTATTGCAGGAAGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTCCCTGCCACAGACAAATTGGACAACAACTCAATTaccac                                                                                                                <  2:588259/133‑1 (MQ=255)           cTGTATTGCAGGAAGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTaccac                                                                                                                <  1:127229/133‑1 (MQ=255)           cTGTATTGCAGGAAGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTaccac                                                                                                                <  1:459202/133‑1 (MQ=255)           cTGTATTGCAGGAAGTGGGTAAAACCAAACCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTaccac                                                                                                                <  1:585342/133‑1 (MQ=255)            tGTATTGCAGGAAGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTaccacc                                                                                                               >  2:589990/1‑133 (MQ=255)              tATTGCAGGAAGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTaccac                                                                                                                <  1:137852/130‑1 (MQ=255)              tATTGCAGGAAGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTaccac                                                                                                                >  2:137852/1‑130 (MQ=255)                ttGCAGGAAGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGa                                                                                                           <  2:523667/133‑1 (MQ=255)                  gCAGGAAGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGata                                                                                                         <  2:494820/133‑1 (MQ=255)                   cAGGAAGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAAATGGACAAGAACCCACTTACCAACAGGatat                                                                                                        >  1:306502/1‑133 (MQ=255)                     gtaaGTGGGCAAAACTAAATCTGACTTGGACGCCACAGGATAATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATAc                                                                                                      <  2:312877/131‑1 (MQ=255)                     ggAAGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATAc                                                                                                      >  2:120885/1‑133 (MQ=255)                       aaGTGGGTAAAACCACATCTGACTAGTATGCCGCAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACtt                                                                                                    <  2:230569/133‑1 (MQ=255)                        aGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGATAACAACTTAATAACGACCAGGATATACTtg                                                                                                   >  1:214005/1‑133 (MQ=255)                        aGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTtg                                                                                                   >  1:332688/1‑133 (MQ=255)                        aGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTtg                                                                                                   >  1:330040/1‑133 (MQ=255)                        aGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACAACCAGGATATACTtg                                                                                                   >  2:205379/1‑133 (MQ=255)                            ggTAAAACTAAATCTAACCTGTCCGCCACAAGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCa                                                                                               <  1:231536/133‑1 (MQ=255)                            ggTAAAACCAAATCTTACTTGTACTCCAGAGGATGATTCCCGGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCa                                                                                               <  2:214005/133‑1 (MQ=255)                               aaaaCCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGaa                                                                                            <  1:97525/133‑1 (MQ=255)                                agaCCAAATGAGACTTGTACGCCTAAGGTTGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAAc                                                                                           <  2:306502/131‑1 (MQ=255)                                aaaCCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAAc                                                                                           <  2:448768/133‑1 (MQ=255)                                 aaCCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAAcc                                                                                          >  1:369235/1‑133 (MQ=255)                                   ccAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCt                                                                                        >  1:487330/1‑133 (MQ=255)                                            aCTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCtatt                                                                                     <  2:81327/127‑1 (MQ=255)                                            aCTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCtatt                                                                                     >  1:81327/1‑127 (MQ=255)                                            aCTTGTACCCCACAGTATCATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTcc                                                                               <  1:185463/133‑1 (MQ=255)                                                  aCGCCACAGGGTGATGCCCAGTATGCAGTGATGGATACAAAAACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACa                                                                         <  2:272544/133‑1 (MQ=255)                                                        caGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGGAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATtaata                                                                   >  2:401558/1‑133 (MQ=255)                                                        caGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATtaata                                                                   >  2:322878/1‑133 (MQ=255)                                                           gatgatTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGTAGAACCCGATTTATTCACAGATATTAATAAAc                                                                >  2:168108/1‑133 (MQ=255)                                                               atTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGa                                                            <  2:369235/133‑1 (MQ=255)                                                                ttCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGAt                                                           >  2:136842/1‑133 (MQ=255)                                                                     aGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGcc                                                      <  1:520421/133‑1 (MQ=255)                                                                         tCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCa                                                  <  2:340217/133‑1 (MQ=255)                                                                            aGGTATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAAca                                               <  2:617516/133‑1 (MQ=255)                                                                              ggatggatACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAAcagc                                             <  2:330040/133‑1 (MQ=255)                                                                                 tggatACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAACAGCAgg                                          >  2:553832/1‑133 (MQ=255)                                                                                   gatACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTGATAAACCGGATAGGGCAGCAAGAGCAGGGa                                        >  1:197077/1‑133 (MQ=255)                                                                                   gatACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAACAGCAGGGa                                        >  2:511415/1‑133 (MQ=255)                                                                                   gatACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATACCCAGACAATAAGAAATCCGATGGGGCAGCAACAGCAGGGa                                        >  2:177886/1‑133 (MQ=255)                                                                                              aGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAACAGCAGGGATAGATGATTTg                             >  2:149008/1‑133 (MQ=255)                                                                                              aGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAACAGCAGGGATAGATGATTTg                             >  1:394387/1‑133 (MQ=255)                                                                                              aGCAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGAGGGGCCAGGAACGGCTGGGATATATGATTTg                             >  1:170316/1‑133 (MQ=255)                                                                                                cAAGTGATTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAACAGCAGGGATAGATGATTTggg                           >  1:71439/1‑133 (MQ=255)                                                                                                      aTTACCTGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAACAGCAGGGATAGATGATTTGGGTGAAAt                     >  1:389679/1‑133 (MQ=255)                                                                                                            tGCCACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAACAGCAGGGATAGATGATTTGGGTGAAATTTTAGc               <  1:322878/133‑1 (MQ=255)                                                                                                              ccACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAACAGCAGGGATAGATGATTTGGGTGAAATTTTAGCTa             <  2:493381/133‑1 (MQ=255)                                                                                                              ccACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAACAGCAGGGATAGATGATTTGGGTGAAATTTTAGCTa             <  1:511415/133‑1 (MQ=255)                                                                                                              ccACAGACAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAACAGCAGGGATAGATGATTTGGGTGAAATTTTAGCTa             <  2:70770/133‑1 (MQ=255)                                                                                                                    aCAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAACAGCAGGGATAGATGATTTGGGTGAAATTTTAGCTAAg           <  2:489363/129‑1 (MQ=255)                                                                                                                    aCAAATTGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAACAGCAGGGATAGATGATTTGGGTGAAATTTTAGCTAAg           >  1:489363/1‑129 (MQ=255)                                                                                                                         ttGGACAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAACAGCAGGGATAGATGATTTGGGTGAAATTTTAGCTAAGATGAAACag  <  2:332688/133‑1 (MQ=255)                                                                                                                             |                                                                                                                                 CGACTACACACTGTATTGCAGGAAGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGCCAACAACTCAATTACCACCAGGATATACTTGTGCAGAACCCTATTTATTCCCAGACATTAATAAACCCGATGGGCCAGCAACAGCAGGGATAGATGATTTGGGTGAAATTTTAGCTAAGATGAAACAG  >  NC_000913/735374‑735626

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