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

AAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTG  >  NC_000913/716265‑716521                                                                                                                            |                                                                                                                                      aaaaaTTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTatcaatc                                                                                                                        >  2:1443293/1‑139 (MQ=255)             tGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                                            >  1:201832/1‑139 (MQ=255)               cGACATCGACAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGTAGACCTCCAGa                                                                                                          >  2:1048150/1‑139 (MQ=255)                       tCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTc                                                                                                              <  2:851960/127‑1 (MQ=255)                       tCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTc                                                                                                              >  1:851960/1‑127 (MQ=255)                                    aaGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCAt                                                                                     <  1:1613475/139‑1 (MQ=255)                                          gCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATa                                                                                                    >  2:876280/1‑118 (MQ=255)                                          gCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATa                                                                                                    <  1:876280/118‑1 (MQ=255)                                           cACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGAt                                                                              <  1:100414/139‑1 (MQ=255)                                                cTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttgt                                                                         >  2:1382781/1‑139 (MQ=255)                                                        tCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCt                                                                 <  1:422752/139‑1 (MQ=255)                                                                 tgtgTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGAt                                                        <  2:1592308/139‑1 (MQ=255)                                                                        aaCTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGc                                                 >  1:1160360/1‑139 (MQ=255)                                                                               tGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCgg                                                                                  >  2:1381037/1‑99 (MQ=255)                                                                               tGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCgg                                                                                  <  1:1381037/99‑1 (MQ=255)                                                                                   tctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGTGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGc                                        >  2:1349814/1‑137 (MQ=255)                                                                                   tctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGTGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGc                                        <  1:634304/137‑1 (MQ=255)                                                                                   tctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGTGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGc                                        <  1:1349814/137‑1 (MQ=255)                                                                                   tctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGTGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGc                                        >  2:634304/1‑137 (MQ=255)                                                                                     tcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGtt                                             <  1:1353119/130‑1 (MQ=255)                                                                                     tcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGtt                                             >  2:1353119/1‑130 (MQ=255)                                                                                     tcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTaaa                                    >  2:414953/1‑139 (MQ=255)                                                                                     tcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTaaa                                    <  1:1443293/139‑1 (MQ=255)                                                                                      cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCCTTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                                   <  2:958750/107‑1 (MQ=255)                                                                                      cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                                   >  1:958750/1‑107 (MQ=255)                                                                                      cGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAAt                                   >  2:1296951/1‑139 (MQ=255)                                                                                          gCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCttt                               <  2:1160360/139‑1 (MQ=255)                                                                                            aaCCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTttat                             <  2:201832/139‑1 (MQ=255)                                                                                             aCCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTttatt                            >  1:1320571/1‑139 (MQ=255)                                                                                                 ggggTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCAt                                                                                     >  2:130323/1‑78 (MQ=255)                                                                                                 ggggTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCAt                                                                                     <  1:130323/78‑1 (MQ=255)                                                                                                            aTCCTGCTATGCGTAATTAACCATTATCAATCAGCAGAcc                                                                                                                >  2:21575/1‑40 (MQ=255)                                                                                                            aTCCTGCTATGCGTAATTAACCATTATCAATCAGCAGAcc                                                                                                                <  1:21575/40‑1 (MQ=255)                                                                                                                 gCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTAt        >  1:664623/1‑139 (MQ=255)                                                                                                                  cTATGCGTAATTAACCATTATCAATCAGCAGACCTCCgg                                                                                                           >  2:1048101/1‑39 (MQ=255)                                                                                                                  cTATGCGTAATTAACCATTATCAATCAGCAGACCTCCgg                                                                                                           <  1:1048101/39‑1 (MQ=255)                                                                                                                  cTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTtgttgt                                                                      <  1:449142/76‑1 (MQ=255)                                                                                                                  cTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTtgttgt                                                                      >  2:449142/1‑76 (MQ=255)                                                                                                                      gCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCAt                                                                                     >  1:807167/1‑57 (MQ=255)                                                                                                                      gCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCAt                                                                                     <  2:807167/57‑1 (MQ=255)                                                                                                                       cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTttatta                           >  2:98654/1‑114 (MQ=255)                                                                                                                       cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTttatta                           <  1:98654/114‑1 (MQ=255)                                                                                                                       cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTtg  >  1:1460765/1‑139 (MQ=255)                                                                                                                            |                                                                                                                                      AAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAGAGTTCGGTATGGTTTG  >  NC_000913/716265‑716521

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