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%	216.5 / NA	67	*223* (TAG→TAA)	ydcL	lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/32);  total (35/32)

TATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGG  >  NC_000913/1502993‑1503251                                                                                                                                     |                                                                                                                               tATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAg                                                                                                                          >  1:515907/1‑139 (MQ=255)  aTCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                         <  2:515907/139‑1 (MQ=255)    cAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcg                                                                                                                       <  1:888349/139‑1 (MQ=255)            ttcgtCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                                     >  2:214162/1‑133 (MQ=255)            ttcgtCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                                     <  1:214162/133‑1 (MQ=255)                    ggCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                                     <  2:1536011/125‑1 (MQ=255)                    ggCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                                     >  1:1536011/1‑125 (MQ=255)                              aaGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTAt                                                                                             >  1:1092533/1‑139 (MQ=255)                                   cTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGt                                                                                        >  2:809761/1‑139 (MQ=255)                                     gAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                    >  2:460134/1‑109 (MQ=255)                                     gAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                    <  1:460134/109‑1 (MQ=255)                                     gAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTgc                                                                                      >  2:17051/1‑139 (MQ=255)                                      aaTAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                    <  1:1646899/108‑1 (MQ=255)                                      aaTAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                    >  2:1646899/1‑108 (MQ=255)                                           cGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTa                                                                                >  1:865917/1‑139 (MQ=255)                                             aaaGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTaaa                                                                              >  2:191434/1‑139 (MQ=255)                                                gTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                    <  2:1487500/98‑1 (MQ=255)                                                gTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCATTAGATGGCGCTTTTTATCTGGTGCGCCGTAAAAcc                                                                           >  2:428166/1‑139 (MQ=255)                                                gTACGCCAATGGCTTTCGAAAATAATAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                    >  1:1487500/1‑98 (MQ=255)                                                  aCGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAAcccc                                                                         <  1:191434/139‑1 (MQ=255)                                                      cAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATcc                                                                     >  2:1644663/1‑139 (MQ=255)                                                         tGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTgc                                                                                      <  2:1384952/119‑1 (MQ=255)                                                         tGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTgc                                                                                      >  1:1384952/1‑119 (MQ=255)                                                             tttCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                              <  1:17051/139‑1 (MQ=255)                                                             tttCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                              <  1:809761/139‑1 (MQ=255)                                                             tttCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                              <  2:865917/139‑1 (MQ=255)                                                             tttCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                              <  1:1644663/139‑1 (MQ=255)                                                               tCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCg                                                            <  2:1092533/139‑1 (MQ=255)                                                                                    ttGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTc                                       >  1:473958/1‑139 (MQ=255)                                                                                          aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTaaaa                                                                             >  1:859622/1‑95 (MQ=255)                                                                                          aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTaaaa                                                                             <  2:859622/95‑1 (MQ=255)                                                                                           cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                    >  1:317340/1‑55 (MQ=255)                                                                                           cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                    <  2:317340/55‑1 (MQ=255)                                                                                            aTGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGAc                               >  1:275954/1‑139 (MQ=255)                                                                                              ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                    <  2:438058/52‑1 (MQ=255)                                                                                              ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                    <  2:1082070/52‑1 (MQ=255)                                                                                              ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                    >  1:1082070/1‑52 (MQ=255)                                                                                              ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                    >  1:438058/1‑52 (MQ=255)                                                                                              ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTa                                                                                >  2:933281/1‑88 (MQ=255)                                                                                              ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTa                                                                                <  1:933281/88‑1 (MQ=255)                                                                                              ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCa                             >  2:373956/1‑139 (MQ=255)                                                                                               tcaaCCGATGCCACCTTGTTTGACGTTAACAAAAACTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTTTCCGATGCGCCTTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAg                            <  1:428166/136‑1 (MQ=255)                                                                                                 gACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                    >  2:1142388/1‑49 (MQ=255)                                                                                                 gACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                    <  1:1142388/49‑1 (MQ=255)                                                                                                  aCCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                              >  1:284603/1‑102 (MQ=255)                                                                                                  aCCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                              <  2:284603/102‑1 (MQ=255)                                                                                                   ccGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGAc                               >  2:120928/1‑132 (MQ=255)                                                                                                   ccGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGAc                               <  1:120928/132‑1 (MQ=255)                                                                                                   ccGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGAc                               <  1:210222/132‑1 (MQ=255)                                                                                                   ccGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGAc                               >  2:210222/1‑132 (MQ=255)                                                                                                      aTGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGc                                                             >  2:364727/1‑99 (MQ=255)                                                                                                      aTGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGc                                                             <  1:364727/99‑1 (MQ=255)                                                                                                           accaTGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCtgtt                >  1:1662707/1‑139 (MQ=255)                                                                                                                gTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCa           <  2:784127/139‑1 (MQ=255)                                                                                                                 tttGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTg                                                                                       >  1:1255953/1‑62 (MQ=255)                                                                                                                 tttGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTg                                                                                       <  2:1255953/62‑1 (MQ=255)                                                                                                                 tttGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGgcgc                                              <  1:160277/103‑1 (MQ=255)                                                                                                                 tttGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGgcgc                                              >  2:160277/1‑103 (MQ=255)                                                                                                                 tttGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCAc                                    >  1:1278744/1‑113 (MQ=255)                                                                                                                 tttGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCAc                                    <  2:1278744/113‑1 (MQ=255)                                                                                                                   tGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTa                                                                                <  1:1282258/67‑1 (MQ=255)                                                                                                                   tGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTa                                                                                >  2:1282258/1‑67 (MQ=255)                                                                                                                   tGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCg                                                            <  2:757372/87‑1 (MQ=255)                                                                                                                   tGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCg                                                            >  1:757372/1‑87 (MQ=255)                                                                                                                         tAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTgg  >  1:1647575/1‑139 (MQ=255)                                                                                                                          aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCg                                                            <  2:247471/80‑1 (MQ=255)                                                                                                                          aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCg                                                            >  1:247471/1‑80 (MQ=255)                                                                                                                                     |                                                                                                                               TATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGG  >  NC_000913/1502993‑1503251

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