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

GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGT  >  NC_000913/803605‑803717                                                          |                                                        ggCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGAtaaaattaaaat                                                  <  1:254131/65‑1 (MQ=255)      tATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGt                                             <  1:262610/65‑1 (MQ=255)      tATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGt                                             <  1:383478/65‑1 (MQ=255)         tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTg                                          >  1:1941246/1‑65 (MQ=255)         tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTg                                          >  1:2904675/1‑65 (MQ=255)            cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                       >  1:2088887/1‑65 (MQ=255)            cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                       >  1:2030559/1‑65 (MQ=255)            cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                       >  1:685872/1‑65 (MQ=255)            cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                       >  1:2363927/1‑65 (MQ=255)            cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                       >  1:434498/1‑65 (MQ=255)              gtgcTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                     <  1:1461639/62‑1 (MQ=255)              ggTCTATTCTCTCTTCGCGTAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                     <  1:353668/65‑1 (MQ=255)              ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                     >  1:984418/1‑65 (MQ=255)              ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                     >  1:2191017/1‑65 (MQ=255)              ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                     >  1:2825922/1‑65 (MQ=255)              ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                     >  1:3066126/1‑65 (MQ=255)              ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                     >  1:1579926/1‑65 (MQ=255)              ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                     >  1:12823/1‑65 (MQ=255)                tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAACATTAAAATTAAGTTTGATTCAtt                                   >  1:2214777/1‑65 (MQ=255)                tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                   >  1:2921558/1‑65 (MQ=255)                tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                   >  1:408665/1‑65 (MQ=255)                   aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                <  1:1877685/65‑1 (MQ=255)                   aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                <  1:685787/65‑1 (MQ=255)                    ttCTCGCGTCGCATAATACGTCCGGAAAGATGGATAAAATTAAAATTAAGTTTGATTCATTAGcc                               <  1:2125095/65‑1 (MQ=255)                     tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:3267530/1‑65 (MQ=255)                     tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:387635/1‑65 (MQ=255)                     tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:639947/1‑65 (MQ=255)                     tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:2381801/1‑65 (MQ=255)                     tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:704338/1‑65 (MQ=255)                     tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:2889878/1‑65 (MQ=255)                     tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:1535873/1‑65 (MQ=255)                     tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:1810888/1‑65 (MQ=255)                     tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:1496353/1‑65 (MQ=255)                     tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:1924796/1‑65 (MQ=255)                     tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:1421108/1‑65 (MQ=255)                      ctctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctc                             <  1:3084457/65‑1 (MQ=255)                      ctctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctc                             <  1:1486397/65‑1 (MQ=255)                        ctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                           >  1:1982310/1‑65 (MQ=255)                          ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCAc                         >  1:44932/1‑65 (MQ=255)                           ttCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGGTTCATTAGCCTCTCACg                        <  1:1828970/65‑1 (MQ=255)                           ttCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACg                        <  1:856987/65‑1 (MQ=255)                            tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                       <  1:73594/65‑1 (MQ=255)                            tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                       <  1:189142/65‑1 (MQ=255)                            gcGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                       <  1:74259/64‑1 (MQ=255)                             cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      <  1:424635/65‑1 (MQ=255)                             cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      <  1:1498127/65‑1 (MQ=255)                             cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      <  1:2043934/65‑1 (MQ=255)                                  aataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAc                 <  1:929860/65‑1 (MQ=255)                                   ataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcc                >  1:2058222/1‑65 (MQ=255)                                   ataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcc                <  1:300192/65‑1 (MQ=255)                                   ataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcc                <  1:766274/65‑1 (MQ=255)                                     aCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcccc              >  1:3296889/1‑63 (MQ=255)                                     aCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCtt              >  1:3081483/1‑65 (MQ=255)                                     aCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCtt              >  1:3066988/1‑65 (MQ=255)                                     aCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCtt              >  1:760348/1‑65 (MQ=255)                                      cGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCttt             >  1:1172110/1‑65 (MQ=255)                                       gTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttg            <  1:221186/65‑1 (MQ=255)                                       gTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttg            <  1:2286067/65‑1 (MQ=255)                                        tCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgt           >  1:310574/1‑65 (MQ=255)                                            ggaaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:3196034/63‑1 (MQ=255)                                            gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:1784420/65‑1 (MQ=255)                                            gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:3270591/65‑1 (MQ=255)                                            gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:3329045/65‑1 (MQ=255)                                            gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:1597337/65‑1 (MQ=255)                                            gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:839480/65‑1 (MQ=255)                                            gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:3202501/65‑1 (MQ=255)                                            gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:794355/65‑1 (MQ=255)                                            gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:3132476/65‑1 (MQ=255)                                            gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:2299940/65‑1 (MQ=255)                                            gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:2475985/65‑1 (MQ=255)                                            gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:1309766/65‑1 (MQ=255)                                              aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt     <  1:2167621/65‑1 (MQ=255)                                              aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt     <  1:1827031/65‑1 (MQ=255)                                              aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt     <  1:1559019/65‑1 (MQ=255)                                              aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTATTGGTAt     <  1:3211587/65‑1 (MQ=255)                                                 atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTGCCTTTGTTGGTATTGt  <  1:1977514/65‑1 (MQ=255)                                                 atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:1018969/65‑1 (MQ=255)                                                 atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:1332851/65‑1 (MQ=255)                                                 atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  >  1:505839/1‑65 (MQ=255)                                                 atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:464523/65‑1 (MQ=255)                                                 atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:1996862/65‑1 (MQ=255)                                                 atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:400873/65‑1 (MQ=255)                                                 atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:3311119/65‑1 (MQ=255)                                                 atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:885501/65‑1 (MQ=255)                                                 atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:3244521/65‑1 (MQ=255)                                                 atGGATAAAATTAAAATTAAGCTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:2525092/65‑1 (MQ=255)                                                          |                                                        GGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGT  >  NC_000913/803605‑803717

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