Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	214,689	T→C	Y133Y (TAT→TAC)	yaeQ →	uncharacterized protein YaeQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	214,689	0	T	C	100.0%	138.4 / NA	40	Y133Y (TAT→TAC)	yaeQ	uncharacterized protein YaeQ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/18);  total (22/18)

GCCTGCACCCAGGCCGCAGAAGTGGCGCTGTTTACCTATAATAGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCC  >  NC_000913/214573‑214818                                                                                                                     |                                                                                                                                  gCCTGCACCCAGACCGCAGAAGTGGCGCTGTTTACCTATAATAGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGaacaac                                                                                                                   >  1:382836/1‑133 (MQ=255)  ccTGCACCCAGACCGCAGAAGTGGCGCTGTTTACCTATAATAGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACt                                                                                                                  >  1:479635/1‑133 (MQ=255)                 cAGAAGTGGCGCTGTTTACCTATAATAGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAgcgc                                                                                                   >  1:102236/1‑133 (MQ=255)                 cAGAAGTGGCGCTGTTTACCTATAATAGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATAAACCACTGGCGAAAGTAAgcgc                                                                                                   >  1:201441/1‑133 (MQ=255)                       tGGCGCTGTTTACCTATAATAGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGc                                                                                             >  1:114059/1‑133 (MQ=255)                        ggCGCTGTTTACCTATAATAGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTCCCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGcc                                                                                            <  2:124070/133‑1 (MQ=255)                           gcTGTTTACCTATAATAGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGAt                                                                                         >  1:507699/1‑133 (MQ=255)                                    cTATAATAGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCCAAAGTAAGCGCATTAGCCGATGGTACCATg                                                                                >  1:131780/1‑133 (MQ=255)                                       taataGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACg                                                                             >  2:20662/1‑133 (MQ=255)                                         ataGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCt                                                                           >  1:416261/1‑133 (MQ=255)                                              cGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAgg                                                                      >  2:294081/1‑133 (MQ=255)                                              cGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAgg                                                                      >  1:262093/1‑133 (MQ=255)                                               gggCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGACTTTGCCCATCGTAACATGACGCTGCAGGc                                                                     >  1:195699/1‑133 (MQ=255)                                                 gcggcgCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAGAGTAAGCGTCTTGGCCAATCGTACCATGACGCTGCAGGCaa                                                                   >  2:159049/1‑133 (MQ=255)                                                 gcggcgCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCaa                                                                   >  2:362990/1‑133 (MQ=255)                                                 gcggcgCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCaa                                                                   <  1:83593/133‑1 (MQ=255)                                                 gcggcgCAAATCTGGAGGCAGCAAGATCCGAGAAAATGTGTCCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCaa                                                                   <  2:195699/133‑1 (MQ=255)                                                  cggcgCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATGGTATCATGACGCTGCAGGCAAc                                                                  >  1:242740/1‑133 (MQ=255)                                                          aTCTGGTGGTAGCAACATCAGAGCCAGTGTGTTCAGTTTGCGAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAgg                                                          <  2:242740/133‑1 (MQ=255)                                                                    agcaAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGAt                                                >  2:474142/1‑133 (MQ=255)                                                                     gcaAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGTGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATc                                               <  1:294081/133‑1 (MQ=255)                                                                        aaaTCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTgg                                            <  1:432985/133‑1 (MQ=255)                                                                               agCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCgg                                     <  2:114059/133‑1 (MQ=255)                                                                               agCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCgg                                     <  1:20662/133‑1 (MQ=255)                                                                               agCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCgg                                     <  2:479635/133‑1 (MQ=255)                                                                               agCAAATGTGTGCAGCTTCCCACTCTTTCCCTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCgg                                     <  1:157018/133‑1 (MQ=255)                                                                                 cAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGgat                                   <  2:102236/133‑1 (MQ=255)                                                                                 cAAATGTGTGCAGTTTGCCAATCTCTCCTTCTGGTACCTGGCCGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGgat                                   <  2:201441/133‑1 (MQ=255)                                                                                   aaTGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGgatga                                 <  2:299263/133‑1 (MQ=255)                                                                                    aTGTGTGCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGgatgat                                <  2:114909/133‑1 (MQ=255)                                                                                         tgCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGgatga                                 <  2:305183/127‑1 (MQ=255)                                                                                         tgCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGgatga                                 >  1:305183/1‑127 (MQ=255)                                                                                         tgCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGaa                           >  1:447108/1‑133 (MQ=255)                                                                                         tgCAGTTTGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGaa                           <  2:382836/133‑1 (MQ=255)                                                                                         tgCAGTTTCCCAATCTTTTCGTTTGGTACCGGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGaa                           <  1:159049/133‑1 (MQ=255)                                                                                                tGCCAATCTTTCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTg                    >  2:321548/1‑133 (MQ=255)                                                                                                          tCCGTCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGTTTATCGGATGATAAGAATAATCTGGAAGTGAACt          >  2:320729/1‑133 (MQ=255)                                                                                                              tCTGGTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAAc      <  1:362990/133‑1 (MQ=255)                                                                                                                  gTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGcc  >  2:60980/1‑133 (MQ=255)                                                                                                                  gTACCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGcc  >  2:82092/1‑133 (MQ=255)                                                                                                                     |                                                                                                                                  GCCTGCACCCAGGCCGCAGAAGTGGCGCTGTTTACCTATAATAGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCC  >  NC_000913/214573‑214818

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: ‑