Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,993,856	G→A	*73* (TAG→TAA)	ygeI →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,993,856	0	G	A	100.0%	104.3 / NA	33	*73* (TAG→TAA)	ygeI	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/17);  total (16/17)

AAACACATCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATAT  >  NC_000913/2993722‑2993987                                                                                                                                       |                                                                                                                                    aaaCACATCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGa                                                                                                                                 <  1:724377/139‑1 (MQ=255)      caTCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGTCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTCACTATAAAATTGAGTGGTAAACGATTGAc                                                                                                                            <  1:1157770/139‑1 (MQ=255)      caTCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGAc                                                                                                                            <  2:1708765/139‑1 (MQ=255)         cAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATGTGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGAACaa                                                                                                                         >  2:1244028/1‑139 (MQ=255)                                    aaTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTataa                                                                                              >  2:1052260/1‑139 (MQ=255)                                      tACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTataata                                                                                            <  2:339150/139‑1 (MQ=255)                                       aCAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAg                                                                                           <  1:1244028/139‑1 (MQ=255)                                          aGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAt                                                                                        >  1:1424932/1‑139 (MQ=255)                                                       tCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACttt                                                                                                   >  2:567091/1‑115 (MQ=255)                                                       tCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACttt                                                                                                   <  1:567091/115‑1 (MQ=255)                                                       tCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCtt                                                                              <  1:1368105/136‑1 (MQ=255)                                                       tCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCtt                                                                              >  2:1368105/1‑136 (MQ=255)                                                               gatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                   >  2:1282485/1‑139 (MQ=255)                                                                   gcgcCTTTTCAAAATGATATGTATTCAGCAGATATTCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcac                                                               >  2:1121089/1‑139 (MQ=255)                                                                                                  aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTa                                >  2:968693/1‑139 (MQ=255)                                                                                                  aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTa                                >  2:1329971/1‑139 (MQ=255)                                                                                                  aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTa                                <  1:1282485/139‑1 (MQ=255)                                                                                                    aaTCCTATCTCACCTTACTATAAAATTGAGTGGTAAACg                                                                                                                                  <  2:1220536/39‑1 (MQ=255)                                                                                                    aaTCCTATCTCACCTTACTATAAAATTGAGTGGTAAACg                                                                                                                                  >  1:1220536/1‑39 (MQ=255)                                                                                                     aTCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAAt                             <  1:968693/139‑1 (MQ=255)                                                                                                         tATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATgcg                                              <  2:669464/118‑1 (MQ=255)                                                                                                         tATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATgcg                                              >  1:669464/1‑118 (MQ=255)                                                                                                           tctcACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCaa                                                                      >  1:276192/1‑92 (MQ=255)                                                                                                           tctcACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCaa                                                                      <  2:276192/92‑1 (MQ=255)                                                                                                           tctcACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTaa                       <  2:137569/139‑1 (MQ=255)                                                                                                           tctcACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTaa                       >  1:137569/1‑139 (MQ=255)                                                                                                                ccTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                   <  2:727242/74‑1 (MQ=255)                                                                                                                ccTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                   >  1:727242/1‑74 (MQ=255)                                                                                                                ccTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTaa                       >  2:1540272/1‑134 (MQ=255)                                                                                                                ccTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTaa                       <  1:1540272/134‑1 (MQ=255)                                                                                                                   tACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTc               <  2:1424932/139‑1 (MQ=255)                                                                                                                             ttGAGTGGTAAGCGATTGACCAAATGAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTTAATTCTTAACCATCTTCACGTTTTGta     <  1:1121089/139‑1 (MQ=255)                                                                                                                                aGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTatat  >  2:1318068/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                    AAACACATCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATAT  >  NC_000913/2993722‑2993987

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