Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,503,125	G→A	*223* (TAG→TAA)	ydcL →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,503,125	0	G	A	100.0%	133.1 / NA	42	*223* (TAG→TAA)	ydcL	lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/21);  total (21/21)

GTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAA  >  NC_000913/1502991‑1503243                                                                                                                                       |                                                                                                                       gTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTcc                                                                                                                    >  2:109143/1‑139 (MQ=255)       aaaGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgc                                                                                                              <  1:109143/139‑1 (MQ=255)        aaGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                             <  1:538410/139‑1 (MQ=255)                   cAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGa                                                                                                  <  1:1071404/139‑1 (MQ=255)                                        aaTAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                            <  1:1373349/108‑1 (MQ=255)                                        aaTAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                            >  2:1373349/1‑108 (MQ=255)                                            aCGAAAGTATGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGt                                                                         <  1:931243/139‑1 (MQ=255)                                                  gTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATc                                                                                    >  2:7508/1‑122 (MQ=255)                                                  gTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATc                                                                                    <  1:7508/122‑1 (MQ=255)                                                  gTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAAcc                                                                   >  1:859515/1‑139 (MQ=255)                                                  gTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAAcc                                                                   >  2:1255032/1‑139 (MQ=255)                                                     cGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCa                                                                >  2:166542/1‑139 (MQ=255)                                                           tGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTc                                                          >  1:288361/1‑139 (MQ=255)                                                               tttCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                      <  2:1029537/139‑1 (MQ=255)                                                                ttCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCa                                                                <  2:1125519/128‑1 (MQ=255)                                                                ttCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCa                                                                >  1:1125519/1‑128 (MQ=255)                                                                             aaCAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGgc                                        >  2:394401/1‑139 (MQ=255)                                                                               cAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGgcgc                                      <  2:859515/139‑1 (MQ=255)                                                                                  gTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAc                                                                                                                >  1:999124/1‑62 (MQ=255)                                                                                  gTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAc                                                                                                                <  2:999124/62‑1 (MQ=255)                                                                                  gTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGt                                   >  1:1074811/1‑139 (MQ=255)                                                                                     aTTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATATGGCGCGGtttt                                >  2:641807/1‑139 (MQ=255)                                                                                            aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAAcccc                                                                 >  2:1392200/1‑99 (MQ=255)                                                                                            aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAAcccc                                                                 <  1:1392200/99‑1 (MQ=255)                                                                                            aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTg                         >  1:634963/1‑139 (MQ=255)                                                                                                ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                 >  1:1123762/1‑47 (MQ=255)                                                                                                ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                 <  2:1123762/47‑1 (MQ=255)                                                                                                 gCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAg                    <  1:166542/139‑1 (MQ=255)                                                                                                 gCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAg                    <  1:1255032/139‑1 (MQ=255)                                                                                                   gACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAAcc                                                                   >  1:607079/1‑90 (MQ=255)                                                                                                   gACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAAcc                                                                   <  2:607079/90‑1 (MQ=255)                                                                                                        aTGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTaaaa                                                                     >  2:1430488/1‑83 (MQ=255)                                                                                                        aTGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTaaaa                                                                     <  1:1430488/83‑1 (MQ=255)                                                                                                           ccaccaTGTTTGACGTTAACAAAAAGTAATTCCAGACTCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGgcg                                       <  1:458045/110‑1 (MQ=255)                                                                                                           ccaccaTGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGCGATATAAGgcg                                       >  2:458045/1‑110 (MQ=255)                                                                                                                   tttGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCaa  <  2:288361/139‑1 (MQ=255)                                                                                                                      gACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCtgt         >  2:908638/1‑129 (MQ=255)                                                                                                                      gACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCtgt         <  1:908638/129‑1 (MQ=255)                                                                                                                        cGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCa   <  1:966076/133‑1 (MQ=255)                                                                                                                        cGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCa   >  2:966076/1‑133 (MQ=255)                                                                                                                             aCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCAc                            >  2:6600/1‑103 (MQ=255)                                                                                                                             aCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCAc                            <  1:6600/103‑1 (MQ=255)                                                                                                                                       |                                                                                                                       GTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAA  >  NC_000913/1502991‑1503243

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