Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,182,227	T→C	T70A (ACG→GCG)	ygiD ←	4,5‑DOPA dioxygenase extradiol

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,182,227	0	T	C	100.0%	90.9 / NA	26	T70A (ACG→GCG)	ygiD	4,5‑DOPA dioxygenase extradiol
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/10);  total (16/10)

CTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGGGTATACAAATTATCTTCC  >  NC_000913/3182092‑3182348                                                                                                                                        |                                                                                                                          cttTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTcc                                                                                                                       >  1:473973/1‑140 (MQ=255)  tttATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCa                                                                                                                      >  2:195847/1‑140 (MQ=255)           cGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGcct                                                                                                                         <  1:520087/128‑1 (MQ=255)           cGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGcct                                                                                                                         >  2:520087/1‑128 (MQ=255)           cGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATc                                                                                                                    <  1:327118/133‑1 (MQ=255)           cGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATc                                                                                                                    >  2:327118/1‑133 (MQ=255)             tCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACt                                                                                                           <  1:195847/140‑1 (MQ=255)                                 aGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCa                                                                                       <  2:63335/140‑1 (MQ=255)                                             cGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAac                                                                           >  1:482951/1‑140 (MQ=255)                                                  tgCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAAccacca                                                                      >  2:276657/1‑140 (MQ=255)                                                  tgCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAAccacca                                                                      >  2:296354/1‑140 (MQ=255)                                                  tgCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAAccacca                                                                      >  2:498968/1‑140 (MQ=255)                                                          cAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCt                                                              >  1:243200/1‑140 (MQ=255)                                                                       aGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCa                                                 >  2:227247/1‑140 (MQ=255)                                                                       aGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCATGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCa                                                 >  2:281306/1‑140 (MQ=255)                                                                            gATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTc                                            <  1:276657/140‑1 (MQ=255)                                                                                 tGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATccc                                       >  1:269665/1‑140 (MQ=255)                                                                                      aTCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACt                                  >  1:248284/1‑140 (MQ=255)                                                                                          tACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTg                              >  1:7034/1‑140 (MQ=255)                                                                                          tACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTg                              <  1:498968/140‑1 (MQ=255)                                                                                          tACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTg                              <  1:296354/140‑1 (MQ=255)                                                                                                 ccTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTg                       >  1:474181/1‑140 (MQ=255)                                                                                                 ccTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTg                       <  2:8619/140‑1 (MQ=255)                                                                                                    gCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGc                      <  2:2969/138‑1 (MQ=255)                                                                                                    gCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGc                      >  1:2969/1‑138 (MQ=255)                                                                                                              accaAAGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGGGTATACAAAt          <  2:7034/140‑1 (MQ=255)                                                                                                                   aGTCATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGGGTATACAAATTATCt     <  2:269665/140‑1 (MQ=255)                                                                                                                      cATGAATCGTGGGCGGCGCCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGGGTATACAAATTATCTTcc  >  1:448654/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                          CTTTATCCAGCGTCACCGGGATCGGCGCTAACAGCTCAACCAGACGCTGTGCCAGCGCAGGCGAACCCGGAGCAGGATAATGCGTATCGTACAGCGCCTGCGGGAAGCCACCAAAGTCATGAATCGTGGGCGGCGTCTCCATCGCGGTCACTCCTGTTCCACGGGTAAACCAGTGAGCCGAAACCACCACAATCGCTTGCGGGCGTGGCAATGTCATCCCCAACTTCTGCCAGCTGCGGGTATACAAATTATCTTCC  >  NC_000913/3182092‑3182348

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