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

AAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGT  >  NC_000913/803602‑803717                                                             |                                                        aaaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGAtaaaattaaaat                                                  >  1:2211627/1‑68 (MQ=255)  aaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATaaaattaaaatt                                                 >  1:1736201/1‑68 (MQ=255)      cACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGt                                             <  1:623316/68‑1 (MQ=255)      cACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGt                                             <  1:2991642/68‑1 (MQ=255)          aTTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGa                                         >  1:640636/1‑68 (MQ=255)          aTTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGa                                         <  1:1652085/68‑1 (MQ=255)             gCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTc                                      <  1:3074422/68‑1 (MQ=255)               cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAt                                    >  1:2638393/1‑68 (MQ=255)               cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAt                                    >  1:1792653/1‑68 (MQ=255)                 ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTa                                  >  1:301172/1‑68 (MQ=255)                 ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTa                                  >  1:233736/1‑68 (MQ=255)                  gTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAg                                 <  1:1519428/68‑1 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                >  1:325343/1‑68 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                >  1:3255825/1‑68 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                <  1:1909601/68‑1 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                >  1:641158/1‑68 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                <  1:1126302/68‑1 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                >  1:2456400/1‑68 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                >  1:454953/1‑68 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                >  1:330006/1‑68 (MQ=255)                   tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGc                                >  1:2855940/1‑68 (MQ=255)                    cTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGcc                               >  1:308419/1‑68 (MQ=255)                    cTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGcc                               <  1:1991305/68‑1 (MQ=255)                    cTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGcc                               <  1:1171208/68‑1 (MQ=255)                     tATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                              <  1:510817/68‑1 (MQ=255)                        tctctcTTCGCATAATACGTCCGGCAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctctg                           >  1:2973285/1‑67 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                           >  1:2553504/1‑68 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                           >  1:2040023/1‑68 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                           >  1:3060309/1‑68 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                           >  1:2283894/1‑68 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                           >  1:353351/1‑68 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                           >  1:755378/1‑68 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                           >  1:911315/1‑68 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                           >  1:1680622/1‑68 (MQ=255)                        tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                           >  1:1700914/1‑68 (MQ=255)                           ctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACg                        <  1:173087/68‑1 (MQ=255)                            tctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                       <  1:2065410/68‑1 (MQ=255)                            tctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                       <  1:829121/68‑1 (MQ=255)                            tctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                       <  1:1356179/68‑1 (MQ=255)                            tctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                       <  1:1681533/68‑1 (MQ=255)                            tctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                       <  1:2563161/68‑1 (MQ=255)                             ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      <  1:3511999/68‑1 (MQ=255)                             ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      >  1:90216/1‑68 (MQ=255)                             ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      <  1:2987791/68‑1 (MQ=255)                             ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      >  1:3391410/1‑68 (MQ=255)                             ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      <  1:1796755/68‑1 (MQ=255)                             ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      <  1:2469369/68‑1 (MQ=255)                             ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      <  1:484394/68‑1 (MQ=255)                             ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      <  1:1983860/68‑1 (MQ=255)                             ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                      >  1:2304233/1‑68 (MQ=255)                                   ataataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcc                <  1:1455984/68‑1 (MQ=255)                                   ataataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcc                <  1:1050568/68‑1 (MQ=255)                                   ataataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAcc                <  1:1802294/68‑1 (MQ=255)                                       taCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttg            <  1:1083486/68‑1 (MQ=255)                                        aCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgt           >  1:1333605/1‑68 (MQ=255)                                        aCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgt           >  1:1097931/1‑68 (MQ=255)                                         cGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgtt          >  1:3400985/1‑68 (MQ=255)                                         cGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgtt          >  1:1236575/1‑68 (MQ=255)                                         cGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgtt          >  1:1162089/1‑68 (MQ=255)                                            ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:709943/68‑1 (MQ=255)                                            ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:621469/68‑1 (MQ=255)                                            ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:954745/68‑1 (MQ=255)                                            ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:311816/68‑1 (MQ=255)                                            ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:1030963/68‑1 (MQ=255)                                            ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:280995/68‑1 (MQ=255)                                            ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:1532497/68‑1 (MQ=255)                                            ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:2589094/68‑1 (MQ=255)                                            ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:2306269/68‑1 (MQ=255)                                            ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:2228306/68‑1 (MQ=255)                                            ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:2150442/68‑1 (MQ=255)                                            ccGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt       <  1:2025016/68‑1 (MQ=255)                                              ggAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt     <  1:289495/68‑1 (MQ=255)                                              ggAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt     <  1:840140/68‑1 (MQ=255)                                              ggAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt     <  1:732203/68‑1 (MQ=255)                                              ggAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt     <  1:1830842/68‑1 (MQ=255)                                              cgAAATATTGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt     <  1:10178/67‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTga  <  1:1056070/68‑2 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:690151/68‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:1784493/68‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:2385537/68‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:2769926/68‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:1343814/68‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:3410343/68‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:2830640/68‑1 (MQ=255)                                                 aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt  <  1:3013462/68‑1 (MQ=255)                                                             |                                                        AAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGT  >  NC_000913/803602‑803717

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