Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,587,370	C→T	T143I (ACT→ATT)  *248* (TAG→TAA)	yhhA → ugpQ ←	DUF2756 family protein glycerophosphodiester phosphodiesterase, cytosolic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,587,370	0	C	T	100.0%	99.5 / NA	32	T143I (ACT→ATT)  *248* (TAG→TAA)	yhhA ugpQ	DUF2756 family protein glycerophosphodiester phosphodiesterase, cytosolic
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/16);  total (16/16)

CAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCT  >  NC_000913/3587260‑3587502                                                                                                               |                                                                                                                                     cAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCa                                                                                                          >  1:736041/1‑139 (MQ=255)            cGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCa                                                                                               <  2:736041/139‑1 (MQ=255)              gATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATc                                                                                             <  1:3211/139‑1 (MQ=255)                        tGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATg                                                                                   <  2:941590/139‑1 (MQ=255)                            tCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAAt                                                                               >  1:806666/1‑139 (MQ=255)                               gcagcaTATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCca                                                                            >  2:59140/1‑139 (MQ=255)                                          tgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgcgc                                                                 >  1:168474/1‑139 (MQ=255)                                          tgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgcgc                                                                 >  1:408785/1‑139 (MQ=255)                                            ccGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCgcag                                                               <  2:861959/139‑1 (MQ=255)                                             cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGt                                                                                                                           >  1:285273/1‑78 (MQ=255)                                             cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGt                                                                                                                           <  2:285273/78‑1 (MQ=255)                                             cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                   <  1:704880/102‑1 (MQ=255)                                             cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAAt                                                                                                   >  2:704880/1‑102 (MQ=255)                                             cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgcgc                                                                 <  2:770831/136‑1 (MQ=255)                                             cGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgcgc                                                                 >  1:770831/1‑136 (MQ=255)                                               gagaggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgcgc                                                                 >  1:313516/1‑134 (MQ=255)                                               gagaggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgcgc                                                                 <  2:313516/134‑1 (MQ=255)                                                      gaAACGGCGACATGCTGAATCAGCCCAGCACGTCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACAccc                                                                       >  1:611496/1‑121 (MQ=255)                                                      gaAACGGCGACATGCTGAATCAGCCCAGCACGTCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACAccc                                                                       <  2:611496/121‑1 (MQ=255)                                                                         tCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  <  1:59140/139‑1 (MQ=255)                                                                             cccAGCACGCCGCAGCCTGACATTGCGTTGAAAATTATTGGGCCGAAAAGTTCGGACCAATCACGTCAATCGCATCGGTAGAGATGCAATCCACACACCAGCGCAGCAACTCTGCTGCGAGCTGGGGTTTATTGACGGt                              >  1:611588/1‑139 (MQ=255)                                                                             cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  <  1:557495/135‑1 (MQ=255)                                                                             cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  >  2:557495/1‑135 (MQ=255)                                                                               cAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGtat                            <  1:495546/139‑1 (MQ=255)                                                                                 gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  >  2:457982/1‑131 (MQ=255)                                                                                 gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  <  1:457982/131‑1 (MQ=255)                                                                                  cACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATaaa                         >  1:488107/1‑139 (MQ=255)                                                                                       cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  >  2:575208/1‑125 (MQ=255)                                                                                       cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                  <  1:575208/125‑1 (MQ=255)                                                                                               gACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAg            >  2:564665/1‑139 (MQ=255)                                                                                                      cGTTGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCg     <  2:488107/139‑1 (MQ=255)                                                                                                         tGAAAATTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCt  <  1:564665/139‑1 (MQ=255)                                                                                                               |                                                                                                                                     CAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCT  >  NC_000913/3587260‑3587502

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