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%	166.6 / NA	51	*73* (TAG→TAA)	ygeI	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/24);  total (27/24)

AAACACATCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATAT  >  NC_000913/2993722‑2993987                                                                                                                                       |                                                                                                                                    aaaCACATCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGa                                                                                                                                 <  2:1509681/139‑1 (MQ=255)               cAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCaaataaat                                                                                                                   >  1:1649258/1‑139 (MQ=255)                 gTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGa                                                                                                                             >  2:1605500/1‑127 (MQ=255)                 gTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGa                                                                                                                             <  1:1605500/127‑1 (MQ=255)                 gTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCa                                                                                                                 >  1:166588/1‑139 (MQ=255)                  tGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCAt                                                                                                                >  2:780225/1‑139 (MQ=255)                      aCTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGAc                                                                                                                            >  1:11070/1‑123 (MQ=255)                      aCTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGAc                                                                                                                            <  2:11070/123‑1 (MQ=255)                         gATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACgaga                                                                                                         >  1:222668/1‑139 (MQ=255)                            aCAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATAc                                                                                                      <  2:1649258/139‑1 (MQ=255)                                     aTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACttt                                                                                                   <  2:689722/133‑1 (MQ=255)                                     aTACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACttt                                                                                                   >  1:689722/1‑133 (MQ=255)                                             aTTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtctc                                                                                     <  2:166588/139‑1 (MQ=255)                                               ttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa                                                                                         >  2:4184/1‑133 (MQ=255)                                               ttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa                                                                                         <  1:4184/133‑1 (MQ=255)                                                       tCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTACAAAATTGAGTGGTAAACGATTGAc                                                                                                                            <  1:1391981/90‑1 (MQ=255)                                                       tCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTAGCTCACCTTACTATAAAATTGAGTGGTAAACGATTGAc                                                                                                                            >  2:1391981/1‑90 (MQ=255)                                                         aaCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcatt                                                                         <  2:800013/139‑1 (MQ=255)                                                               gatggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                   >  2:1816251/1‑139 (MQ=255)                                                                   gcgcCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATAc                                                                                                      <  2:1227271/100‑1 (MQ=255)                                                                   gcgcCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATAc                                                                                                      >  1:1227271/1‑100 (MQ=255)                                                                   gcgcCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcac                                                               <  1:529192/139‑1 (MQ=255)                                                                     gcCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGa                                                                                                                             <  2:1955171/75‑1 (MQ=255)                                                                     gcCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGa                                                                                                                             >  1:1955171/1‑75 (MQ=255)                                                                         tttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGa                                                         >  1:399932/1‑139 (MQ=255)                                                                         tttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTAT‑AAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaa                                                        >  1:973344/1‑139 (MQ=255)                                                                                           tCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCa                                       <  2:399932/139‑1 (MQ=255)                                                                                             agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        <  1:1559642/72‑1 (MQ=255)                                                                                             agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        >  2:1559642/1‑72 (MQ=255)                                                                                             agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaaa                                     >  1:1505792/1‑139 (MQ=255)                                                                                                agaTAATCCTATCTCACCTTACTAT‑AAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGtt                                 <  2:973344/139‑1 (MQ=255)                                                                                                  aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAg                                                                                           >  2:430710/1‑80 (MQ=255)                                                                                                  aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAg                                                                                           <  1:430710/80‑1 (MQ=255)                                                                                                  aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                   <  2:1457013/88‑1 (MQ=255)                                                                                                  aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTa                                <  1:1816251/139‑1 (MQ=255)                                                                                                  aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTAACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                   >  1:1457013/1‑88 (MQ=255)                                                                                                    aaTCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        <  1:516742/65‑1 (MQ=255)                                                                                                    aaTCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        >  2:516742/1‑65 (MQ=255)                                                                                                    aaTCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTaaa                              >  2:1415051/1‑139 (MQ=255)                                                                                                    aaTCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTaaa                              >  2:437370/1‑139 (MQ=255)                                                                                                           tctcACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTaa                       >  1:1231938/1‑139 (MQ=255)                                                                                                              cACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCa                    >  1:2056223/1‑139 (MQ=255)                                                                                                                  ttACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCtt                <  1:437370/139‑1 (MQ=255)                                                                                                                        taAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtt          <  2:1505792/139‑1 (MQ=255)                                                                                                                        taAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtt          <  1:780225/139‑1 (MQ=255)                                                                                                                                aGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCt                                            <  2:1532949/97‑1 (MQ=255)                                                                                                                                aGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCt                                            >  1:1532949/1‑97 (MQ=255)                                                                                                                                aGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGt      <  1:793662/135‑1 (MQ=255)                                                                                                                                aGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGt      >  2:793662/1‑135 (MQ=255)                                                                                                                                aGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTatat  >  2:288544/1‑139 (MQ=255)                                                                                                                                aGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTatat  >  2:932952/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                    AAACACATCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATAT  >  NC_000913/2993722‑2993987

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