Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,051,617	C→T	P91S (CCT→TCT)	yihI →	Der GTPase‑activating protein YihI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,051,617	0	C	T	100.0%	86.7 / NA	26	P91S (CCT→TCT)	yihI	Der GTPase‑activating protein YihI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/13);  total (13/13)

AGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGCAGGCGAAACCCTGAGTGCCGAAGAGCAATCCTGGGTGGATGCTAAACTG  >  NC_000913/4051484‑4051745                                                                                                                                      |                                                                                                                                 aGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGc                                                                                                                            >  1:404659/1‑140 (MQ=255)     ggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTc                                                                                                                        >  2:53370/1‑140 (MQ=255)             caccacGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAgg                                                                                                                <  1:324990/140‑1 (MQ=255)             caccacGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAgg                                                                                                                <  1:448582/140‑1 (MQ=255)                   gTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGt                                                                                                          <  2:404659/140‑1 (MQ=255)                    tCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGtt                                                                                                         >  1:36963/1‑140 (MQ=255)                          agcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAAtt                                                                                                   >  1:359186/1‑140 (MQ=255)                                      aaCGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGa                                                                                       >  2:255627/1‑140 (MQ=255)                                       aCGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGAt                                                                                      <  2:359186/140‑1 (MQ=255)                                        cGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATg                                                                                     <  1:344082/140‑1 (MQ=255)                                           aCCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGc                                                                                  >  2:414000/1‑140 (MQ=255)                                                                          ccTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGc                                                   <  2:293684/140‑1 (MQ=255)                                                                                   ttGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGCAGGCg                                              >  2:349429/1‑136 (MQ=255)                                                                                   ttGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGCAGGCg                                              <  1:349429/136‑1 (MQ=255)                                                                                   ttGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGCAGGCGAAAc                                          >  2:153238/1‑140 (MQ=255)                                                                                    tGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGCAGGCGAAAcc                                         <  2:36963/140‑1 (MQ=255)                                                                                        cGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGCAGGCGAAACCCTGa                                     <  1:53370/140‑1 (MQ=255)                                                                                        cGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGCAGGCGAAACCCTGa                                     >  1:116591/1‑140 (MQ=255)                                                                                          tGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGCAGGCg                                              >  1:140787/1‑129 (MQ=255)                                                                                          tGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGCAGGCg                                              <  2:140787/129‑1 (MQ=255)                                                                                           gACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGCAGGCGAAACCCTGAGTg                                  >  1:34498/1‑140 (MQ=255)                                                                                                   aaGTCACCAAACAGCACAAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGCAGGCGAAACCCTGAGTGCCGAagag                          <  2:34498/140‑1 (MQ=255)                                                                                                                   caAACCGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGCAGGCGAAACCCTGAGTGCCGAAGAGCAATCCTGGGTGGATg          <  1:383562/140‑1 (MQ=255)                                                                                                                        cGAAGAGTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGCAGGCGAAACCCTGAGTGCCGAAGAGCAATCCTGGGTGGATGCTaaa     <  2:477583/140‑1 (MQ=255)                                                                                                                           agagTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGTAGGCGAAACCCTGAGTGCCGAAGAGCAATCCTGGGTGGATGCTAAACTg  >  2:7327/1‑140 (MQ=255)                                                                                                                           agagTGAGAAATCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGCAGGCGAAACCCTGAGTGCCGAAGAGCAATCCTGGGTGGATGCTAAACTg  >  1:294542/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                 AGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACGGATGAGCGTCTGGATGCGCTGCTGGAACGTCTGGAAGCAGGCGAAACCCTGAGTGCCGAAGAGCAATCCTGGGTGGATGCTAAACTG  >  NC_000913/4051484‑4051745

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