Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,612,540	T→C	W226R (TGG→CGG)	yhhT →	UPF0118 family putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,612,540	0	T	C	100.0%	102.6 / NA	30	W226R (TGG→CGG)	yhhT	UPF0118 family putative transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/14);  total (16/14)

GAATAATCCACAGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCATTGTTTTTAGTGGTCCATATGGT  >  NC_000913/3612392‑3612685                                                                                                                                                     |                                                                                                                                                  gAATAATCCACAGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCgg                                                                                                                                                 >  2:550889/1‑151 (MQ=255)     aatCCACAGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCa                                                                                                                                              <  2:385299/150‑1 (MQ=255)         cacaGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACt                                                                                                                                         <  1:350675/151‑1 (MQ=255)         cacaGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACt                                                                                                                                         >  1:453432/1‑151 (MQ=255)                   caTCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTg                                                                                                                               >  1:54921/1‑151 (MQ=255)                     tCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCt                                                                                                                             >  2:300647/1‑151 (MQ=255)                             ttACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACg                                                                                                                     >  2:290411/1‑151 (MQ=255)                                            aaaGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGgcg                                                                                                      >  1:656710/1‑151 (MQ=255)                                                 cGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTa                                                                                                 <  1:524912/151‑1 (MQ=255)                                                  gTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTa                                                                                                 >  1:141484/1‑150 (MQ=255)                                                          cTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTccgc                                                                                         >  1:644748/1‑150 (MQ=255)                                                            aTCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTccgcc                                                                                        <  2:483923/149‑1 (MQ=255)                                                               ttGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTAcc                                                                                   >  2:55593/1‑151 (MQ=255)                                                                      gAAGACGCTACTCAGTTTATGGACAGGTGCCATCTTCTGGCGGGGGCTGGAGTTGATGGGGGTGCAGTTTGCGTTGATGCGGGCAGTACTGGCTTTTTGGCTCAACTACGTGCCCAATACCGGCGGGGTAATTTCCGCCGTACCGCCAAt                                                                             <  1:400496/150‑1 (MQ=255)                                                                              tACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGt                                                                    >  1:494366/1‑151 (MQ=255)                                                                                  cAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGtgct                                                                 <  1:580338/150‑1 (MQ=255)                                                                                  cAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGtgct                                                                 <  2:44918/150‑1 (MQ=255)                                                                                   aGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGtgctgc                                                               >  2:64565/1‑151 (MQ=255)                                                                                        aTGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGttt                                                          <  1:605645/151‑1 (MQ=255)                                                                                        aTGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGATTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGttt                                                          >  1:350940/1‑151 (MQ=255)                                                                                                  gtCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGttt                                                <  1:32422/151‑1 (MQ=255)                                                                                                    cATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTAc                                              <  2:181191/151‑1 (MQ=255)                                                                                                        gTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAAt                                          >  1:51455/1‑151 (MQ=255)                                                                                                              ctggggctggAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTc                                    >  2:592503/1‑151 (MQ=255)                                                                                                                     tggAGCTGATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCgg                             >  2:630829/1‑151 (MQ=255)                                                                                                                            gATGGGGGTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCATTg                      <  1:577777/151‑1 (MQ=255)                                                                                                                               gggggTGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCATTGtt                    <  1:630819/150‑1 (MQ=255)                                                                                                                                    tGCAGTTTGCGCTGATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCATTGTTTTTAGt              >  2:606976/1‑151 (MQ=255)                                                                                                                                             cgcTGATGCGGGCAGTACTGGCGTTGTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCATTGTTTTTAGTGGTCCata      <  2:133708/150‑1 (MQ=255)                                                                                                                                                 gATGCGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCATTGTTTTTAGTGGTCCATATGGt  <  2:430579/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GAATAATCCACAGATTCACATCGCGGGATTACACCGCGCACTTAAAGGCGTTTCGCACTATCTTGCATTGAAGACGCTACTCAGTTTATGGACAGGTGTCATCGTCTGGCTGGGGCTGGAGCTGATGGGGGTGCAGTTTGCGCTGATGTGGGCAGTACTGGCGTTTTTGCTCAACTACGTGCCCAATATCGGCGCGGTAATTTCCGCCGTACCGCCAATGATTCAGGTGCTGCTGTTTAATGGTGTTTACGAATGTATTCTGGTCGGCGCATTGTTTTTAGTGGTCCATATGGT  >  NC_000913/3612392‑3612685

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