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

AAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGA  >  NC_000913/803602‑803722                                                             |                                                             aaaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTaa                                                          >  1:1751017/1‑65 (MQ=255)    ggCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGAtaaaattaaaat                                                       >  1:2623144/1‑65 (MQ=255)        cTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAg                                                   <  1:780888/65‑1 (MQ=255)          aTTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGtt                                                 >  1:162539/1‑65 (MQ=255)          aTTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGtt                                                 >  1:3078581/1‑65 (MQ=255)          aTTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGtt                                                 >  1:2705882/1‑65 (MQ=255)          aTTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGtt                                                 >  1:2113601/1‑65 (MQ=255)            tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTg                                               >  1:1692185/1‑65 (MQ=255)            tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTg                                               >  1:1744914/1‑65 (MQ=255)            tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTg                                               >  1:90130/1‑65 (MQ=255)            tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTg                                               >  1:2909591/1‑65 (MQ=255)               cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAATATTAAGTTTGAtt                                            >  1:3076118/1‑65 (MQ=255)               cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                            >  1:903012/1‑65 (MQ=255)               cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                            >  1:2167640/1‑65 (MQ=255)                tGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTc                                           <  1:3108988/65‑1 (MQ=255)                 ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                          <  1:1705628/65‑1 (MQ=255)                 ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                          <  1:2021478/65‑1 (MQ=255)                 ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                          <  1:2320186/65‑1 (MQ=255)                 ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                          <  1:236901/65‑1 (MQ=255)                 ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                          <  1:1317109/65‑1 (MQ=255)                 ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                          >  1:86653/1‑65 (MQ=255)                 ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                          <  1:2536334/65‑1 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                        >  1:2182208/1‑65 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                        >  1:2148449/1‑65 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                        >  1:805393/1‑65 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                        >  1:1017710/1‑65 (MQ=255)                       ttCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGcc                                    <  1:2125567/65‑1 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:2267911/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:2589969/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:1075863/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:230383/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:1160581/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:2234902/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:2192623/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:1731206/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   <  1:2170647/65‑1 (MQ=255)                         ctctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctc                                  <  1:2564487/65‑1 (MQ=255)                         ctctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctc                                  <  1:580739/65‑1 (MQ=255)                         ctctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctc                                  <  1:2663968/65‑1 (MQ=255)                          tctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctct                                 <  1:1520349/65‑1 (MQ=255)                          tctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctct                                 <  1:607181/65‑1 (MQ=255)                           ctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                >  1:303872/1‑65 (MQ=255)                             ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCAc                              >  1:1734805/1‑65 (MQ=255)                             ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCAc                              >  1:3223855/1‑65 (MQ=255)                             ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCAc                              >  1:184268/1‑65 (MQ=255)                               tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                            <  1:1733603/65‑1 (MQ=255)                               tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                            <  1:2706309/65‑1 (MQ=255)                               tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                            <  1:2357910/65‑1 (MQ=255)                               tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                            <  1:840395/65‑1 (MQ=255)                               tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                            <  1:2202722/65‑1 (MQ=255)                               tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                            <  1:2188899/65‑1 (MQ=255)                                cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                           <  1:3149115/65‑1 (MQ=255)                                cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                           <  1:705960/65‑1 (MQ=255)                                cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                           <  1:2447882/65‑1 (MQ=255)                                cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                           <  1:2631043/65‑1 (MQ=255)                                   ataataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATAtt                        <  1:3162123/65‑1 (MQ=255)                                      ataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcc                     <  1:2821034/65‑1 (MQ=255)                                        aCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCtt                   >  1:1857781/1‑65 (MQ=255)                                         cGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCttt                  >  1:2388608/1‑65 (MQ=255)                                         cGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCttt                  >  1:2927793/1‑65 (MQ=255)                                         cGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCttt                  >  1:671398/1‑65 (MQ=255)                                          gTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttg                 <  1:3018096/65‑1 (MQ=255)                                          gTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttg                 <  1:2824637/65‑1 (MQ=255)                                           tCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgt                >  1:246270/1‑65 (MQ=255)                                           tCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTACCCTCTCACGATATTACCTttgt                >  1:1149359/1‑65 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTggc            <  1:33024/65‑2 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:1080299/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:1085561/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:431967/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:2740873/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:2690525/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:118076/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:2488242/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:700921/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:1376941/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:1624264/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:1758796/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:1969760/65‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:2818017/65‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:2679937/65‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:684879/65‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:569564/65‑1 (MQ=255)                                                  atatGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAtt         >  1:1965977/1‑65 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:1024418/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:478696/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:3037153/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:3032963/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:2953130/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:2662059/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:2302930/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:2096856/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:1297186/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:1108230/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTGTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:2315313/65‑1 (MQ=255)                                                    agagATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:1806489/62‑1 (MQ=255)                                                         taaaattaaaatTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGa  <  1:2825074/65‑1 (MQ=255)                                                         taaaattaaaatTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGa  <  1:2412749/65‑1 (MQ=255)                                                         taaaattaaaatTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGa  <  1:1210765/65‑1 (MQ=255)                                                             |                                                             AAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGA  >  NC_000913/803602‑803722

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