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

TCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATG  >  NC_000913/518591‑518856                                                                                                                                   |                                                                                                                                        tCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAAt                                                                                                                                 <  1:293354/139‑1 (MQ=255)   aGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATAc                                                                                                                               <  1:101802/139‑1 (MQ=255)           tttAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCa                                                                                                                       <  1:192421/139‑1 (MQ=255)            ttAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCAc                                                                                                                      <  1:706854/139‑1 (MQ=255)                         ggCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGGTGTAATCACAATACGCCCTTCACCGTGCGGTAACAt                                                                                                         <  1:1366009/139‑1 (MQ=255)                           cgcAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATcg                                                                                                       <  1:1031830/139‑1 (MQ=255)                                   cATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAAc                                                                                               <  1:541222/139‑1 (MQ=255)                                           gTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCAt                                                                                       >  2:1007859/1‑139 (MQ=255)                                            tGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATg                                                                                      <  2:1267111/139‑1 (MQ=255)                                                                   cATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATc                                                                                                        <  1:1333544/98‑1 (MQ=255)                                                                   cATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATc                                                                                                        >  2:1333544/1‑98 (MQ=255)                                                                          aCTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCg                                                        >  2:1432321/1‑139 (MQ=255)                                                                              gCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCg                                                                                                                    >  1:958348/1‑75 (MQ=255)                                                                              gCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCg                                                                                                                    <  2:958348/75‑1 (MQ=255)                                                                              gCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCg                                                                                                                    >  2:853994/1‑75 (MQ=255)                                                                              gCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCg                                                                                                                    <  1:853994/75‑1 (MQ=255)                                                                              gCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGaaa                                                    >  2:1397024/1‑139 (MQ=255)                                                                              gCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGaaa                                                    >  2:907956/1‑139 (MQ=255)                                                                               ccGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGACGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGtgtg                                                                         <  2:1005154/117‑1 (MQ=255)                                                                               acgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGtgtg                                                                         >  1:1005154/2‑117 (MQ=255)                                                                                cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTctc                                                                                                                      <  1:967408/71‑3 (MQ=255)                                                                                cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCAc                                                                                                                      >  2:967408/1‑71 (MQ=255)                                                                                cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCTTGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAc                                                  >  1:793252/1‑139 (MQ=255)                                                                                cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAc                                                  <  1:1397024/139‑1 (MQ=255)                                                                                cgcgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAc                                                  >  1:808194/1‑139 (MQ=255)                                                                                  cgTAAGCTCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGaa                                                     >  2:1419476/1‑134 (MQ=255)                                                                                  cgTAAGCTCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGaa                                                     <  1:1419476/134‑1 (MQ=255)                                                                                  cgTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAActg                                                >  2:1503211/1‑139 (MQ=255)                                                                                   gTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAActgc                                               >  2:650101/1‑139 (MQ=255)                                                                                     aaGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTATCAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAACActgctg                                             <  2:985513/139‑1 (MQ=255)                                                                                             cGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGt                                                             >  1:1450272/1‑115 (MQ=255)                                                                                             cGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGt                                                             <  2:1450272/115‑1 (MQ=255)                                                                                                   ggCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCAc                               >  1:539239/1‑139 (MQ=255)                                                                                                   ggCGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCAc                               >  1:492188/1‑139 (MQ=255)                                                                                                     cGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTgc                                   >  1:487144/1‑133 (MQ=255)                                                                                                     cGTGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTgc                                   <  2:487144/133‑1 (MQ=255)                                                                                                       tGGAGATTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCt                           <  2:580094/139‑1 (MQ=255)                                                                                                          agaTTAATCCCATCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAt                        >  1:832761/1‑139 (MQ=255)                                                                                                                      tCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAActgc                                               <  2:64241/104‑1 (MQ=255)                                                                                                                      tCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAActgc                                               >  1:64241/1‑104 (MQ=255)                                                                                                                      tCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAct                                                 >  1:1437861/1‑102 (MQ=255)                                                                                                                      tCACCGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAAct                                                 <  2:1437861/102‑1 (MQ=255)                                                                                                                          cGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAg                                                                                 <  2:1138795/66‑1 (MQ=255)                                                                                                                          cGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAg                                                                                 >  1:1138795/1‑66 (MQ=255)                                                                                                                          cGATGATGTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCa        <  1:745678/139‑1 (MQ=255)                                                                                                                            atgatgTAATCACAATACGCCCTTCACCGTGCGGTAACa                                                                                                          <  1:130537/39‑1 (MQ=37)                                                                                                                            atgatgTAATCACAATACGCCCTTCACCGTGCGGTAACa                                                                                                          >  2:130537/1‑39 (MQ=37)                                                                                                                            atgatgTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAg                                                                                 <  2:526843/64‑1 (MQ=255)                                                                                                                            atgatgTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAg                                                                                 >  1:526843/1‑64 (MQ=255)                                                                                                                            atgatgTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATa      >  2:114447/1‑139 (MQ=255)                                                                                                                               atgTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGgc                                                                                           <  2:870815/51‑1 (MQ=255)                                                                                                                               atgTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGgc                                                                                           >  1:870815/1‑51 (MQ=255)                                                                                                                                tgTAATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATg  <  2:1087225/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                        TCAGGCTGACTTTAATTCCGCTGTGGCGCAGCTCCATGCGCAGTGCATCTGACCACGCCTCCAGCGCATATTTACTGGCCGCGTAAGCGCCACGACCCGGCGTGGAGATTAATCCCATCACCGATGATGTCATCACAATACGCCCTTCACCGTGCGGTAACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATG  >  NC_000913/518591‑518856

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