Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	775,120	G→A	A123A (GCG→GCA)	ybgC →	acyl‑CoA thioester hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	775,120	0	G	A	100.0%	118.8 / NA	37	A123A (GCG→GCA)	ybgC	acyl‑CoA thioester hydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/20);  total (17/20)

CAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCGCTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGC  >  NC_000913/775002‑775250                                                                                                                       |                                                                                                                                   cAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgc                                                                                                                <  1:152156/139‑1 (MQ=255)  aaTGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgcg                                                                                                               >  1:362104/1‑139 (MQ=255)  aaTGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgcg                                                                                                               >  1:668895/1‑139 (MQ=255)   aTGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCgg                                                                                                              <  2:362104/139‑1 (MQ=255)    tGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGa                                                                                                             >  1:33845/1‑139 (MQ=255)       gTGGCACCTCTTTGGTTTTCAAGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGATCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGtt                                                                                                          <  2:449797/139‑1 (MQ=255)         ggCACCTCTTTGGTTGTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTa                                                                                                        <  2:746002/139‑1 (MQ=255)          gCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTaa                                                                                                       >  1:172113/1‑139 (MQ=255)                    ggTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtgactg                                                                                             <  2:172113/139‑1 (MQ=255)                      ttttCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtgactgac                                                                                           <  2:668895/139‑1 (MQ=255)                        ttCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtgact                                                                                              >  2:763606/1‑134 (MQ=255)                        ttCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtgact                                                                                              <  1:763606/134‑1 (MQ=255)                                           aCGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGc                                                                                                     >  1:136028/1‑108 (MQ=255)                                           aCGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGc                                                                                                     <  2:136028/108‑1 (MQ=255)                                           aCGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTg                                                                         <  1:569181/136‑1 (MQ=255)                                           aCGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTg                                                                         >  2:569181/1‑136 (MQ=255)                                           aCGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTc                                                                      >  2:643062/1‑139 (MQ=255)                                            cGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTcc                                                                     >  1:522764/1‑139 (MQ=255)                                             gCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCt                                                                    <  2:50250/139‑1 (MQ=255)                                                          tGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGaa                                                                                      >  1:74235/1‑108 (MQ=255)                                                          tGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGaa                                                                                      <  2:74235/108‑1 (MQ=255)                                                                atgaaGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgc                                                                                                                <  1:650647/76‑1 (MQ=255)                                                                atgaaGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgc                                                                                                                >  2:650647/1‑76 (MQ=255)                                                                 tgaaGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtga                                                                                                <  1:389690/91‑1 (MQ=255)                                                                 tgaaGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtga                                                                                                >  2:389690/1‑91 (MQ=255)                                                                              cTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGt                                                  <  2:218832/124‑1 (MQ=255)                                                                              cTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGt                                                  >  1:218832/1‑124 (MQ=255)                                                                              cTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGtt                                   >  1:510335/1‑139 (MQ=255)                                                                                     tttGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttgattttg                            >  2:274253/1‑139 (MQ=255)                                                                                                      aaaaTGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGc           <  2:33845/139‑1 (MQ=255)                                                                                                        aaTGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttgattttg                            >  1:213462/1‑120 (MQ=255)                                                                                                        aaTGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttgattttg                            <  2:213462/120‑1 (MQ=255)                                                                                                          tGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCt       <  2:228131/139‑1 (MQ=255)                                                                                                            aaGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTg     <  2:510335/139‑1 (MQ=255)                                                                                                              gCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTggg   >  2:331747/1‑139 (MQ=255)                                                                                                               ccTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGc  <  1:331747/139‑1 (MQ=255)                                                                                                               ccTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGc  <  1:682751/139‑1 (MQ=255)                                                                                                                       |                                                                                                                                   CAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCGCTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGC  >  NC_000913/775002‑775250

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: ‑