Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	518,721	C→A	M143I (ATG→ATT)	ybbO ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	518,721	0	C	A	100.0%	127.4 / NA	41	M143I (ATG→ATT)	ybbO	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (23/18);  total (23/18)

TCGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACAT  >  NC_000913/518587‑518855                                                                                                                                       |                                                                                                                                       tCGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGtaatca                                                                                                                                    >  2:1154520/1‑139 (MQ=255)   gATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATcaca                                                                                                                                  >  1:1126528/1‑139 (MQ=255)               tttAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCa                                                                                                                      <  2:1126528/139‑1 (MQ=255)                ttAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCAc                                                                                                                     <  1:247163/139‑1 (MQ=255)                  aaTTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCg                                                                                                                   <  1:13666/139‑1 (MQ=255)                         cTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTa                                                                                                            >  2:1258846/1‑139 (MQ=255)                            tgGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACa                                                                                                         <  2:316561/139‑1 (MQ=255)                            tgGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACa                                                                                                         <  2:1415315/139‑1 (MQ=255)                               cgcAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATcg                                                                                                      <  2:816654/139‑1 (MQ=255)                                    cTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGt                                                                                                 >  1:339536/1‑139 (MQ=255)                                      ccATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTaa                                                                                               >  1:342931/1‑139 (MQ=255)                                              aGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGgcgc                                                                                        >  2:259918/1‑138 (MQ=255)                                              aGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGgcgc                                                                                        <  1:259918/138‑1 (MQ=255)                                              aGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCa                                                                                       <  1:566178/139‑1 (MQ=255)                                                   aTCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACCATACGCCCTTCACCg                                                                                                                   <  1:981326/106‑1 (MQ=255)                                                   aTCTGACCACGCCTCCAGCCCATATTTACGGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCg                                                                                                                   >  2:981326/1‑106 (MQ=255)                                                                  cAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCg                                                                   >  1:591140/1‑139 (MQ=255)                                                                       cATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGt                                                                                   >  2:1409566/1‑118 (MQ=255)                                                                       cATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGt                                                                                   <  1:1409566/118‑1 (MQ=255)                                                                               cTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCgg                                                      >  2:109147/1‑139 (MQ=255)                                                                                 ggCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGaa                                                    >  1:871283/1‑139 (MQ=255)                                                                                      cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGtaat                                                                                                                                      <  1:1171845/52‑1 (MQ=255)                                                                                      cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGtaat                                                                                                                                      >  2:1171845/1‑52 (MQ=255)                                                                                      cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGt                                                                                                             <  2:246376/77‑1 (MQ=255)                                                                                      cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGt                                                                                                             >  1:246376/1‑77 (MQ=255)                                                                                      cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAActg                                               >  1:1145842/1‑139 (MQ=255)                                                                                      cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAActg                                               >  2:888962/1‑139 (MQ=255)                                                                                                 cGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCt                                    >  1:918022/1‑139 (MQ=255)                                                                                                       ggCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTgcg                                 <  1:904055/136‑1 (MQ=255)                                                                                                       agCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTACTGTTCCATCTgcg                                 >  2:904055/2‑136 (MQ=255)                                                                                                        gCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACg                             >  2:620511/1‑139 (MQ=255)                                                                                                          gTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGc                           <  2:497307/139‑1 (MQ=255)                                                                                                            ggAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGgcg                                                                                         <  2:849940/75‑1 (MQ=255)                                                                                                            ggAGATTAATCCCATCACCGATGATGTAAGCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGgcg                                                                                         >  1:849940/1‑75 (MQ=255)                                                                                                              agaTTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAt                       <  2:339536/139‑1 (MQ=255)                                                                                                              agaTTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAt                       >  2:1361690/1‑139 (MQ=255)                                                                                                                  tAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTgcg                                                                      <  2:847354/88‑1 (MQ=255)                                                                                                                  tAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAtggt                    <  2:156088/138‑1 (MQ=255)                                                                                                                  tAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAtggt                    >  1:156088/1‑138 (MQ=255)                                                                                                                  tAATCCCATCACCGATGATGTAATCACAAGACGCCCTTCACCGTGAGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTgcg                                                                      >  1:847354/1‑88 (MQ=255)                                                                                                                   aaTCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCAGAAAAAGTTGGCGGAAAACTGCTGTTACATCTGCGCACGGCTGAtggtgg                  >  2:1088526/1‑139 (MQ=255)                                                                                                                                atgatgTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAct                                                <  1:1383671/96‑1 (MQ=255)                                                                                                                                atgatgTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAct                                                >  2:1383671/1‑96 (MQ=255)                                                                                                                                   atgTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACAt  >  1:1421118/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                       TCGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACAT  >  NC_000913/518587‑518855

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