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%	184.6 / NA	57	M143I (ATG→ATT)	ybbO	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (32/25);  total (32/25)

GTTCGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATA  >  NC_000913/518585‑518852                                                                                                                                         |                                                                                                                                    gTTCGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGtaat                                                                                                                                   <  1:1282818/139‑1 (MQ=255)   tCGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGtaatca                                                                                                                                 <  1:589270/139‑1 (MQ=255)       tCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAAt                                                                                                                             <  2:190361/139‑1 (MQ=255)           gCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGc                                                                                                                         >  2:66140/1‑139 (MQ=255)           gCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGc                                                                                                                         >  1:248224/1‑139 (MQ=255)               aCTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCtt                                                                                                                     <  1:807195/139‑1 (MQ=255)                 tttAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCa                                                                                                                   <  1:1281638/139‑1 (MQ=255)                                gcgcAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATc                                                                                                    >  2:721268/1‑139 (MQ=255)                                 cgcAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATcg                                                                                                   <  1:721268/139‑1 (MQ=255)                                           tGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATa                                                                                                                            >  2:657587/1‑104 (MQ=255)                                           tGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATa                                                                                                                            <  1:657587/104‑1 (MQ=255)                                                aGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCa                                                                                    <  1:83905/139‑1 (MQ=255)                                                  tGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGgcg                                                                                      >  2:700308/1‑135 (MQ=255)                                                  tGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGgcg                                                                                      <  1:700308/135‑1 (MQ=255)                                                            cACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCg                                                                                                                >  2:1503282/1‑99 (MQ=255)                                                            cACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCg                                                                                                                <  1:1503282/99‑1 (MQ=255)                                                            cACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGgt                                                                        >  1:1782905/1‑139 (MQ=255)                                                            cACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGgt                                                                        >  1:776391/1‑139 (MQ=255)                                                               gCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATAc                                                                                                                           <  1:185275/85‑1 (MQ=255)                                                               gCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATAc                                                                                                                           >  2:185275/1‑85 (MQ=255)                                                                    cAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCa                                                                                                                   <  2:529014/88‑1 (MQ=255)                                                                    cAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCa                                                                                                                   >  1:529014/1‑88 (MQ=255)                                                                    cAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCg                                                                >  1:316254/1‑139 (MQ=255)                                                                         cATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACa                                                                                                      <  1:1394026/96‑1 (MQ=255)                                                                         cATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACa                                                                                                      >  2:1394026/1‑96 (MQ=255)                                                                         cATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATg                                                                                  >  2:483396/1‑116 (MQ=255)                                                                         cATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATg                                                                                  <  1:483396/116‑1 (MQ=255)                                                                            atTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGtt                                                        >  2:744071/1‑139 (MQ=255)                                                                            atTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGtt                                                        >  2:1378097/1‑139 (MQ=255)                                                                                    gCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATcaca                                                                                                                               <  2:1093852/60‑1 (MQ=255)                                                                                    gCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGaaa                                                >  2:706257/1‑139 (MQ=255)                                                                                    gCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGaa                                                 <  1:605646/138‑1 (MQ=255)                                                                                    gCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGaa                                                 >  2:605646/1‑138 (MQ=255)                                                                                    gCCGCGTAAGCGCCACGACCCGGCGTGCAGATTAATCCCATCACCGATGATGTAATcaca                                                                                                                               >  1:1093852/1‑60 (MQ=255)                                                                                      cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAg                                                                             >  1:206010/1‑108 (MQ=255)                                                                                      cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAg                                                                             <  2:206010/108‑1 (MQ=255)                                                                                      cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGt                                                         >  1:1800263/1‑128 (MQ=255)                                                                                      cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGt                                                         <  2:1800263/128‑1 (MQ=255)                                                                                      cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAc                                              >  2:252985/1‑139 (MQ=255)                                                                                      cgcgAAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTACTCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAc                                              >  1:1064390/1‑139 (MQ=255)                                                                                        cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCg                                                                                                                <  2:1415171/71‑1 (MQ=255)                                                                                        cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCg                                                                                                                >  1:1415171/1‑71 (MQ=255)                                                                                        cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAActg                                            <  1:1378097/139‑1 (MQ=255)                                                                                        cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAActg                                            >  2:920111/1‑139 (MQ=255)                                                                                        cgTAAGCGCCACGACCCGGCGTGCAGATTAATCCCATCACCGATGTTGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACACGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAATctg                                            >  2:1056674/1‑139 (MQ=255)                                                                                             gcgcCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGtt                                       <  2:776391/139‑1 (MQ=255)                                                                                                   cGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCAt                                   <  2:975158/137‑1 (MQ=255)                                                                                                   cGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCAt                                   >  1:975158/1‑137 (MQ=255)                                                                                                           cGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCa                            >  2:1610480/1‑136 (MQ=255)                                                                                                           cGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCa                            <  1:1610480/136‑1 (MQ=255)                                                                                                              ggAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTg                      <  2:316254/139‑1 (MQ=255)                                                                                                                     aaTCCCATCACCGATGATGTATCCCCAAAACGCCCTTCCCCGTGCGGTAACATCGCGGGAAACAGCCGCATGGTGAGCTGTTGTGCGCCGAAAAAGTTGGCGGAAAACTGCGGTTCCACCTCCGAAGGGCGGAtggggg               >  1:684016/1‑139 (MQ=255)                                                                                                                         ccATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAg           <  1:920111/139‑1 (MQ=255)                                                                                                                           aTCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAct                                             <  2:510905/103‑1 (MQ=255)                                                                                                                           aTCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAct                                             >  1:510905/1‑103 (MQ=255)                                                                                                                            tCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCAc                           >  1:970081/1‑120 (MQ=255)                                                                                                                            tCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCAc                           <  2:970081/120‑1 (MQ=255)                                                                                                                                cGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCa    >  2:1252924/1‑139 (MQ=255)                                                                                                                                  atgatgTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATa  >  2:4436/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                    GTTCGATCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATA  >  NC_000913/518585‑518852

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