Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	609,496	A→C	L33F (TTA→TTC)	ybdZ →	enterobactin biosynthesis protein YbdZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	609,496	0	A	C	100.0%	148.4 / NA	40	L33F (TTA→TTC)	ybdZ	enterobactin biosynthesis protein YbdZ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (18/22);  total (18/22)

CTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAG  >  NZ_CP009273/609363‑609621                                                                                                                                      |                                                                                                                              ctctGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGc                                                                                                                          >  1:290221/1‑139 (MQ=255)      gCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCt                                                                                                                     >  1:530600/1‑139 (MQ=255)          ccACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGa                                                                                                                 <  1:241189/139‑1 (MQ=255)           cACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGAc                                                                                                                <  2:474510/139‑1 (MQ=255)              ttttCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACAtt                                                                                                             >  2:554738/1‑139 (MQ=255)                          ggAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCg                                                                                                 <  2:378157/139‑1 (MQ=255)                          ggAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCg                                                                                                 <  2:530600/139‑1 (MQ=255)                                  atatGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTcaca                                                                                         <  2:435365/139‑1 (MQ=255)                                     tGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAgg                                                                                                                       <  1:4816/106‑1 (MQ=255)                                     tGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAgg                                                                                                                       >  2:4816/1‑106 (MQ=255)                                           tCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTg                                                                                >  2:325212/1‑139 (MQ=255)                                                       tCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCt                                                                    >  1:15662/1‑139 (MQ=255)                                                       tCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCt                                                                    >  2:313098/1‑139 (MQ=255)                                                        cGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCt                                                                    <  1:554738/138‑1 (MQ=255)                                                           tgatCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTgg                                                                  >  1:12608/1‑137 (MQ=255)                                                                cGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCa                                                           <  1:650835/139‑1 (MQ=255)                                                                cGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCa                                                           <  2:54339/139‑1 (MQ=255)                                                                cGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCa                                                           <  1:430072/139‑1 (MQ=255)                                                                cGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCa                                                           <  2:326781/139‑1 (MQ=255)                                                                cGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCa                                                           <  2:668642/139‑1 (MQ=255)                                                                      gAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCg                                                     <  1:55574/139‑1 (MQ=255)                                                                      gAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCg                                                     <  2:457114/139‑1 (MQ=255)                                                                         cGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTcaca                                                                                         <  2:228617/100‑1 (MQ=255)                                                                         cGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTcaca                                                                                         >  1:228617/1‑100 (MQ=255)                                                                                     tGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACcgac                                      >  1:214501/1‑139 (MQ=255)                                                                                     tGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATACGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTAATCTGACACcgac                                      >  2:433430/1‑139 (MQ=255)                                                                                               gcAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTAcc                            <  2:172871/139‑1 (MQ=255)                                                                                               gcAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTAcc                            <  2:15662/139‑1 (MQ=255)                                                                                                        aTTCAGTCTGTGGCCGCAACAATGCGTCTTCCCgg                                                                                                                           >  1:588685/1‑35 (MQ=255)                                                                                                        aTTCAGTCTGTGGCCGCAACAATGCGTCTTCCCgg                                                                                                                           <  2:588685/35‑1 (MQ=255)                                                                                                        aTTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCag                   >  1:71629/1‑139 (MQ=255)                                                                                                        aTTCAGTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCag                   <  1:675611/139‑1 (MQ=255)                                                                                                            aGTCTGTGGCCGCAACAAGGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCaggagg               <  2:38576/139‑1 (MQ=255)                                                                                                             gTCTGTGGCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCaggagg               >  2:187371/1‑138 (MQ=255)                                                                                                                   ggCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTAcc                            <  2:669015/119‑1 (MQ=255)                                                                                                                   ggCCGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAAt         >  1:727958/1‑138 (MQ=255)                                                                                                                   ggCCGCAACAATGCGTCTTCCCGGCAGGCTGGCACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTAcc                            >  1:669015/1‑119 (MQ=255)                                                                                                                      cGCAACAATGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGc     <  2:615578/139‑1 (MQ=255)                                                                                                                        caacaaTGCGTCTTCCCGGCAGGCTGGGACAATGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGc     >  2:766274/1‑137 (MQ=255)                                                                                                                          acaaTGCGTCTTCCCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGCGGCACAATGAGCCAg  >  2:494824/1‑138 (MQ=255)                                                                                                                                      |                                                                                                                              CTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGGCAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAG  >  NZ_CP009273/609363‑609621

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑