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

CACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAG  >  NC_000913/716252‑716505                                                                                                                                         |                                                                                                                      caccaTTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGtaatta                                                                                                                     <  1:263353/139‑1 (MQ=255)              aaaaaTTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTatca                                                                                                           <  2:789680/136‑1 (MQ=255)              aaaaaTTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTatca                                                                                                           >  1:789680/1‑136 (MQ=255)                        ccTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTc                                                                                              >  2:525620/1‑139 (MQ=255)                          tGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                            <  2:163871/139‑1 (MQ=255)                          tGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                            <  2:504411/139‑1 (MQ=255)                          tGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCg                                                                                            >  2:796595/1‑139 (MQ=255)                           gCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCgg                                                                                           <  1:91685/139‑1 (MQ=255)                           gCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCgg                                                                                           <  2:222940/139‑1 (MQ=255)                           gCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCgg                                                                                           <  2:300895/139‑1 (MQ=255)                           gCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCgg                                                                                           <  1:25723/139‑1 (MQ=255)                            cGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGa                                                                                          >  2:932449/1‑139 (MQ=255)                                    tCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGa                                                                                  <  1:932449/139‑1 (MQ=255)                                             ggTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCAc                                                                         >  1:519940/1‑139 (MQ=255)                                                  aGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATc                                                                    >  1:532140/1‑139 (MQ=255)                                                    cccGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCgg                                                                  >  1:107534/1‑139 (MQ=255)                                                       gCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGa                                                               >  2:358515/1‑139 (MQ=255)                                                        cACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGtaattaac                                                                                                                   >  2:1049340/1‑86 (MQ=255)                                                        cACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGtaattaac                                                                                                                   <  1:1049340/86‑1 (MQ=255)                                                        cACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGAt                                                              >  2:88648/1‑139 (MQ=255)                                                            ccTGCGAATTCAGGGCGTTGTGTCCATCTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTATCCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttg                                                          <  2:394107/139‑1 (MQ=255)                                                              tGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttgtt                                                        <  1:935951/139‑1 (MQ=255)                                                               gCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttgttg                                                       <  1:88648/139‑1 (MQ=255)                                                                cGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTtgttgt                                                      <  2:107534/139‑1 (MQ=255)                                                                   aTTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCGGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                   <  2:147199/139‑1 (MQ=255)                                                                   aTTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                   <  1:238598/139‑1 (MQ=255)                                                                   aTTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                   <  2:532140/139‑1 (MQ=255)                                                                     tCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCt                                                 <  1:796595/139‑1 (MQ=255)                                                                       aGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGt                                               <  2:362802/139‑1 (MQ=255)                                                                         ggCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTaa                                             >  2:146546/1‑139 (MQ=255)                                                                              tgtgTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttgttg                                                       >  1:504078/1‑124 (MQ=255)                                                                              tgtgTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttgttg                                                       <  2:504078/124‑1 (MQ=255)                                                                                               aTCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGTGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCt                       <  1:164368/139‑1 (MQ=255)                                                                                                 ctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                   <  2:913899/109‑1 (MQ=255)                                                                                                 ctcGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAc                                                   >  1:913899/1‑109 (MQ=255)                                                                                                    gTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCt                       >  1:467725/1‑134 (MQ=255)                                                                                                    gTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCt                       <  2:467725/134‑1 (MQ=255)                                                                                                     tGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAAGCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAAGTGATAAAGGGGGGTTTCGGCTAAAGCt                 >  2:846525/1‑139 (MQ=255)                                                                                                          aCCGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTttatt            >  2:305026/1‑139 (MQ=255)                                                                                                            cGGGGGTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTaa          >  1:287086/1‑139 (MQ=255)                                                                                                              ggggTGACCAGATCCTGCTATGCGTAATTAACCATTatca                                                                                                           >  1:1055098/1‑40 (MQ=255)                                                                                                              ggggTGACCAGATCCTGCTATGCGTAATTAACCATTatca                                                                                                           <  2:1055098/40‑1 (MQ=255)                                                                                                                 gTGACCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGt     <  2:287086/139‑1 (MQ=255)                                                                                                                    aCCAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAag  >  2:541225/1‑139 (MQ=255)                                                                                                                      cAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCgg                                                                                           <  2:473202/48‑1 (MQ=255)                                                                                                                      cAGATCCTGCTATGCGTAATTAACCATTATCAATCAGCAGACCTCCgg                                                                                           >  1:473202/1‑48 (MQ=255)                                                                                                                                  tGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTAGTTGTCAc                                                   >  1:993391/1‑76 (MQ=255)                                                                                                                                  tGCGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTAGTTGTCAc                                                   <  2:993391/76‑1 (MQ=255)                                                                                                                                    cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTttat             >  2:478616/1‑112 (MQ=255)                                                                                                                                    cGTAATTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTttat             <  1:478616/112‑1 (MQ=255)                                                                                                                                         |                                                                                                                      CACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAAAGTGCGGTTTCGGCTAAATCTTTATTAACCGGTAAG  >  NC_000913/716252‑716505

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