Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,051,495	A→G	N50S (AAC→AGC)	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,495	0	A	G	100.0%	155.8 / NA	45	N50S (AAC→AGC)	yihI	Der GTPase‑activating protein YihI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/23);  total (22/23)

AACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTC  >  NC_000913/4051365‑4051627                                                                                                                                   |                                                                                                                                     aaCTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTc                                                                                                                             >  2:90051/1‑140 (MQ=255)  aCTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCa                                                                                                                            <  2:793574/140‑1 (MQ=255)      aCGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGgca                                                                                                                        >  1:522299/1‑140 (MQ=255)        gcagcaAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGgcagc                                                                                                                      >  2:581866/1‑140 (MQ=255)         cagcaAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGgcagca                                                                                                                     <  1:90051/140‑1 (MQ=255)              aaGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGt                                                                                                                              >  1:257314/1‑126 (MQ=255)              aaGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGt                                                                                                                              <  2:257314/126‑1 (MQ=255)                   cATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGaa                                                                                                           <  2:10474/140‑1 (MQ=255)                   cATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGaa                                                                                                           <  2:356898/140‑1 (MQ=255)                     tGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACg                                                                                                         >  1:590724/1‑140 (MQ=255)                      gCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAgg                                                                                                                          <  2:491185/122‑1 (MQ=255)                      gCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAgg                                                                                                                          >  1:491185/1‑122 (MQ=255)                           aGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCaa                                                                                                   >  1:636347/1‑140 (MQ=255)                           aGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCaa                                                                                                   >  2:894007/1‑140 (MQ=255)                           aGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCaa                                                                                                   >  2:26894/1‑140 (MQ=255)                           aGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCaa                                                                                                   >  2:258223/1‑140 (MQ=255)                            gcgcGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCaaa                                                                                                  <  1:894007/140‑1 (MQ=255)                                 tCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATc                                                                                             <  1:699243/140‑1 (MQ=255)                                 tCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATc                                                                                             <  1:811903/140‑1 (MQ=255)                                   gAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCa                                                                                           <  2:805478/140‑1 (MQ=255)                                                  agCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTaaa                                                                            <  1:851686/140‑1 (MQ=255)                                                      ggATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAAccc                                                                        <  1:258223/140‑1 (MQ=255)                                                         tCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTa                                                                     <  2:522299/140‑1 (MQ=255)                                                          cAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTAt                                                                    <  2:636347/140‑1 (MQ=255)                                                            ggAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTc                                                                  >  1:400166/1‑140 (MQ=255)                                                            ggAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTc                                                                  >  2:416346/1‑140 (MQ=255)                                                             gAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAAc                                                                          >  1:15470/1‑131 (MQ=255)                                                             gAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAAc                                                                          <  2:15470/131‑1 (MQ=255)                                                              aaGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTAt                                                                    <  1:688969/136‑1 (MQ=255)                                                              aaGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTAt                                                                    >  2:688969/1‑136 (MQ=255)                                                              aaGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCa                                                                >  1:121037/1‑140 (MQ=255)                                                                 cTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTg                                                             >  1:315454/1‑140 (MQ=255)                                                                       aTCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTg                                                       >  2:629745/1‑140 (MQ=255)                                                                           cAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTg                                                   <  1:629745/140‑1 (MQ=255)                                                                            aaGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGa                                                  >  1:728249/1‑140 (MQ=255)                                                                                    gAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCa                                          <  1:652405/140‑1 (MQ=255)                                                                                         aaCGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCaaa                                     <  1:733525/140‑1 (MQ=255)                                                                                         aaCGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCaaa                                     <  2:121037/140‑1 (MQ=255)                                                                                               gtGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGcac                               <  1:416346/140‑1 (MQ=255)                                                                                                       gcgcCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGaa                       <  1:581866/140‑1 (MQ=255)                                                                                                       gcgcCGGGCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGaa                       >  2:950073/1‑140 (MQ=255)                                                                                                             ggCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTga                 >  1:418376/1‑140 (MQ=255)                                                                                                              gCAGCCGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCTAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTgag                <  2:100298/140‑1 (MQ=255)                                                                                                                  ccGCGCAGCGGGCGGAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAc            >  1:609097/1‑140 (MQ=255)                                                                                                                          cgggcggAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCtt    <  2:400166/140‑1 (MQ=255)                                                                                                                            ggcggAAGCACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTc  >  1:470569/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                     AACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTC  >  NC_000913/4051365‑4051627

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