Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	734,845	A→G	pseudogene (626/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	734,845	0	A	G	100.0%	77.8 / NA	25	pseudogene (626/1434 nt)	ybfO	conserved protein, rhs‑like YbfO
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/14);  total (11/14)

GAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTG  >  NC_000913/734786‑734981                                                            |                                                                                                                                         gagtgagGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTg                                                          <  2:170358/140‑1 (MQ=18)    tgagGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCa                                                       >  2:535771/1‑140 (MQ=21)    tgagGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCa                                                       >  1:866481/1‑140 (MQ=21)        gAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCAc                                                   >  2:641780/1‑140 (MQ=21)         aaaGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACg                                                  >  1:390713/1‑140 (MQ=21)         aaaGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACg                                                  >  2:590783/1‑140 (MQ=21)            gccgccGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGaa                                               <  1:502535/140‑1 (MQ=255)             ccgccgCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAg                                              >  1:443754/1‑140 (MQ=21)              cgccgcTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAgg                                             <  1:628710/140‑1 (MQ=25)                ccgcTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAGAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAgggg                                           >  1:361716/1‑140 (MQ=255)                ccgcTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAgggg                                           <  1:191492/140‑1 (MQ=25)                  gctggctgGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCa                                         >  2:5468/1‑140 (MQ=255)                       ctggcATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGa                                    <  1:413272/140‑1 (MQ=255)                        tggcATCGTGCGGCCTGACGGTGGCGGAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGaa                                   <  1:57235/140‑1 (MQ=255)                         ggcATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAAt                                  <  2:669793/140‑1 (MQ=255)                               gTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcg                            <  1:640013/140‑1 (MQ=255)                                tGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcgc                           >  2:724461/1‑140 (MQ=255)                                tGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcgc                           <  1:641780/140‑1 (MQ=255)                                  cGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAg                         >  1:915383/1‑140 (MQ=255)                                  cGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAg                         >  1:915382/1‑140 (MQ=255)                                        gACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACg                   <  1:535771/140‑1 (MQ=255)                                        gACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACg                   <  1:110260/140‑1 (MQ=255)                                                 gCAGATGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATggg          <  1:535776/140‑1 (MQ=255)                                                      tGCAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAGAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAAc     <  2:361716/140‑1 (MQ=255)                                                        cAAAGCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACct   <  1:273791/140‑1 (MQ=255)                                                         aaagCCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACctg  <  2:111617/140‑1 (MQ=255)                                                            |                                                                                                                                         GAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTG  >  NC_000913/734786‑734981

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