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

CACATCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATAT  >  NC_000913/2993725‑2993987                                                                                                                                    |                                                                                                                                    cacaTCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTg                                                                                                                              >  2:345218/1‑139 (MQ=255)  acaTCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGa                                                                                                                             <  2:523573/139‑1 (MQ=255)     tCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCa                                                                                                                          >  1:645316/1‑139 (MQ=255)      cAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCaa                                                                                                                         >  2:549727/1‑139 (MQ=255)        acacCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCaaat                                                                                                                       <  2:139377/139‑1 (MQ=255)           ccAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCaaa                                                                                                                        >  1:67201/1‑135 (MQ=255)           ccAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCaaa                                                                                                                        <  2:67201/135‑1 (MQ=255)                      gATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACgaga                                                                                                         <  1:593588/139‑1 (MQ=255)                      gATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACgaga                                                                                                         <  1:241410/139‑1 (MQ=255)                       aTACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        <  1:730410/139‑1 (MQ=255)                         aCAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        <  1:670758/137‑1 (MQ=255)                         aCAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        >  2:670758/1‑137 (MQ=255)                         aCAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATAc                                                                                                      <  1:549727/139‑1 (MQ=255)                         aCAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATAc                                                                                                      <  2:249128/139‑1 (MQ=255)                             aTCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATAc                                                                                                      <  2:361280/135‑1 (MQ=255)                             aTCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATAc                                                                                                      >  1:361280/1‑135 (MQ=255)                             aTCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACtttt                                                                                                  >  2:214456/1‑139 (MQ=255)                              tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGTCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTta                                                                                                 >  1:535425/1‑139 (MQ=255)                              tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        <  1:398575/132‑1 (MQ=255)                              tCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        >  2:398575/1‑132 (MQ=255)                                     cAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGt                                                                                          >  1:571939/1‑139 (MQ=255)                                       aGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAt                                                                                        >  2:141374/1‑139 (MQ=255)                                          aTTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGAc                                                                                                                            <  2:530197/100‑1 (MQ=255)                                          aTTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGAc                                                                                                                            >  1:530197/1‑100 (MQ=255)                                            ttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        <  1:371434/118‑1 (MQ=255)                                            ttATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGAt                                                                                                        >  2:371434/1‑118 (MQ=255)                                             tataGTTTCAACATGGATGGTGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGtt                                                                                  >  2:675880/1‑139 (MQ=255)                                                    tCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttca                                                                           >  2:64171/1‑139 (MQ=255)                                                    tCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttca                                                                           >  2:353907/1‑139 (MQ=255)                                                     cAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcat                                                                          <  1:141374/139‑1 (MQ=255)                                                      aaCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATGAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcatt                                                                         >  2:726511/1‑139 (MQ=255)                                                      aaCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcatt                                                                         >  1:316683/1‑139 (MQ=255)                                                       aCATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcattc                                                                        <  2:645316/139‑1 (MQ=255)                                                              tggCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAATATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTc                                                                 >  1:631197/1‑139 (MQ=255)                                                                                       ttCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCAcc                                        >  2:521093/1‑139 (MQ=255)                                                                                        tCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCa                                       >  2:600164/1‑139 (MQ=255)                                                                                        tCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCa                                       >  1:494534/1‑139 (MQ=255)                                                                                        tCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCa                                       >  2:43866/1‑139 (MQ=255)                                                                                          agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaaa                                     >  2:802504/1‑139 (MQ=255)                                                                                            aagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCt                                                                  >  1:765257/1‑108 (MQ=255)                                                                                            aagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCt                                                                  <  2:765257/108‑1 (MQ=255)                                                                                            aagaGAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATTAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAg                                   <  1:226325/139‑1 (MQ=255)                                                                                               aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTa                                >  1:293105/1‑139 (MQ=255)                                                                                                  aTCCTATCTCACCTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAAt                             <  2:316683/139‑1 (MQ=255)                                                                                                              cTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCt                 <  2:494534/139‑1 (MQ=255)                                                                                                              cTTACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCt                 >  2:307910/1‑139 (MQ=255)                                                                                                                tACTATAAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTc               <  1:802504/139‑1 (MQ=255)                                                                                                                    ataAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaa                                                        <  2:413131/94‑1 (MQ=255)                                                                                                                    ataAAATTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaa                                                        >  1:413131/1‑94 (MQ=255)                                                                                                                      aaaaTTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcatt                                                                         >  2:354609/1‑75 (MQ=255)                                                                                                                      aaaaTTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcatt                                                                         <  1:354609/75‑1 (MQ=255)                                                                                                                      aaaaTTGAGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGttt         >  1:699034/1‑139 (MQ=255)                                                                                                                             aGTGGTAAACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAACATCTTCACGTTTTGTatat  >  2:378509/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                    CACATCAACACCAGTGCTACTGATACAGATCAAATACAAGCATTTATAGTTTCAACATGGATGGCGCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATAT  >  NC_000913/2993725‑2993987

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