Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,136,921	(G)5→4	coding (814/906 nt)	yijE →	cystine exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,136,917	0	G	.	100.0%	174.7 / NA	41	coding (810/906 nt)	yijE	cystine exporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (15/26);  total (15/26)

GCCAGTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCG  >  NC_000913/4136837‑4136993                                                                                 |                                                                             gccAGTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCg                                                                              >  2:771446/1‑81 (MQ=255)  ccAGTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCgg                                                                             <  2:771395/81‑1 (MQ=255)     gTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCggggc                                                                          <  1:771397/81‑2 (MQ=255)     gTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCggggc                                                                          >  1:771462/1‑80 (MQ=255)       aTTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCggggccg                                                                        >  2:771440/1‑78 (MQ=255)        ttGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCggggccgt                                                                       >  2:771444/1‑77 (MQ=255)         tGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGtt                                                                     >  2:771457/1‑81 (MQ=255)          gCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGtt                                                                     >  2:771451/1‑80 (MQ=255)            cAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTg                                                                    >  1:771446/1‑79 (MQ=255)               cTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGt                                                               <  1:771398/81‑1 (MQ=255)                    gCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCgg                                                          <  1:771400/81‑1 (MQ=255)                    gCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCgg                                                          <  2:771398/81‑1 (MQ=255)                     cacaCTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGt                                                         <  2:771400/81‑1 (MQ=255)                              cGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTattgtgctg                                                <  2:771401/81‑1 (MQ=255)                               gTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTattgtgctga                                               >  1:771449/1‑81 (MQ=255)                               gTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTattgtgctga                                               <  2:771402/81‑1 (MQ=255)                                   ccGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTattgtgctgattgt                                           >  2:771453/1‑81 (MQ=255)                                   ccGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTattgtgctgattg                                            <  2:771403/80‑1 (MQ=255)                                         gcggcgTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTgg                                      >  1:771450/1‑80 (MQ=255)                                          cggcgTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGc                                     <  1:771406/80‑1 (MQ=255)                                          cggcgTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGc                                     <  1:771407/80‑1 (MQ=255)                                          cggcgTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGc                                     <  2:771404/80‑1 (MQ=255)                                           ggcgTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCAc                                   >  2:771447/1‑81 (MQ=255)                                            gcgTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACt                                  <  2:771407/81‑1 (MQ=255)                                             cgTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTg                                 <  2:771408/81‑1 (MQ=255)                                                    ttttCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCt                            <  2:771409/79‑1 (MQ=255)                                                    ttttCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTgg                          >  1:771466/1‑81 (MQ=255)                                                      ttCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGt                         <  1:771412/80‑1 (MQ=255)                                                      ttCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGt                         <  2:771410/80‑1 (MQ=255)                                                       tCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTg                        <  1:771416/80‑1 (MQ=255)                                                       tCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTg                        <  1:771414/80‑1 (MQ=255)                                                       tCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTg                        <  2:771412/80‑1 (MQ=255)                                                       tCCTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTg                        <  1:771415/80‑1 (MQ=255)                                                        ccTGGTGGCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGa                       <  1:771417/80‑1 (MQ=255)                                                           ggtggCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCg                   >  1:771469/1‑81 (MQ=255)                                                           ggtggCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCg                   >  2:771455/1‑81 (MQ=255)                                                            gtggCTGCTCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGcc                    <  2:771417/79‑1 (MQ=255)                                                               gctgctCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAg               <  1:771419/81‑1 (MQ=255)                                                                   ctCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGaaaa           <  2:771419/81‑1 (MQ=255)                                                                   ctCGGCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGaa             <  2:771420/79‑1 (MQ=255)                                                                      ggCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGaaaa           >  1:771465/1‑78 (MQ=255)                                                                      ggCGAGAATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGaaaa           >  1:771463/1‑78 (MQ=255)                                                                           gaATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGcc   >  2:771456/1‑81 (MQ=255)                                                                           gaATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGcc   <  1:771422/81‑1 (MQ=255)                                                                           gaATCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGcc   >  2:771459/1‑81 (MQ=255)                                                                            aaTCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCg  <  1:771423/81‑1 (MQ=255)                                                                            aaTCC‑GGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCg  >  2:771458/1‑81 (MQ=255)                                                                                 |                                                                             GCCAGTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCG  >  NC_000913/4136837‑4136993

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