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

AGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGA  >  NC_000913/803604‑803722                                                           |                                                             aGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTaaaa                                                        >  1:3371661/1‑65 (MQ=255)     aCTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAaattaaaattaa                                                    >  1:1789442/1‑65 (MQ=255)     aCTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAaattaaaattaa                                                    <  1:720571/65‑1 (MQ=255)      cTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAg                                                   <  1:155280/65‑1 (MQ=255)        aTTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGtt                                                 >  1:2257736/1‑65 (MQ=255)          tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTg                                               >  1:229143/1‑65 (MQ=255)          tGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTg                                               >  1:818976/1‑65 (MQ=255)             cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                            >  1:2768983/1‑65 (MQ=255)             cTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                            >  1:170750/1‑65 (MQ=255)             cTGGTCTATTCTCTCTTCGCATAATACATCCGGAAATATGGATAAAATTAAAATTAAGTTTGAtt                                            >  1:2384405/1‑65 (MQ=255)              tGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTc                                           <  1:2947883/65‑1 (MQ=255)              tGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTc                                           <  1:2724367/65‑1 (MQ=255)              tGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTc                                           <  1:1753627/65‑1 (MQ=255)               ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                          <  1:3419686/65‑1 (MQ=255)               ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                          >  1:315573/1‑65 (MQ=255)               ggTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCa                                          >  1:3780117/1‑65 (MQ=255)                 tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                        >  1:2504626/1‑65 (MQ=255)                 tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                        >  1:3386759/1‑65 (MQ=255)                 tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                        >  1:464643/1‑65 (MQ=255)                 tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                        >  1:3555377/1‑65 (MQ=255)                 tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                        >  1:2719019/1‑65 (MQ=255)                 tCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCAtt                                        >  1:487861/1‑65 (MQ=255)                     ttCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGcc                                    <  1:2540448/65‑1 (MQ=255)                      tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:3209696/1‑65 (MQ=255)                      tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:2115972/1‑65 (MQ=255)                      tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:3194569/1‑65 (MQ=255)                      tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:306410/1‑65 (MQ=255)                      tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:103578/1‑65 (MQ=255)                      tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:3034736/1‑65 (MQ=255)                      tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:297133/1‑65 (MQ=255)                      tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   <  1:2132989/65‑1 (MQ=255)                      tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:3688981/1‑65 (MQ=255)                      tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:1672876/1‑65 (MQ=255)                      tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:61037/1‑65 (MQ=255)                      tctctcTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCct                                   >  1:925122/1‑65 (MQ=255)                        tctctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCctct                                 <  1:2358003/65‑1 (MQ=255)                           ctTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAACTTTGATTCATTAGCCTCTCAc                              >  1:638494/1‑65 (MQ=255)                            ttCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACg                             <  1:1884030/65‑1 (MQ=255)                            ttCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACg                             <  1:658276/65‑1 (MQ=255)                             tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                            <  1:3864680/65‑1 (MQ=255)                             tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                            <  1:1430388/65‑1 (MQ=255)                             tCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGa                            <  1:989/65‑1 (MQ=255)                              cGCATAATACGTCCGGAAATATTGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                           <  1:3341234/65‑1 (MQ=255)                              cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                           <  1:3058237/65‑1 (MQ=255)                              cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                           <  1:3520783/65‑1 (MQ=255)                              cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                           <  1:559935/65‑1 (MQ=255)                              cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGat                           <  1:2056389/65‑1 (MQ=255)                                  taataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTa                       <  1:1688484/65‑1 (MQ=255)                                   aataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAc                      <  1:1650385/65‑1 (MQ=255)                                   aataCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTAc                      <  1:3816825/65‑1 (MQ=255)                                      aCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCtt                   >  1:555166/1‑65 (MQ=255)                                      aCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCtt                   <  1:589338/65‑1 (MQ=255)                                      aCGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCtt                   >  1:3326860/1‑65 (MQ=255)                                       cGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCttt                  >  1:45873/1‑65 (MQ=255)                                       cGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCttt                  >  1:3375367/1‑65 (MQ=255)                                         tCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgt                >  1:2986234/1‑65 (MQ=255)                                         tCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTttgt                >  1:3228407/1‑65 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:2019228/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:3565408/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:187480/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:1654907/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:3452529/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:593090/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:338388/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:1606163/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:1472484/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:2878011/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:3294966/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:325604/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:2797202/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:2700098/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:699162/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:1425580/65‑1 (MQ=255)                                             gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:958003/65‑1 (MQ=255)                                             agaaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGt            <  1:1927008/63‑1 (MQ=255)                                               aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:2043530/65‑1 (MQ=255)                                               aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:3003752/65‑1 (MQ=255)                                               aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:3099210/65‑1 (MQ=255)                                               aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:3219404/65‑1 (MQ=255)                                               aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:3246241/65‑1 (MQ=255)                                               aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:3306769/65‑1 (MQ=255)                                               aaTATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:2301822/65‑1 (MQ=255)                                               aaTATGGATAAAATTAAAATTAAGTTGGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:1072382/65‑1 (MQ=255)                                               aaTATGGATAAAATTAAAATGAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTAt          <  1:2727255/65‑1 (MQ=255)                                                 tatGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTg        >  1:2121775/1‑65 (MQ=255)                                                  atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:3792704/65‑1 (MQ=255)                                                  atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:1961487/65‑1 (MQ=255)                                                  atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:3514265/65‑1 (MQ=255)                                                  atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:2179598/65‑1 (MQ=255)                                                  atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:3352459/65‑1 (MQ=255)                                                  atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:2331549/65‑1 (MQ=255)                                                  atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:711302/65‑1 (MQ=255)                                                  atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:3172001/65‑1 (MQ=255)                                                  atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:1249922/65‑1 (MQ=255)                                                  atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:2727713/65‑1 (MQ=255)                                                  atGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGt       <  1:2890168/65‑1 (MQ=255)                                                      aTAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAg   >  1:1297705/1‑65 (MQ=255)                                                       taaaattaaaatTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGa  <  1:3338555/65‑1 (MQ=255)                                                           |                                                             AGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGA  >  NC_000913/803604‑803722

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