Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	854,765	G→A	*373* (TAG→TAA)	ybiR →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	854,765	0	G	A	100.0%	85.8 / NA	28	*373* (TAG→TAA)	ybiR	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/13);  total (15/13)

TGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGC  >  NC_000913/854676‑854896                                                                                          |                                                                                                                                    tggtggCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAg                                                                                    >  2:587224/1‑139 (MQ=255)  ggtggCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTa                                                                                                          <  2:650910/116‑1 (MQ=255)  ggtggCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTa                                                                                                          >  1:650910/1‑116 (MQ=255)                    aTTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGa                                                                 >  2:125858/1‑139 (MQ=255)                          tACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGt                                                                                                                                  <  2:505876/68‑1 (MQ=255)                          tACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGt                                                                                                                                  >  1:505876/1‑68 (MQ=255)                             cGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGc                                                        >  1:236920/1‑139 (MQ=255)                               aTGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCAt                                                      >  2:769768/1‑139 (MQ=255)                               aTGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCAt                                                      >  2:71460/1‑139 (MQ=255)                               aTGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCAt                                                      >  1:43079/1‑139 (MQ=255)                                   tgttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCg                                                  <  1:402023/139‑1 (MQ=255)                                      tgtgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAg                                               <  1:71460/139‑1 (MQ=255)                                      tgtgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAg                                               <  2:236920/139‑1 (MQ=255)                                               cgTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTa                                      >  2:474237/1‑139 (MQ=255)                                                                   ttagttaTACTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGaacaac                  >  1:455127/1‑139 (MQ=255)                                                                        tataCTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGAt                                                                          >  1:243530/1‑78 (MQ=255)                                                                        tataCTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGAt                                                                          <  2:243530/78‑1 (MQ=255)                                                                        tataCTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGt             <  2:43079/139‑1 (MQ=255)                                                                        tataCTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGt             >  2:607469/1‑139 (MQ=255)                                                                        tataCTCCCGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGt             <  1:474237/139‑1 (MQ=255)                                                                               ccGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                     <  2:256661/60‑1 (MQ=255)                                                                               ccGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                     >  1:256661/1‑60 (MQ=255)                                                                               ccGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATc      <  1:772228/139‑1 (MQ=255)                                                                                cGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGaacaac                  >  2:580247/1‑126 (MQ=255)                                                                                cGGCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGaacaac                  <  1:580247/126‑1 (MQ=255)                                                                                 ggCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTAt                                                                                                         >  1:496650/1‑38 (MQ=255)                                                                                 ggCCAACTAAGGTCTGGCAAATAAAAAAGGCGGATTAt                                                                                                         <  2:496650/38‑1 (MQ=255)                                                                                   ccAACTAAGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGc  <  1:277727/139‑1 (MQ=255)                                                                                          |                                                                                                                                    TGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGC  >  NC_000913/854676‑854896

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