Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	477,025	G→A	*118* (TAG→TAA)	ybaA →	DUF1428 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	477,025	0	G	A	100.0%	124.7 / NA	39	*118* (TAG→TAA)	ybaA	DUF1428 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/21);  total (18/21)

CCCTTCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAGGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACG  >  NC_000913/476890‑477148                                                                                                                                        |                                                                                                                            cccTTCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTa                                                                                                                          <  1:467924/139‑1 (MQ=255)    ttCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCg                                                                                                                       >  2:1106016/1‑139 (MQ=255)      cAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTg                                                                                                                     <  2:565522/139‑1 (MQ=255)              gTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCtg                                                                                                             >  2:497141/1‑139 (MQ=255)                        ctgctAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGc                                                                                                   >  2:1655215/1‑139 (MQ=255)                                        gatgTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAAtt                                                                                   >  2:340009/1‑139 (MQ=255)                                                    aCGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCtgtg                                                                                                           >  2:486960/1‑103 (MQ=255)                                                    aCGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCtgtg                                                                                                           <  1:486960/103‑1 (MQ=255)                                                     cGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTcac                                                                                <  2:359653/129‑1 (MQ=255)                                                     cGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTcac                                                                                >  1:359653/1‑129 (MQ=255)                                                     cGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                      <  1:1655215/139‑1 (MQ=255)                                                     cGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                      <  2:1138565/139‑1 (MQ=255)                                                                gTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGAGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGt                                                           <  1:1106016/139‑1 (MQ=255)                                                                      cGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTgc                                                                                                          >  2:97561/1‑86 (MQ=255)                                                                      cGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTgc                                                                                                          <  1:97561/86‑1 (MQ=255)                                                                                    ttGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTaa                                       >  2:835527/1‑139 (MQ=255)                                                                                               cGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGt                            <  1:340009/139‑1 (MQ=255)                                                                                                     aTCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGagcc                      <  1:835527/139‑1 (MQ=255)                                                                                                         aTGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATg                                                                           >  1:576583/1‑82 (MQ=255)                                                                                                         aTGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATg                                                                           <  2:576583/82‑1 (MQ=255)                                                                                                               gATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAt                                                    <  2:199614/99‑1 (MQ=255)                                                                                                               gATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAt                                                    >  1:199614/1‑99 (MQ=255)                                                                                                               gATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGAcc            >  1:1421534/1‑139 (MQ=255)                                                                                                                aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCaaaaa                                                                                     >  2:143281/1‑65 (MQ=255)                                                                                                                aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCaaaaa                                                                                     <  1:143281/65‑1 (MQ=255)                                                                                                                aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCt                                  >  2:1331504/1‑116 (MQ=255)                                                                                                                aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCt                                  <  1:1331504/116‑1 (MQ=255)                                                                                                                  tCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                 <  1:310370/83‑1 (MQ=255)                                                                                                                  tCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                 >  2:310370/1‑83 (MQ=255)                                                                                                                         aaTCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACg  <  1:497141/139‑1 (MQ=255)                                                                                                                               cgacgaATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                 <  2:1578319/70‑1 (MQ=255)                                                                                                                               cgacgaATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                 >  1:1578319/1‑70 (MQ=255)                                                                                                                                 acgaATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATg                                                                           >  1:1031648/1‑58 (MQ=255)                                                                                                                                 acgaATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATg                                                                           <  2:1031648/58‑1 (MQ=255)                                                                                                                                 acgaATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGt       <  2:765971/126‑1 (MQ=255)                                                                                                                                 acgaATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGt       >  1:765971/1‑126 (MQ=255)                                                                                                                                     aTAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                 <  2:270224/64‑1 (MQ=255)                                                                                                                                     aTAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                 >  1:270224/1‑64 (MQ=255)                                                                                                                                     aTAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGccc                 <  2:1572144/112‑1 (MQ=255)                                                                                                                                        |                                                                                                                            CCCTTCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAGGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACG  >  NC_000913/476890‑477148

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