Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	716,122	T→G	pseudogene (476/651 nt)	ybfG ←	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	716,122	0	T	G	100.0%	213.2 / NA	56	pseudogene (476/651 nt)	ybfG	pseudogene
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (34/22);  total (34/22)

TACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAA  >  NC_000913/715975‑716267                                                                                                                                                    |                                                                                                                                                  tACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCAt                                                                                                                                                >  2:587842/1‑151 (MQ=255)    tGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCAt                                                                                                                                                >  2:190906/1‑148 (MQ=255)     gCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAg                                                                                                                                              <  2:101611/149‑1 (MQ=255)      cATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGa                                                                                                                                             >  1:352279/1‑149 (MQ=255)         ttCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGt                                                                                                                                         >  2:91643/1‑150 (MQ=255)           cccAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTAt                                                                                                                                       <  2:528632/150‑1 (MQ=255)            ccAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTTCGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCt                                                                                                                                     <  1:359390/151‑1 (MQ=255)              agagTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTa                                                                                                                                    >  1:18324/1‑150 (MQ=255)              agagTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTa                                                                                                                                    >  2:607688/1‑150 (MQ=255)                    cAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCAcc                                                                                                                             >  1:527608/1‑151 (MQ=255)                       tAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCgg                                                                                                                           >  1:809460/1‑150 (MQ=255)                       tAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCgg                                                                                                                           >  2:488437/1‑150 (MQ=255)                        aaGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGAt                                                                                                                         >  1:560615/1‑151 (MQ=255)                           cGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAAt                                                                                                                      >  2:469638/1‑151 (MQ=255)                                tAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGccc                                                                                                                  >  2:477345/1‑150 (MQ=255)                                 aaCTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCa                                                                                                                 <  1:37291/150‑1 (MQ=255)                                  aCTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACt                                                                                                               >  1:93679/1‑151 (MQ=255)                                    ttttCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTg                                                                                                              >  2:300986/1‑150 (MQ=255)                                      ttcttgttcttgCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTga                                                                                                             >  1:633588/1‑149 (MQ=255)                                            ttcttgCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATACCCACTGAAGAAACt                                                                                                      >  2:742471/1‑150 (MQ=255)                                                tgCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGa                                                                                                  <  2:428411/150‑1 (MQ=255)                                                 gCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAat                                                                                                <  2:383991/151‑1 (MQ=255)                                                   ttcttcACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAatat                                                                                              >  1:706188/1‑151 (MQ=255)                                                            ttGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACtt                                                                                       <  2:429436/149‑1 (MQ=255)                                                            ttGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTc                                                                                      <  1:339379/150‑1 (MQ=255)                                                                aCCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGAc                                                                                 >  1:408488/1‑151 (MQ=255)                                                                 cccATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACt                                                                                >  1:67388/1‑151 (MQ=255)                                                                 cccATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACt                                                                                >  1:447674/1‑151 (MQ=255)                                                                    aTAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTc                                                                             <  2:122395/151‑1 (MQ=255)                                                                            gggTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCa                                                                     <  1:162514/151‑1 (MQ=255)                                                                              gTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCGATTCAGc                                                                   <  1:818471/151‑1 (MQ=255)                                                                              gTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGc                                                                   >  2:581842/1‑151 (MQ=255)                                                                               tAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCg                                                                  >  2:264200/1‑151 (MQ=255)                                                                                aGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGt                                                                 <  2:207884/151‑1 (MQ=255)                                                                                     gTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTg                                                            <  1:867049/151‑1 (MQ=255)                                                                                      tGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTgg                                                           >  1:619982/1‑151 (MQ=255)                                                                                         aGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGgcgc                                                        >  1:453344/1‑151 (MQ=255)                                                                                         aGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGATGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGgcgc                                                        >  2:518261/1‑151 (MQ=255)                                                                                           gCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCa                                                       >  2:560217/1‑150 (MQ=255)                                                                                                  ttGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAaataat                                                >  1:818897/1‑150 (MQ=255)                                                                                                        aTAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCa                                         <  2:126284/151‑1 (MQ=255)                                                                                                        aTAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCa                                         <  1:867497/151‑1 (MQ=255)                                                                                                        aTAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCa                                         <  1:306388/151‑1 (MQ=255)                                                                                                                 ttaCGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATc                                >  2:371796/1‑151 (MQ=255)                                                                                                                 ttaCGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATc                                >  2:388554/1‑151 (MQ=255)                                                                                                                   aCGTCTTCTCGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCtt                              <  2:66437/151‑1 (MQ=255)                                                                                                                           tcGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCa                         >  2:312492/1‑148 (MQ=255)                                                                                                                            cGCCATTGCTGCCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGa                      <  2:18325/150‑1 (MQ=255)                                                                                                                                  tgctgcCATCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGccac               >  1:603184/1‑151 (MQ=255)                                                                                                                                         aTCGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAag        <  2:786705/151‑1 (MQ=255)                                                                                                                                           cGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAag        >  1:697770/1‑149 (MQ=255)                                                                                                                                           cGTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAag        >  1:864406/1‑149 (MQ=255)                                                                                                                                            gTACATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAaga       <  2:237400/149‑1 (MQ=255)                                                                                                                                              aCATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCa    <  1:653396/150‑1 (MQ=255)                                                                                                                                              aCATTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCa    <  2:352291/150‑1 (MQ=255)                                                                                                                                                aTTTGCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCaaa  >  1:371402/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TACTGCATTTCCCAGAGTTCAGTAAGCGAAATAACTTTTCTTGTTCTTGCTTCTTCACTTTGTACCCATAAGAAAGGGTAGAACGTGAAGGCTTCCGTTGCCGATAATTTTATTACGTCTTCTCGCCATTGCTGCCATCGTACATTTTCATAGAAGGTATCTAAATCACCGGATAATGCCCACTGAAGAAACTCAGAATATCCAACTTCCAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAA  >  NC_000913/715975‑716267

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑