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%	316.6 / NA	83	L54I (CTT→ATT)	ybhJ	putative hydratase YbhJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (30/53);  total (30/53)

AAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGA  >  NC_000913/803602‑803722                                                             |                                                             aaaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTaa                                                          >  1:2752146/1‑65 (MQ=255) aaaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTaa                                                          >  1:141575/1‑65 (MQ=255)  aaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTaaa                                                         >  1:1116325/1‑65 (MQ=255)  aaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTaaa                                                         >  1:757423/1‑65 (MQ=255)   aGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTaaaa                                                        <  1:1553679/65‑1 (MQ=255)   aGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTaaaa                                                        <  1:628319/65‑1 (MQ=255)      cACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAaaattaaaatta                                                     >  1:469459/1‑65 (MQ=255)         tATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGt                                                  <  1:2829708/65‑1 (MQ=255)         tATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGt                                                  <  1:888479/65‑1 (MQ=255)         tATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATAGATAAAATTAAAATTAAGt                                                  <  1:1364106/65‑1 (MQ=255)          aTTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGtt                                                 >  1:2335112/1‑65 (MQ=255)            tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTg                                               >  1:1715850/1‑65 (MQ=255)            tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTg                                               <  1:237055/65‑1 (MQ=255)            tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTg                                               >  1:391204/1‑65 (MQ=255)            tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTg                                               >  1:1283863/1‑65 (MQ=255)               cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                            >  1:2871884/1‑65 (MQ=255)               cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                            >  1:2258952/1‑65 (MQ=255)                tGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTc                                           <  1:2336601/65‑1 (MQ=255)                  gTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAt                                         <  1:1965801/65‑1 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAGATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                        >  1:2502495/1‑65 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                        >  1:62613/1‑65 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                        >  1:1388692/1‑65 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                        >  1:1257453/1‑65 (MQ=255)                      aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                     <  1:592840/65‑1 (MQ=255)                      aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                     <  1:1313623/65‑1 (MQ=255)                      aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                     <  1:2897687/65‑1 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:2593943/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:2671173/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:2482347/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:2481788/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:174854/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:467588/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:616866/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:1232017/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:1210088/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:782294/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:810109/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:1104313/1‑65 (MQ=255)                         ctctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctc                                  >  1:3063528/1‑65 (MQ=255)                         ctctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctc                                  <  1:2718802/65‑1 (MQ=255)                         ctctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctc                                  <  1:589271/65‑1 (MQ=255)                         ctctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctc                                  >  1:893849/1‑65 (MQ=255)                           ccctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                >  1:2198579/3‑65 (MQ=255)                              ttCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACg                             <  1:961657/65‑1 (MQ=255)                               tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                            <  1:270553/65‑1 (MQ=255)                                cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                           <  1:2478316/65‑1 (MQ=255)                                cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                           <  1:2709108/65‑1 (MQ=255)                                cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                           <  1:2194426/65‑1 (MQ=255)                                cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                           <  1:1399717/65‑1 (MQ=255)                                cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                           <  1:2512553/65‑1 (MQ=255)                                      ataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcc                     >  1:3144951/1‑65 (MQ=255)                                        aCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCtt                   >  1:1442591/1‑65 (MQ=255)                                          gTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttg                 <  1:1036289/65‑1 (MQ=255)                                            ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgtt               <  1:1606741/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTCCCTTTGTTGGt            <  1:3079171/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:140026/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:140697/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:3149441/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:3078032/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:1388926/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:1191222/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:2988727/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:869438/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:1995728/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:1116481/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:2078453/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:95126/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:2198073/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:2562024/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:2509732/65‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCCCGATATTGCCTTTGTTGGTAt          <  1:445853/65‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:793675/65‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:1033730/65‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:3234436/65‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:2598970/65‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:2374742/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:1442361/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:594469/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:1576983/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:1278009/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:3013155/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:1822397/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:1124287/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:1961314/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:2047244/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATCCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:88336/65‑1 (MQ=255)                                                        aTAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAg   >  1:1800887/1‑65 (MQ=255)                                                         taaaattaaaatTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGa  <  1:2790254/65‑1 (MQ=255)                                                             |                                                             AAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGA  >  NC_000913/803602‑803722

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