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%	154.2 / NA	43	S97P (TCA→CCA)	ybaY	outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/24);  total (19/24)

CGTCTGGATCCGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGTCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGA  >  NC_000913/475537‑475776                                                                                                                                   |                                                                                                              cGTCTGGATCCGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTc                                                                                                       >  1:828085/1‑139 (MQ=255)   tCTGGATCCGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAg                                                                                                     >  2:479878/1‑139 (MQ=255)      ggATCCGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTGAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCtt                                                                                                  >  2:891851/1‑139 (MQ=255)        aTCCGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTg                                                                                                >  2:530133/1‑139 (MQ=255)          ccGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGtt                                                                                              <  2:1437817/139‑1 (MQ=255)          ccGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGtt                                                                                              <  1:479878/139‑1 (MQ=255)           cGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTc                                                                                             <  1:580448/139‑1 (MQ=255)           cGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTc                                                                                             <  1:530133/139‑1 (MQ=255)             tCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCAtt                                                                                                        >  1:1011291/1‑126 (MQ=255)             tCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCAtt                                                                                                        <  2:1011291/126‑1 (MQ=255)                  aaGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCAtt                                                                                      <  2:828085/139‑1 (MQ=255)                     tCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTaa                                                                                   <  2:106501/139‑1 (MQ=255)                                               aCCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAAcgcgc                                                         <  1:1327222/139‑1 (MQ=255)                                                 cGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGt                                                       <  1:891851/139‑1 (MQ=255)                                                    gACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTAtt                                                    >  2:1357620/1‑139 (MQ=255)                                                            cTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTa                                                      <  2:1325642/129‑1 (MQ=255)                                                            cTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTa                                                      >  1:1325642/1‑129 (MQ=255)                                                                 gcgTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTgcgg                                       >  2:748438/1‑139 (MQ=255)                                                                    tcgtTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAAcgc                                                           >  1:954588/1‑116 (MQ=255)                                                                    tcgtTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAAcgc                                                           <  2:954588/116‑1 (MQ=255)                                                                       ttAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTa                                 <  1:421638/139‑1 (MQ=255)                                                                         aGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTAcc                               <  1:748438/139‑1 (MQ=255)                                                                              aTGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGcct                                                                <  2:939238/101‑2 (MQ=255)                                                                              aTGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTAGTTCTGTCATTTAACCCGGCAGATGTTCAGCCg                                                                >  1:939238/1‑101 (MQ=255)                                                                               tGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGa                                    >  2:168918/1‑128 (MQ=255)                                                                               tGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGa                                    <  1:168918/128‑1 (MQ=255)                                                                                 cACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGa                       <  1:1096296/139‑1 (MQ=255)                                                                                   ccGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACa                     <  1:1357620/139‑1 (MQ=255)                                                                                         aaaGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAAc                                                                                  >  1:542996/1‑72 (MQ=255)                                                                                         aaaGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAAc                                                                                  <  2:542996/72‑1 (MQ=255)                                                                                          aaGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGt              >  2:901151/1‑139 (MQ=255)                                                                                                    gcAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGAt                                   <  1:480534/108‑1 (MQ=255)                                                                                                    gcAGAAAGCGGTGCGTACTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGAt                                   >  2:480534/1‑108 (MQ=255)                                                                                                                    aCTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGa  <  2:1394417/125‑1 (MQ=255)                                                                                                                    aCTGAAGGTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGa  >  1:1394417/1‑125 (MQ=255)                                                                                                                          ggTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGt                                                       >  2:180421/1‑66 (MQ=255)                                                                                                                          ggTAAACAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGt                                                       <  1:180421/66‑1 (MQ=255)                                                                                                                               aCAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAAcgc                                                           <  2:339858/57‑1 (MQ=255)                                                                                                                               aCAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAAcgc                                                           >  1:339858/1‑57 (MQ=255)                                                                                                                               aCAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTgcgg                                       <  2:345041/77‑1 (MQ=255)                                                                                                                               aCAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTgcgg                                       >  1:345041/1‑77 (MQ=255)                                                                                                                                cAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGaa                                                              <  2:1078760/53‑1 (MQ=255)                                                                                                                                cAGCCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGaa                                                              >  1:1078760/1‑53 (MQ=255)                                                                                                                                   |                                                                                                              CGTCTGGATCCGTCAGAAAGTCGCACTGCCGCCTGATGCTGTGCTGACCGTGACACTTTCTGACGCGTCGTTAGCCGATGCACCGTCAAAAGTGCTGGCGCAGAAAGCGGTGCGTACTGAAGGTAAACAGTCACCATTCAGCTTTGTTCTGTCATTTAACCCGGCAGATGTTCAGCCGAACGCGCGTATTCTGTTGAGTGCGGCGATTACCGTGAATGACAAACTGGTATTTATCACCGA  >  NC_000913/475537‑475776

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