Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,538,483	T→A	K95* (AAA→TAA)	nanM ←	N‑acetylneuraminate mutarotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,483	0	T	A	100.0%	74.4 / NA	21	K95* (AAA→TAA)	nanM	N‑acetylneuraminate mutarotase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (7/14);  total (7/14)

TCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCC  >  NC_000913/4538405‑4538552                                                                               |                                                                      tCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTAGc                                                                     <  2:913019/81‑1 (MQ=255)    aTTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTAGCCaa                                                                  >  2:913055/1‑81 (MQ=255)         aCCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTAGCCAATgccg                                                             <  1:913021/81‑1 (MQ=255)         aCCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTAGCCAATgcc                                                              <  1:913020/80‑1 (MQ=255)           ccAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTAGCCAATgccgcc                                                           <  1:913023/81‑1 (MQ=255)           ccAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTAGCCAATgccgcc                                                           <  2:913022/81‑1 (MQ=255)            cAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTAGCCAATGCCGCCa                                                          <  1:913024/81‑1 (MQ=255)                           gggTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTAGCCAATGCCGCCAAACACATACAGATTg                                           >  1:913056/1‑81 (MQ=255)                                 tATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTAGCCAATGCCGCCAAACACATACAGATTGCCATCa                                     <  2:913025/81‑1 (MQ=255)                                   tttGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTAGCCAATGCCGCCAAACACATACAGATTGCCATCAAt                                   <  2:913026/81‑1 (MQ=255)                                         tACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTAGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATg                              >  1:913053/1‑80 (MQ=255)                                         tACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTAGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATg                              >  1:913060/1‑80 (MQ=255)                                                 tAAATACCTGAGTCAAGCCCTCGCTGTTTTAGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTg                     >  1:913051/1‑81 (MQ=255)                                                  aaaTACCTGAGTCAAGCCCTCGCTGTTTTAGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGc                    >  1:913061/1‑81 (MQ=255)                                                    aTACCTGAGTCAAGCCCTCGCTGTTTTAGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCtt                  <  1:913027/81‑1 (MQ=255)                                                       ccTGAGTCAAGCCCTCGCTGTTTTAGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGAt               <  1:913029/81‑1 (MQ=255)                                                       ccTGAGTCAAGCCCTCGCTGTTTTAGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGAt               <  2:913028/81‑1 (MQ=255)                                                        cTGAGTCAAGCCCTCGCTGTTTTAGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGAtc              <  1:913030/81‑1 (MQ=255)                                                          gAGTCAAGCCCTCGCTGTTTTAGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGAtctc            >  2:913058/1‑81 (MQ=255)                                                               aaGCCCTCGCTGTTTTAGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTgg        <  1:913031/80‑1 (MQ=255)                                                                    cTCGCTGTTTTAGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGcc  <  1:913032/81‑1 (MQ=255)                                                                               |                                                                      TCAATTTAACCCAACTATTGGTTTTGGGGTTGTATTTGTGTACGTCATTAAATACCTGAGTCAAGCCCTCGCTGTTTTTGCCAATGCCGCCAAACACATACAGATTGCCATCAATAAATGCAGAGGTTGCTTGATCTCTTGGTCCGCC  >  NC_000913/4538405‑4538552

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