Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	613,456	C→T	R240R (CGC→CGT)	fes →	enterochelin esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	613,456	0	C	T	100.0%	168.4 / NA	48	R240R (CGC→CGT)	fes	enterochelin esterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/26);  total (22/26)

GCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTT  >  NC_000913/613321‑613592                                                                                                                                        |                                                                                                                                         gCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGccc                                                                                                                                      >  1:677486/1‑140 (MQ=255) gCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGccc                                                                                                                                      <  2:1051529/140‑1 (MQ=255)   tGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCAc                                                                                                                                    >  1:489931/1‑140 (MQ=255)      cAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGaa                                                                                                                                 <  2:677486/140‑1 (MQ=255)               tCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCgtgt                                                                                                                        >  2:426569/1‑140 (MQ=255)                   tGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATg                                                                                                                    >  1:187906/1‑140 (MQ=255)                   tGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATg                                                                                                                    >  2:349474/1‑140 (MQ=255)                   tGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATg                                                                                                                    >  2:880363/1‑140 (MQ=255)                        aaTTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGAt                                                                                                               >  1:960993/1‑140 (MQ=255)                             tGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTg                                                                                                          <  2:750483/140‑1 (MQ=255)                                gcgcAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCt                                                                                                       <  2:470996/140‑1 (MQ=255)                                gcgcAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCt                                                                                                       <  1:9229/140‑1 (MQ=255)                                 cgcAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTc                                                                                                      <  1:981783/140‑1 (MQ=255)                                  gcAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCg                                                                                                     <  1:349474/140‑1 (MQ=255)                                   cAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGc                                                                                                    <  2:489931/140‑1 (MQ=255)                                          tGCCCGTCTGGTCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAg                                                                                             <  2:333925/140‑1 (MQ=255)                                           gCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGc                                                                                            <  1:1082749/140‑1 (MQ=255)                                             ccGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCaa                                                                                          <  1:373493/140‑1 (MQ=255)                                              cGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAag                                                                                         <  1:560782/140‑1 (MQ=255)                                              cGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAag                                                                                         <  2:375106/140‑1 (MQ=255)                                               gTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGTCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAaga                                                                                        <  2:675401/140‑1 (MQ=255)                                                  tGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCaa                                                                                          >  2:598287/1‑135 (MQ=255)                                                  tGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCaa                                                                                          <  1:598287/135‑1 (MQ=255)                                                  tGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGtt                                                                                     <  1:104413/140‑1 (MQ=255)                                                  tGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGtt                                                                                     >  2:104413/1‑140 (MQ=255)                                                    gCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGttat                                                                                   <  2:187906/140‑1 (MQ=255)                                                    gCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGttat                                                                                   <  2:960993/140‑1 (MQ=255)                                                         tGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGc                                                                                            >  1:669085/1‑126 (MQ=255)                                                         tGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGc                                                                                            <  2:669085/126‑1 (MQ=255)                                                              aCTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTAccc                                                                              <  1:593079/135‑1 (MQ=255)                                                              aCTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTAccc                                                                              >  2:593079/1‑135 (MQ=255)                                                              aCTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTAccc                                                                              >  2:1085939/1‑135 (MQ=255)                                                              aCTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTAccc                                                                              <  1:1085939/135‑1 (MQ=255)                                                                          cATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGc                                                             <  2:766235/140‑1 (MQ=255)                                                                                        tcctccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATc                                               >  2:293558/1‑140 (MQ=255)                                                                                           tccCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTg                                            <  1:293558/140‑1 (MQ=255)                                                                                                   tgtATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGa                                                 <  1:617236/127‑1 (MQ=255)                                                                                                   tgtATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGa                                                 >  2:617236/1‑127 (MQ=255)                                                                                                   tgtATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGAt                                       <  2:896807/137‑1 (MQ=255)                                                                                                   tgtATGTGTTGATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGAt                                       >  1:896807/1‑137 (MQ=255)                                                                                                             gATCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTg                          >  1:748021/1‑140 (MQ=255)                                                                                                              aTCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGc                         >  1:610708/1‑140 (MQ=255)                                                                                                               tCGACGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGcc                        >  2:511813/1‑140 (MQ=255)                                                                                                                  aCGCTATCGACACCACGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCggg                     >  2:893350/1‑140 (MQ=255)                                                                                                                             accacGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtgg          >  1:631627/1‑140 (MQ=255)                                                                                                                              ccacGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggt         >  1:828952/1‑140 (MQ=255)                                                                                                                                acGCACCGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggtgg       >  1:446237/1‑140 (MQ=255)                                                                                                                                     ccGTGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCttt  <  2:748021/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                         GCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTT  >  NC_000913/613321‑613592

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