Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	475,667	T→C	S97P (TCA→CCA)	ybaY →	outer membrane lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	475,667	0	T	C	100.0%	144.0 / NA	41	S97P (TCA→CCA)	ybaY	outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/21);  total (20/21)

GTACCGTCTGGATCCGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGTCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATC  >  NC_000913/475533‑475795                                                                                                                                       |                                                                                                                                 gTACCGTCTGGATCCGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCAcc                                                                                                                              >  1:74518/1‑139 (MQ=255)    ccGTCTGGATCCGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCAtt                                                                                                                           <  2:88341/139‑1 (MQ=255)     cGTCTGGATCCGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTCCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTc                                                                                                                          >  1:782011/1‑139 (MQ=255)       tCTGGATCCGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCaa                                                                                                                        >  1:319281/1‑138 (MQ=255)       tCTGGATCCGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAg                                                                                                                        <  1:547639/139‑1 (MQ=255)       tCTGGATCCGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAg                                                                                                                        >  2:629842/1‑139 (MQ=255)       tCTGGATCCGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAg                                                                                                                        >  1:756921/1‑139 (MQ=255)               cGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACAGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTc                                                                                                                <  1:628292/139‑1 (MQ=255)                 tCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGTCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTg                                                                                                              >  2:957227/1‑139 (MQ=255)                                              tgctgACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGaa                                                                                 <  2:74518/139‑1 (MQ=255)                                                   aCCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAAcgcgc                                                                            >  1:360453/1‑139 (MQ=255)                                                   aCCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAAcgcgc                                                                            <  2:557964/139‑1 (MQ=255)                                                             tttCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGtt                                                                  <  2:771167/139‑1 (MQ=255)                                                               tCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGa                                                                >  1:96442/1‑139 (MQ=255)                                                                       gtcgtTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAAcgcgc                                                                            >  1:112063/1‑119 (MQ=255)                                                                       gtcgtTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAAcgcgc                                                                            <  2:112063/119‑1 (MQ=255)                                                                               ccGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAACTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGt                                                <  2:319281/139‑1 (MQ=255)                                                                                   tGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAg                                                                                      <  2:58787/97‑1 (MQ=255)                                                                                   tGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAg                                                                                      >  1:58787/1‑97 (MQ=255)                                                                                   tGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAAt                                            >  1:402916/1‑139 (MQ=255)                                                                                   tGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAAt                                            >  2:155188/1‑139 (MQ=255)                                                                                     cACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGa                                          <  2:324004/139‑1 (MQ=255)                                                                                       ccGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACa                                        <  2:173879/139‑1 (MQ=255)                                                                                          tCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAAcgc                                                                              >  1:302889/1‑98 (MQ=255)                                                                                          tCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAAcgc                                                                              <  2:302889/98‑1 (MQ=255)                                                                                                    tGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTAt                           <  1:957227/139‑1 (MQ=255)                                                                                                          aGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGa                     <  2:360453/139‑1 (MQ=255)                                                                                                                cGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAAc                                                                                                     >  1:776602/1‑53 (MQ=255)                                                                                                                cGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAAc                                                                                                     <  2:776602/53‑1 (MQ=255)                                                                                                                          tGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTg     <  2:402916/139‑1 (MQ=255)                                                                                                                           gAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGa    <  2:96442/139‑1 (MQ=255)                                                                                                                            aaGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGtt                                                                  <  2:641034/76‑1 (MQ=255)                                                                                                                            aaGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGtt                                                                  <  2:102060/76‑1 (MQ=255)                                                                                                                            aaGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGtt                                                                  >  1:641034/1‑76 (MQ=255)                                                                                                                            aaGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGtt                                                                  >  1:102060/1‑76 (MQ=255)                                                                                                                            aaGGTAAACAGCCACCATTCAGCTTTGTTATGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATGTATCACCGATACCGTTCAGCCGGTGAt   >  1:582311/1‑139 (MQ=255)                                                                                                                             aGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATc  >  2:224862/1‑139 (MQ=255)                                                                                                                                 aaaCAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGtt              <  1:596459/123‑1 (MQ=255)                                                                                                                                 aaaCAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGtt              >  2:596459/1‑123 (MQ=255)                                                                                                                                 aaaCAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGc          >  2:133070/1‑127 (MQ=255)                                                                                                                                 aaaCAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGc          <  1:133070/127‑1 (MQ=255)                                                                                                                                       |                                                                                                                                 GTACCGTCTGGATCCGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGTCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGATACCGTTCAGCCGGTGATC  >  NC_000913/475533‑475795

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