Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	716,388	C→T	pseudogene (210/651 nt)	ybfG ←	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	716,388	0	C	T	100.0%	244.4 / NA	76	pseudogene (210/651 nt)	ybfG	pseudogene
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/36);  total (40/36)

CATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTG  >  NC_000913/716255‑716523                                                                                                                                      |                                                                                                                                        caTTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAAcc                                                                                                                                    >  1:99767/1‑139 (MQ=255)          cAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTatcaat                                                                                                                           <  2:99767/139‑1 (MQ=255)                        gCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCgg                                                                                                             >  1:217893/1‑139 (MQ=255)                        gCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCgg                                                                                                             >  2:492758/1‑139 (MQ=255)                         cGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGa                                                                                                            >  2:50609/1‑139 (MQ=255)                            cATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGcc                                                                                                         <  2:217893/139‑1 (MQ=255)                                 tCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                                              <  1:1449202/129‑1 (MQ=255)                                 tCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                                              >  2:1449202/1‑129 (MQ=255)                                        tGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAgcgc                                                                                             <  1:50609/139‑1 (MQ=255)                                              aaGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCAt                                                                                       <  2:90171/139‑1 (MQ=255)                                                     cACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGAt                                                                                >  1:1421518/1‑139 (MQ=255)                                                     cACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGAt                                                                                >  1:1443055/1‑139 (MQ=255)                                                              gAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTc                                                                       <  1:492758/139‑1 (MQ=255)                                                                aTTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                                     <  2:1421518/139‑1 (MQ=255)                                                                aTTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                                     <  2:1442091/139‑1 (MQ=255)                                                                     gggCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCgg                                                                                    >  2:1360205/1‑119 (MQ=255)                                                                     gggCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCgg                                                                                    <  1:1360205/119‑1 (MQ=255)                                                                             tgtCCAACTTACTGGATCTCTTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAg                                                                                                           <  2:1407771/88‑1 (MQ=255)                                                                             tgtCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAg                                                                                                           >  1:1407771/1‑88 (MQ=255)                                                                                              ctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTatca                                                                                                                             <  1:703788/53‑1 (MQ=255)                                                                                              ctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTatca                                                                                                                             >  2:1284896/1‑53 (MQ=255)                                                                                              ctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTatca                                                                                                                             <  1:1284896/53‑1 (MQ=255)                                                                                              ctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTatca                                                                                                                             >  2:703788/1‑53 (MQ=255)                                                                                              ctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAg                                                                                                           >  2:206333/1‑71 (MQ=255)                                                                                              ctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAg                                                                                                           <  1:206333/71‑1 (MQ=255)                                                                                              ctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCa                                                                                            <  2:1609480/86‑1 (MQ=255)                                                                                              ctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCa                                                                                            >  1:1609480/1‑86 (MQ=255)                                                                                              ctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGc                                          >  1:1398267/1‑136 (MQ=255)                                                                                              ctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGc                                          <  2:1398267/136‑1 (MQ=255)                                                                                              ctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTaa                                       >  1:869641/1‑139 (MQ=255)                                                                                                cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCCTTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTACTTGATCAAGTGCGGTTTCGGCTAAAt                                     >  1:464386/1‑139 (MQ=255)                                                                                                cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                                              <  2:1449628/66‑1 (MQ=255)                                                                                                cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                                              >  1:1449628/1‑66 (MQ=255)                                                                                                cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGa                                                                                                    <  1:1284692/76‑1 (MQ=255)                                                                                                cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGa                                                                                                    >  2:1284692/1‑76 (MQ=255)                                                                                                cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttg                                                                            <  2:17550/100‑1 (MQ=255)                                                                                                cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttg                                                                            >  1:17550/1‑100 (MQ=255)                                                                                                cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAAt                                     >  2:550837/1‑139 (MQ=255)                                                                                                 gTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAAGCATTATCAATCAGCAGACCTCCg                                                                                                              >  2:1042976/1‑65 (MQ=255)                                                                                                 gTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                                              <  1:1042976/65‑1 (MQ=255)                                                                                                  tGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                                              <  2:1030831/64‑1 (MQ=255)                                                                                                  tGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCt                                   >  1:398165/1‑139 (MQ=255)                                                                                                  tGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAAACATTATCAATCAGCAGACCTCCg                                                                                                              >  1:1030831/1‑64 (MQ=255)                                                                                                   ggCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAgcg                                                                                              <  2:939461/79‑1 (MQ=255)                                                                                                   ggCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAgcg                                                                                              >  1:939461/1‑79 (MQ=255)                                                                                                      aaCCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATcag                                                                                                                        >  1:1667873/1‑50 (MQ=255)                                                                                                      aaCCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATcag                                                                                                                        <  2:1667873/50‑1 (MQ=255)                                                                                                      aaCCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                                              >  1:927134/1‑60 (MQ=255)                                                                                                      aaCCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                                              <  2:927134/60‑1 (MQ=255)                                                                                                       aCCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTttatt                              >  1:1304010/1‑139 (MQ=255)                                                                                                        ccGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAAt                                                              >  2:1681832/1‑106 (MQ=255)                                                                                                        ccGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAAt                                                              <  1:1681832/106‑1 (MQ=255)                                                                                                           ggggTGACCAGATCCTGCTATTCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCgg                                                                                    <  2:714697/81‑1 (MQ=255)                                                                                                           ggggTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCgg                                                                                    >  1:714697/1‑81 (MQ=255)                                                                                                                gACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttgtt                                                                          >  1:1297187/1‑86 (MQ=255)                                                                                                                gACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttgtt                                                                          <  2:1297187/86‑1 (MQ=255)                                                                                                                  ccAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAaga                   <  2:1304010/139‑1 (MQ=255)                                                                                                                   cAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAagag                  <  2:398165/139‑1 (MQ=255)                                                                                                                      aTCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCgg                                                                                    <  2:946718/70‑1 (MQ=255)                                                                                                                      aTCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCgg                                                                                    >  1:946718/1‑70 (MQ=255)                                                                                                                           gCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCAcc                                                                                          >  1:832488/1‑59 (MQ=255)                                                                                                                           gCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCAcc                                                                                          <  2:832488/59‑1 (MQ=255)                                                                                                                           gCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                                     <  1:654922/80‑1 (MQ=255)                                                                                                                           gCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                                     >  2:654922/1‑80 (MQ=255)                                                                                                                                gCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCa                                                                                            >  1:817306/1‑52 (MQ=255)                                                                                                                                gCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCa                                                                                            <  2:817306/52‑1 (MQ=255)                                                                                                                                gCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTg                                                    >  2:385090/1‑92 (MQ=255)                                                                                                                                gCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTg                                                    <  1:385090/92‑1 (MQ=255)                                                                                                                                gCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGttt     >  2:1427969/1‑139 (MQ=255)                                                                                                                                gCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGttt     >  2:843309/1‑139 (MQ=255)                                                                                                                                 cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCAt                                                                                       >  1:978578/1‑56 (MQ=255)                                                                                                                                 cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCAt                                                                                       <  2:978578/56‑1 (MQ=255)                                                                                                                                 cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGtt                                               >  1:126587/1‑96 (MQ=255)                                                                                                                                 cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGtt                                               <  2:126587/96‑1 (MQ=255)                                                                                                                                 cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTaaa                                      <  1:449478/105‑1 (MQ=255)                                                                                                                                 cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTaaa                                      >  2:449478/1‑105 (MQ=255)                                                                                                                                   taattaacCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTtgtg  >  1:1437663/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                        CATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTGTG  >  NC_000913/716255‑716523

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