Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,068,157	T→C	intergenic (‑77/+16)	yggE ← / ← argO	uncharacterized protein YggE/L‑arginine exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,068,157	0	T	C	100.0%	118.6 / NA	34	intergenic (‑77/+16)	yggE/argO	uncharacterized protein YggE/L‑arginine exporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/18);  total (16/18)

AATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTTTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAAT  >  NC_000913/3068039‑3068265                                                                                                                       |                                                                                                             aaTACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAgcg                                                                                                <  1:81819/133‑1 (MQ=255) aaTACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAgcg                                                                                                <  2:367937/133‑1 (MQ=255) aaTACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAgcg                                                                                                >  1:484890/1‑133 (MQ=255)       cATTAATGGCGCCAGGGCGATATCTTTGAAATTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAAc                                                                                          <  1:229487/133‑1 (MQ=255)           aaTGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGaa                                                                                      <  2:108807/133‑1 (MQ=255)             tGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACa                                                                                    >  2:444668/1‑133 (MQ=255)             tGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACa                                                                                    <  1:337850/133‑1 (MQ=255)               ccgccAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGaa                                                                                      <  2:391990/129‑1 (MQ=255)               ccgccAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGaa                                                                                      >  1:391990/1‑129 (MQ=255)                cgccAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAgg                                                                                 <  2:484890/133‑1 (MQ=255)                  ccAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCt                                                                               >  1:76329/1‑133 (MQ=255)                      ggCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACa                                                                                    <  2:495725/124‑1 (MQ=255)                      ggCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACa                                                                                    >  1:495725/1‑124 (MQ=255)                             aCTTTGAACTCCACGTTAGTCTCTCCGTCTTACTTGACTTGACCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGc                                                                    <  2:137095/133‑1 (MQ=255)                                  gAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTGGTGCATCAGCAATAc                                                               >  1:171508/1‑133 (MQ=255)                                                    tcCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCa                                             >  2:422605/1‑133 (MQ=255)                                                    tcCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCa                                             >  2:300723/1‑133 (MQ=255)                                                       gtgtTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAgg                                          <  1:82213/133‑1 (MQ=255)                                                             aTTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATaa                                    <  2:107290/133‑1 (MQ=255)                                                                    tGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACAt                             >  1:340515/1‑133 (MQ=255)                                                                    tGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACAt                             >  2:105904/1‑133 (MQ=255)                                                                      cccAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCaaaaaa                                   <  2:81593/125‑1 (MQ=255)                                                                      cccAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATaaa                                   >  1:81593/1‑125 (MQ=255)                                                                             ttCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATcc                    <  1:105904/133‑1 (MQ=255)                                                                               cAGGACCTGTGGAGGCAAGCTTGGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCca                  <  1:170051/133‑1 (MQ=255)                                                                                 ggACCTCTTGAGGCACGCTTAGCGCCTCTGCTTTATTTCTCCGTCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCcaca                <  1:167258/133‑1 (MQ=255)                                                                                   aCCTCTGGTGGCAAGTTTTGCGCCTCTGTTCTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAAc              <  1:134635/133‑1 (MQ=255)                                                                                   aCCTCTGGAGTCAAGTTTAGCGCCACTGGTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAAc              <  2:171508/133‑1 (MQ=255)                                                                                    ccTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCTAGGAAATAAACCACATGACACATCCCACAAcc             >  2:134635/1‑133 (MQ=255)                                                                                    ccTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAAcc             <  1:300723/133‑1 (MQ=255)                                                                                        tGGAGGCAAGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCGAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAGTAAACCACATATCACATCCCATAACCAGAt         >  1:212150/1‑133 (MQ=255)                                                                                               aaGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAGTAAACCACATAACACATCCCACAACCAGATTGATAAt  >  1:71884/1‑133 (MQ=255)                                                                                               aaGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAGTAAACCACATAACACATCCCACAACCAGATTGATAAt  <  2:71884/133‑1 (MQ=255)                                                                                               aaGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAAt  >  2:263487/1‑133 (MQ=255)                                                                                               aaGCTTAGCGCCTCTGTTTTATTTCTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAAt  >  2:313420/1‑133 (MQ=255)                                                                                                                       |                                                                                                             AATACCCATTAATGCCGCCAGGGCGATAACTTTGAACTTCACGTTAGTCTCTCCGTGTTAATTGACTTGCCCAAAATTCAGGACCTCTGGAGGCAAGCTTAGCGCCTCTGTTTTATTTTTCCATCAGATAGCGCCTAACTGAACAAGGCTTGTGCATGAGCAATACCGTCTCTCGCCAGCTGCAAGGCAATAAACCACATAACACATCCCACAACCAGATTGATAAT  >  NC_000913/3068039‑3068265

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