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%	176.6 / NA	55	pseudogene (210/651 nt)	ybfG	pseudogene
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/30);  total (25/30)

AGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGT  >  NC_000913/716262‑716513                                                                                                                               |                                                                                                                              agCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTatca                                                                                                                   <  2:1068680/139‑1 (MQ=255)   cAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGTCCAGATCCTGCTATGCGTAATTAACCATTatcaat                                                                                                                 >  1:636804/1‑139 (MQ=255)   cAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTatcaat                                                                                                                 <  1:381450/139‑1 (MQ=255)    aaaaaTTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTatcaatc                                                                                                                <  1:413075/139‑1 (MQ=255)              ccTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTc                                                                                                      >  1:882200/1‑139 (MQ=255)                   gACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAg                                                                                                 <  1:772949/139‑1 (MQ=255)                          tCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGa                                                                                          <  1:744193/139‑1 (MQ=255)                               gATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATCAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAgc                                                                                     >  2:869471/1‑139 (MQ=255)                                 tGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAgcga                                                                                   >  2:335517/1‑138 (MQ=255)                                  gggTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCa                                                                                  <  2:282576/139‑1 (MQ=255)                                   ggTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCAc                                                                                 >  1:915308/1‑139 (MQ=255)                                     tAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCAccc                                                                               <  2:882200/139‑1 (MQ=255)                                        aGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGTCCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATc                                                                            <  2:636804/139‑1 (MQ=255)                                          cccGCACGCCTGCGAATTCAGGGGGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCgg                                                                          <  2:627589/139‑1 (MQ=255)                                                   cTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttgt                                                                 >  1:1061048/1‑139 (MQ=255)                                                         aTTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                           <  1:946751/139‑1 (MQ=255)                                                         aTTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                           <  2:1283546/139‑1 (MQ=255)                                                                          cAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTg                                          <  1:328440/139‑1 (MQ=255)                                                                                  tGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTg                                          <  2:1324286/131‑1 (MQ=255)                                                                                  tGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTg                                          >  1:1324286/1‑131 (MQ=255)                                                                                  tGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCg                                  >  1:654663/1‑139 (MQ=255)                                                                                      tctcGTGCCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTa                              >  2:1006371/1‑139 (MQ=255)                                                                                       ctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                           <  2:1296762/109‑1 (MQ=255)                                                                                       ctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                           >  1:1296762/1‑109 (MQ=255)                                                                                        tcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGtaattaac                                                                                                                           >  1:1384911/1‑44 (MQ=255)                                                                                        tcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGtaattaac                                                                                                                           <  2:1384911/44‑1 (MQ=255)                                                                                         cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCa                                                                                              <  2:1322360/72‑1 (MQ=255)                                                                                         cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCa                                                                                              >  1:1322360/1‑72 (MQ=255)                                                                                         cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAAt                                                    <  1:767602/114‑1 (MQ=255)                                                                                         cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAAt                                                    <  1:1409076/114‑1 (MQ=255)                                                                                         cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAAt                                                    >  2:1409076/1‑114 (MQ=255)                                                                                         cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAAt                                                    >  2:767602/1‑114 (MQ=255)                                                                                         cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAAt                           <  2:654663/139‑1 (MQ=255)                                                                                         cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAAt                           >  2:135702/1‑139 (MQ=255)                                                                                         cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAAt                           >  2:736413/1‑139 (MQ=255)                                                                                           tGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttgttg                                                               <  1:1013644/101‑1 (MQ=255)                                                                                           tGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttgttg                                                               >  2:1013644/1‑101 (MQ=255)                                                                                              cAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTg                                          <  2:785450/119‑1 (MQ=255)                                                                                              cAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTg                                          >  1:785450/1‑119 (MQ=255)                                                                                                aCCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTttatt                    >  1:1262192/1‑139 (MQ=255)                                                                                                 ccGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAgcg                                                                                    >  2:633148/1‑74 (MQ=255)                                                                                                 ccGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAgcg                                                                                    <  1:633148/74‑1 (MQ=255)                                                                                                  cGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTaa                  <  2:224188/139‑1 (MQ=255)                                                                                                      ggTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCgg              <  1:1290868/139‑1 (MQ=255)                                                                                                            cAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAagag        <  1:1006371/139‑1 (MQ=255)                                                                                                               aTCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTg                                                                        <  1:763100/72‑1 (MQ=255)                                                                                                               aTCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTg                                                                        >  2:763100/1‑72 (MQ=255)                                                                                                                     cTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                           <  1:186286/79‑1 (MQ=255)                                                                                                                     cTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                           >  2:186286/1‑79 (MQ=255)                                                                                                                         gCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAAt                                                    <  2:815613/82‑1 (MQ=255)                                                                                                                         gCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAAt                                                    >  1:815613/1‑82 (MQ=255)                                                                                                                         gCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCt                               <  1:593172/103‑1 (MQ=255)                                                                                                                         gCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGAGTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCt                               >  2:593172/1‑103 (MQ=255)                                                                                                                          cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCgg                                                                          >  1:148504/1‑59 (MQ=255)                                                                                                                          cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCgg                                                                          <  2:148504/59‑1 (MQ=255)                                                                                                                          cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGt  <  2:95882/131‑1 (MQ=255)                                                                                                                          cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGt  >  1:95882/1‑131 (MQ=255)                                                                                                                               |                                                                                                                              AGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGT  >  NC_000913/716262‑716513

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