Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	803,662	C→A	L54I (CTT→ATT)	ybhJ →	putative hydratase YbhJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	803,662	0	C	A	100.0%	288.1 / NA	75	L54I (CTT→ATT)	ybhJ	putative hydratase YbhJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (30/45);  total (30/45)

CCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGC  >  NC_000913/803597‑803726                                                                   |                                                                 ccAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAAtt                                                                >  1:536597/1‑68 (MQ=255) ccAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAAtt                                                                >  1:1466296/1‑68 (MQ=255)     aaaaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTaaaa                                                            <  1:276916/68‑1 (MQ=255)       aaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑Taaaattaaaatt                                                          >  1:2405829/1‑68 (MQ=255)           cACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGt                                                      >  1:2647085/1‑68 (MQ=255)           cACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGt                                                      <  1:1393658/68‑1 (MQ=255)                 tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGAtt                                                >  1:1044354/1‑68 (MQ=255)                  gCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTc                                               <  1:1014343/68‑1 (MQ=255)                  gCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTc                                               <  1:2248442/68‑1 (MQ=255)                   ccTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCa                                              <  1:1890268/68‑1 (MQ=255)                   ccTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCa                                              <  1:3333318/68‑1 (MQ=255)                   ccTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCa                                              <  1:3224038/68‑1 (MQ=255)                    cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCAt                                             <  1:2256065/68‑1 (MQ=255)                    cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCAt                                             >  1:2172792/1‑68 (MQ=255)                    cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCAt                                             >  1:2708582/1‑68 (MQ=255)                    cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCAt                                             >  1:1524477/1‑68 (MQ=255)                      ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTa                                           >  1:1098089/1‑68 (MQ=255)                        tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGc                                         >  1:3112738/1‑68 (MQ=255)                        tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGc                                         >  1:2065042/1‑68 (MQ=255)                        tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGc                                         >  1:467451/1‑68 (MQ=255)                        tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGc                                         >  1:1786276/1‑68 (MQ=255)                        tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGc                                         <  1:1175575/68‑1 (MQ=255)                         cTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGcc                                        <  1:223507/68‑1 (MQ=255)                           aTTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCctc                                      <  1:2831140/68‑1 (MQ=255)                            ttCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCctct                                     <  1:393600/68‑1 (MQ=255)                            ttCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCctct                                     <  1:2493097/68‑1 (MQ=255)                             tctctcTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:3230040/1‑68 (MQ=255)                             tctctcTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:1141118/1‑68 (MQ=255)                             tctctcTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:849806/1‑68 (MQ=255)                             tctctcTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:1238363/1‑68 (MQ=255)                             tctctcTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:1868927/1‑68 (MQ=255)                             tctctcTTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                    >  1:1576953/1‑68 (MQ=255)                                ctctTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACg                                 <  1:1933981/68‑1 (MQ=255)                                ctctTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACg                                 <  1:1848441/68‑1 (MQ=255)                                ctctTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACg                                 <  1:3134453/68‑1 (MQ=255)                                ctctTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACg                                 >  1:3355252/1‑68 (MQ=255)                                 tctTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                                <  1:2894596/68‑1 (MQ=255)                                  ctTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAGAATTAAGTTTGATTCATTCGCCTCTCACGat                               >  1:203965/1‑68 (MQ=255)                                  ctTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                               <  1:982194/68‑1 (MQ=255)                                  ctTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                               <  1:2416866/68‑1 (MQ=255)                                  ctTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                               >  1:825856/1‑68 (MQ=255)                                  ctTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                               <  1:3215011/68‑1 (MQ=255)                                  ctTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                               <  1:136732/68‑1 (MQ=255)                                  ctTCGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                               <  1:141834/68‑1 (MQ=255)                                  ctTCGCATAATACGTCCGGAAATATGGA‑TAAAATCAAAATTAAGTTTGATTCATTAGCCTCTCACGat                               <  1:959135/68‑1 (MQ=255)                                  cagcGCATAATACGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                               <  1:1190867/65‑1 (MQ=255)                                        ataataCGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcc                         <  1:2153692/68‑1 (MQ=255)                                        ataataCGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcc                         <  1:2534867/68‑1 (MQ=255)                                           ataCGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCttt                      >  1:2876569/1‑68 (MQ=255)                                            taCGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttg                     >  1:2209381/1‑68 (MQ=255)                                            taCGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttg                     <  1:2550204/68‑1 (MQ=255)                                             aCGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgt                    >  1:1018128/1‑68 (MQ=255)                                             aCGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgt                    >  1:2283025/1‑68 (MQ=255)                                             aCGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgt                    >  1:2009988/1‑68 (MQ=255)                                              cGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgtt                   >  1:1665307/1‑68 (MQ=255)                                              cGTCCGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgtt                   >  1:1281424/1‑68 (MQ=255)                                                 ccGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTGCCTTTGTTGGt                <  1:1289746/68‑1 (MQ=255)                                                 ccGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:2931314/68‑1 (MQ=255)                                                 ccGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:3016737/68‑1 (MQ=255)                                                 ccGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:1590815/68‑1 (MQ=255)                                                 ccGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:3175378/68‑1 (MQ=255)                                                 ccGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:609918/68‑1 (MQ=255)                                                 ccGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:1735040/68‑1 (MQ=255)                                                 ccGGAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt                <  1:2047617/68‑1 (MQ=255)                                                   ggAAATATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt              <  1:857191/68‑1 (MQ=255)                                                     aaaTATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTg            >  1:1422984/1‑68 (MQ=255)                                                      aaaatGGATTAAAATTAAATTTAA‑TTTGAATCAATAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:447644/65‑1 (MQ=255)                                                      aaTATGGA‑TAAAATTAAAATTAATTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:2167495/68‑1 (MQ=255)                                                      aaTATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:1277347/68‑1 (MQ=255)                                                      aaTATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:2639669/68‑1 (MQ=255)                                                      aaTATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:3179415/68‑1 (MQ=255)                                                      aaTATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:1155063/68‑1 (MQ=255)                                                      aaTATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:838985/68‑1 (MQ=255)                                                      aaTATGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt           <  1:1148295/68‑1 (MQ=255)                                                      aaTATGGA‑TAAAATCAAAATTAAGTTTGATTCATTAGCCTCTCACGATGTTACCTTTGTTGGTATTGt           <  1:479872/68‑1 (MQ=255)                                                       atatGGA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTa          >  1:1624800/1‑68 (MQ=255)                                                           ggA‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGa      <  1:74858/68‑1 (MQ=255)                                                             a‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTGCCTTTGTTGGTATTGTACAGAcc    <  1:241580/68‑1 (MQ=255)                                                             a‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGAcc    <  1:2692614/68‑1 (MQ=255)                                                             a‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGAcc    <  1:1281419/68‑1 (MQ=255)                                                             a‑TAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGAcc    <  1:1630577/68‑1 (MQ=255)                                                                aaaattaaaattAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGc  <  1:89213/68‑1 (MQ=255)                                                                   |                                                                 CCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGA‑TAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGC  >  NC_000913/803597‑803726

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑