Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,002,788	C→T	R59W (CGG→TGG)	xdhB →	putative xanthine dehydrogenase FAD‑binding subunit XdhB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,002,788	0	C	T	100.0%	99.8 / NA	29	R59W (CGG→TGG)	xdhB	putative xanthine dehydrogenase FAD‑binding subunit XdhB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/11);  total (18/11)

GTACTGATTCAGCTCCACCATCACAATGACCGTTATCGCCATATTGTTGATATTCATAATCTGGCGGAGCTGCGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTGCAACGACATTTACCCAGCTAATAGAAGATCCT  >  NC_000913/3002716‑3002865                                                                         |                                                                              gTACTGATTCAGCTCCACCATCACAATGACCGTTATCGCCATATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAAt                                                                        >  2:519393/1‑80 (MQ=255)  tACTGATTCAGCTCCACCATCACAATGACCGTTATCGCCATATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTa                                                                      >  1:519394/1‑81 (MQ=255)    cTGATTCAGCTCCACCATCACAATGACCGTTATCGCCATATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTAc                                                                     <  1:519351/80‑1 (MQ=255)             cTCCACCATCACAATGACCGTTATCGCCATATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCgg                                                             >  1:519385/1‑79 (MQ=255)                caccaTCACAATGACCGTTATCGCCATATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGAt                                                        >  1:519392/1‑81 (MQ=255)                       acaATGACCGTTATCGCCATATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGc                                                 <  2:519352/81‑1 (MQ=255)                         aaTGACCGTTATCGCCTTATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTa                                               >  1:519389/1‑81 (MQ=255)                          aTGACCGTTATCGCCATATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTAc                                              <  2:519353/81‑1 (MQ=255)                             aCCGTTATCGCCATATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGt                                            <  2:519354/80‑1 (MQ=255)                             aCCGTTATCGCCATATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTa                                           >  2:519390/1‑81 (MQ=255)                             aCCGTTATCGCCATATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTa                                           >  1:519388/1‑81 (MQ=255)                               cGTTATCGCCATATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATc                                         <  1:519355/81‑1 (MQ=255)                                  tATCGCCATATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGc                                      >  1:519386/1‑81 (MQ=255)                                      gCCATATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTg                                  >  1:519381/1‑81 (MQ=255)                                      gCCATATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTg                                  >  1:519405/1‑81 (MQ=255)                                      gCCATATTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTg                                  >  2:519391/1‑81 (MQ=255)                                         atatTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGAGGAAGATGGCTCGCTACGTATCGGCTCTGCaa                               <  1:519356/81‑1 (MQ=255)                                           atTGTTGATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTGCAAc                              >  2:519400/1‑80 (MQ=255)                                              gttgATATTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTGCAACGAc                           <  2:519357/80‑1 (MQ=255)                                                  atatTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTGCAACGACAttt                       <  1:519358/80‑1 (MQ=255)                                                  atatTCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTGCAACGACAttt                       <  2:519359/80‑1 (MQ=255)                                                     ttCATAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTGCAACGACATTTAcc                    >  1:519406/1‑80 (MQ=255)                                                         tAATCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTGCAACGACATTTACCCAGc                >  2:519397/1‑80 (MQ=255)                                                           aTCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTGCAACGACATTTACCCAGCTaa             >  1:519402/1‑81 (MQ=255)                                                           aTCTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTGCAACGACATTTACCCAGCTa              <  2:519360/80‑1 (MQ=255)                                                             cTGGCGGAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTGCAACGACATTTACCCAGCTaa             >  2:519396/1‑79 (MQ=255)                                                                  ggAGCTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTGCAACGACATTTACCCAGCTAATag          >  2:519395/1‑77 (MQ=255)                                                                      cTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTGCAACGACATTTACCCAGCTAATAGAAGATcc   >  2:519401/1‑80 (MQ=255)                                                                      cTGTGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTGCAACGACATTTACCCAGCTAATAGAAGATCCt  <  1:519361/81‑1 (MQ=255)                                                                         |                                                                              GTACTGATTCAGCTCCACCATCACAATGACCGTTATCGCCATATTGTTGATATTCATAATCTGGCGGAGCTGCGGGGAATTACGCTGGCGGAAGATGGCTCGCTACGTATCGGCTCTGCAACGACATTTACCCAGCTAATAGAAGATCCT  >  NC_000913/3002716‑3002865

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