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

AAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGT  >  NC_000913/803602‑803717                                                             |                                                        aaaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTaa                                                     >  1:21111/1‑65 (MQ=255)  aaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTaaa                                                    >  1:1199472/1‑65 (MQ=255)  aaGGCACTATTGCCTGGTCTATTCTCTCTTAGCATAATACGTCCGGAAATATGGATAAAATTaaa                                                    >  1:3074458/1‑65 (MQ=255)    ggCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGAtaaaattaaaat                                                  >  1:1719827/1‑65 (MQ=255)       aCTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAaattaaaattaa                                               >  1:361071/1‑65 (MQ=255)       aCTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAaattaaaattaa                                               >  1:874008/1‑65 (MQ=255)         cagtGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGt                                             <  1:2969070/62‑1 (MQ=255)            tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTg                                          >  1:825048/1‑65 (MQ=255)              ccTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAt                                        >  1:2923436/1‑65 (MQ=255)               cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                       >  1:580059/1‑65 (MQ=255)               cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                       >  1:5252/1‑65 (MQ=255)               cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                       >  1:494743/1‑65 (MQ=255)                 ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                     >  1:3177040/1‑65 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                   >  1:2986893/1‑65 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                   >  1:183385/1‑65 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                   >  1:764163/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:2411532/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:496212/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:702194/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:2607066/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:2574916/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:3153740/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:2535038/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:2344248/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:863998/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:1394639/1‑65 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:91963/1‑65 (MQ=255)                        tctctcTTCACATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              >  1:989731/1‑65 (MQ=255)                         ctctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctc                             <  1:1416270/65‑1 (MQ=255)                           ctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                           >  1:2261502/1‑65 (MQ=255)                              ttCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACg                        <  1:2141102/65‑1 (MQ=255)                               tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                       <  1:2015067/65‑1 (MQ=255)                               tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                       <  1:1954981/65‑1 (MQ=255)                                cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      <  1:60892/65‑1 (MQ=255)                                cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      <  1:460164/65‑1 (MQ=255)                                cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      <  1:925254/65‑1 (MQ=255)                                cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      <  1:1893513/65‑1 (MQ=255)                                cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      <  1:2546882/65‑1 (MQ=255)                                cGCATAATACGTCCGGAAATATGGATAAAATTAAAATGAAGTTTGATTCATTAGCCTCTCACGat                      <  1:847343/65‑1 (MQ=255)                                  cATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGatat                    >  1:2781619/1‑65 (MQ=255)                                     aataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAc                 <  1:3214870/65‑1 (MQ=255)                                      ataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcc                >  1:1334926/1‑65 (MQ=255)                                      ataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcc                <  1:2739010/65‑1 (MQ=255)                                      ataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcc                <  1:2425050/65‑1 (MQ=255)                                       taCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCt               <  1:396803/65‑1 (MQ=255)                                       taCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTAATTAGCCTCTCACGATATTACCt               <  1:642692/65‑1 (MQ=255)                                        aCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCtc              >  1:204619/1‑64 (MQ=255)                                         cGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTCGCCTCTCACGATATTACCttt             >  1:26067/1‑65 (MQ=255)                                         cGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCttt             >  1:2093613/1‑65 (MQ=255)                                          gTCCGGAAATATGGATAAAATTAACATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttg            <  1:2809222/65‑1 (MQ=255)                                           tCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgt           >  1:1528188/1‑65 (MQ=255)                                           tCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgt           >  1:1596924/1‑65 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:1490221/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:1477948/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:1650823/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:172104/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:1796128/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:1050986/65‑1 (MQ=255)                                               gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:128202/65‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt     <  1:2696259/65‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt     <  1:519690/65‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt     <  1:2724734/65‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt     <  1:117940/65‑1 (MQ=255)                                                   tatGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTg   <  1:2251152/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTGCCTTTGTTGGTATTGt  <  1:2839887/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:436720/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:3138476/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:2850434/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:2568319/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:2279269/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:873976/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:1945485/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:876906/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:902748/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:1905485/65‑1 (MQ=255)                                                    atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  >  1:1844902/1‑65 (MQ=255)                                                             |                                                        AAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGT  >  NC_000913/803602‑803717

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