Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,045,596	A→G	E466E (GAA→GAG)	yihF →	uncharacterized protein YihF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,045,596	0	A	G	100.0%	171.8 / NA	48	E466E (GAA→GAG)	yihF	uncharacterized protein YihF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/21);  total (27/21)

ATAAACAAATCAGCGGGATGATGACATTAGGTCAGATGTTTCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGC  >  NC_000913/4045460‑4045727                                                                                                                                         |                                                                                                                                    aTAAACAAATCAGCGGGATGATGACATTAGGTCAGATGTTTCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGttt                                                                                                                                  >  2:575328/1‑140 (MQ=255) aTAAACAAATCAGCGGGATGATGACATTAGGTCAGATGTTTCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGttt                                                                                                                                  <  1:370925/140‑1 (MQ=255)               gggATGATGACATTAGGTCAGATGTTTCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCgg                                                                                                                    <  2:143433/140‑1 (MQ=255)                   tgatgaCATTAGGTCAGATGTTTCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGt                                                                                                                >  2:464487/1‑140 (MQ=255)                         aaTTAGGTCAGATGTTTCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGtt                                                                                                          <  1:16244/139‑1 (MQ=255)                           ttAGGTCAGATGTTTCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCa                                                                                                        <  1:336702/140‑1 (MQ=255)                           ttAGGTCAGATGTTTCACTTAACCACGATTGTCAATAATTCCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCa                                                                                                        <  2:234207/140‑1 (MQ=255)                                  aGATGTTTCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAAt                                                                                                 <  2:134055/140‑1 (MQ=255)                                   gATGTTTCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATa                                                                                                >  2:45582/1‑140 (MQ=255)                                    aTGTTTCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATac                                                                                               >  1:548545/1‑140 (MQ=255)                                    aTGTTTCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATac                                                                                               >  2:152119/1‑140 (MQ=255)                                    aTGTTTCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATac                                                                                               >  1:618170/1‑140 (MQ=255)                                        ttCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATactcac                                                                                           >  2:84597/1‑140 (MQ=255)                                         tCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTg                                                                                                            >  2:7649/1‑122 (MQ=255)                                         tCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTg                                                                                                            <  1:7649/122‑1 (MQ=255)                                         tCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATactcact                                                                                          <  1:524630/140‑1 (MQ=255)                                         tCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATactcact                                                                                          <  2:364815/140‑1 (MQ=255)                                             ttAATCACGATTGACAACAATCCCGCCTCGCTGCAGCTGCGTTATACACCTGGTAAAGTTGTTTTTAACAGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGcc                                                                                      <  1:382383/140‑1 (MQ=255)                                             ttAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGcc                                                                                      <  2:618170/140‑1 (MQ=255)                                               aaTCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCctc                                                                                    <  1:31824/140‑1 (MQ=255)                                               aaTCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCctc                                                                                    <  1:83317/140‑1 (MQ=255)                                                  cACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGc                                                                                 <  1:329773/140‑1 (MQ=255)                                                                    cGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTaa                                                               <  2:28765/140‑1 (MQ=255)                                                                      ccTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTaaaa                                                             <  2:624009/140‑1 (MQ=255)                                                                      ccTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTaaaa                                                             >  2:588635/1‑140 (MQ=255)                                                                         cGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAactt                                                          >  2:612265/1‑140 (MQ=255)                                                                         cGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTAGTGCCGGACGTTTTGTCCATTAAAATACTCACTCGCCTCTGCACATTTAGAGGCGAGTAAAactt                                                          >  2:232961/1‑140 (MQ=255)                                                                            tGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAacttact                                                       <  1:464487/140‑1 (MQ=255)                                                                              cAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCTTTCAGAGGCAAGTAAAActtactta                                                     >  1:217654/1‑140 (MQ=255)                                                                                gCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACCCGCGTGTCCTCATTCAGAGGCGAGTAAAACGTACTTAtt                                                   >  2:324004/1‑140 (MQ=255)                                                                                           acaggggTAAAGTTGTTTTTAACGGTCAGGAGATGAGAGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCtgtg                                        <  1:263133/136‑1 (MQ=255)                                                                                              cGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGtgatg                                     >  2:428124/1‑140 (MQ=255)                                                                                              cGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGtgatg                                     >  1:363508/1‑140 (MQ=255)                                                                                              cGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGtgatg                                     >  1:522317/1‑140 (MQ=255)                                                                                              cGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGtgatg                                     >  1:641280/1‑140 (MQ=255)                                                                                              cGGGCAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATACAGAGGCGAGTAAAACTTACTGATTCTGACGCTCtgatg                                     >  2:169935/1‑140 (MQ=255)                                                                                                ggTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGAc                                   >  1:638457/1‑140 (MQ=255)                                                                                                          tttttAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAAAGGCGAGTAAAATTTATTTATTCTGTCGCTGTGATGACATgagactgt                         >  2:263133/1‑136 (MQ=255)                                                                                                            tttAACGGACAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTc                       >  2:350775/1‑140 (MQ=255)                                                                                                                    aCAGGAGATGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTAGTCTGACGGTGTGATGACATGAGAGACGTCAGAAGCCg               >  1:256676/1‑140 (MQ=255)                                                                                                                        gagaTGAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCt           >  2:372088/1‑140 (MQ=255)                                                                                                                             gAGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCAGTGATGAGATGGGATACGTCAGAAGCCGGTCTTTTTc      >  1:170875/1‑140 (MQ=255)                                                                                                                              aGCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTcc     >  1:2705/1‑140 (MQ=255)                                                                                                                               gCGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCt    <  1:152119/140‑1 (MQ=255)                                                                                                                                cGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCCGAAGCCGGTCTTTTTCCTg   >  2:415521/1‑140 (MQ=255)                                                                                                                                cGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTg   >  2:334445/1‑140 (MQ=255)                                                                                                                                cGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTg   <  1:612265/140‑1 (MQ=255)                                                                                                                                cGAAGAAGAGTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTg   <  1:334445/140‑1 (MQ=255)                                                                                                                                 gaagaagagTTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGc  <  2:548545/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                    ATAAACAAATCAGCGGGATGATGACATTAGGTCAGATGTTTCAGTTAATCACGATTGACAACAATACCGCCTCGCTGCAGCTGCGTTATACACCGGGTAAAGTTGTTTTTAACGGACAGGAGATGAGCGAAGAAGAATTTATGTCTCGTGCCGGACGTTTTGTTCATTAAAATACTCACTCGCCTCTGCTCATTCAGAGGCGAGTAAAACTTACTTATTCTGACGCTGTGATGACATGAGAGACGTCAGAAGCCGGTCTTTTTCCTGC  >  NC_000913/4045460‑4045727

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