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%	96.7 / NA	31	*118* (TAG→TAA)	ybaA	DUF1428 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/14);  total (17/14)

CTTCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAGGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCA  >  NC_000913/476892‑477146                                                                                                                                      |                                                                                                                          cTTCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGc                                                                                                                      >  2:545929/1‑139 (MQ=255)               cgcgACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTgc                                                                                                        <  1:464391/139‑1 (MQ=255)                                    atgatgTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCaaaaa                                                                                   <  2:87855/139‑1 (MQ=255)                                        tggtggAACCACGGAAGAAAGGGTTCTGCGAGTTCAAGGCGTTTGACCGCAAACCAAATATCTAAGGGCGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTca                                                                               <  1:1016648/135‑1 (MQ=255)                                        tgTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTgc                                                                                                        >  2:396593/1‑114 (MQ=255)                                        tgTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTgc                                                                                                        <  1:396593/114‑1 (MQ=255)                                                   cGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                    <  2:380019/139‑1 (MQ=255)                                                       tGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGc                                                                <  2:1015332/139‑1 (MQ=255)                                                       tGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGc                                                                >  1:878170/1‑139 (MQ=255)                                                        gAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTca                                                                               >  1:258429/1‑123 (MQ=255)                                                        gAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTca                                                                               <  2:258429/123‑1 (MQ=255)                                                        gAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                               <  1:68734/139‑1 (MQ=255)                                                        gAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                               <  2:878170/139‑1 (MQ=255)                                                                tCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTaa                                                       >  2:479365/1‑139 (MQ=255)                                                                ccGGCGAGTCCATGCCGTTTCACGGCAAGCGAATGATCTATGGCGGCTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCGGTGCGCCGTGCCGTCAGCCCAAAAATCCACAAATGTTGCTGGGCGCAGCCTaa                                                       >  2:567670/2‑139 (MQ=255)                                                                    cGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGaaa                                                   >  1:1009623/1‑139 (MQ=255)                                                                         ccATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAtcttc                                              >  1:102346/1‑139 (MQ=255)                                                                                 tttGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTa                                      >  2:974995/1‑139 (MQ=255)                                                                                      cGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGcc                                 <  1:473315/139‑1 (MQ=255)                                                                                         cAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGGCGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCaaaa                                                                                    <  1:526088/85‑1 (MQ=255)                                                                                         cAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCaaaa                                                                                    >  2:526088/1‑85 (MQ=255)                                                                                           aGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCt                            >  1:347522/1‑139 (MQ=255)                                                                                                      tATGGCGGATTCGAGTCAATCATCGACGAATAAGTATCGTGACGGGCTGTGCGGCGTGCCGGCAGCCCAAAAATTCAGAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCCTGCTGTAGagccagc                 >  2:583075/1‑139 (MQ=255)                                                                                                            ggATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGAc           >  2:226158/1‑139 (MQ=255)                                                                                                             gATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGagc                     >  2:49209/1‑128 (MQ=255)                                                                                                             gATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGagc                     <  1:49209/128‑1 (MQ=255)                                                                                                              aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCg         >  1:65942/1‑139 (MQ=255)                                                                                                               ttCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGt        >  2:317581/1‑139 (MQ=255)                                                                                                                    gTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGc   <  1:982384/139‑1 (MQ=255)                                                                                                                    gTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGc   >  1:474901/1‑139 (MQ=255)                                                                                                                     tCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCa  <  2:65942/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                          CTTCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAGGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCA  >  NC_000913/476892‑477146

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