Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,246,382	T→G	I130L (ATC→CTC)	cirA ←	colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,246,382	0	T	G	100.0%	195.1 / NA	52	I130L (ATC→CTC)	cirA	colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (26/26);  total (26/26)

TCGCGATGTTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGATGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGCTGCTGTCCAGACCACGAATACTAACGCCCTTACGGTTATCCCCTTCGTT  >  NC_000913/2246234‑2246529                                                                                                                                                     |                                                                                                                                                    tcgcgatGTTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAggg                                                                                                                                                   >  2:163919/1‑151 (MQ=255)  cgcgatGTTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAggg                                                                                                                                                   >  1:2270/1‑150 (MQ=255)  cgcgatGTTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGa                                                                                                                                                  <  1:490781/151‑1 (MQ=255)      atGTTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCg                                                                                                                                              <  2:105412/151‑1 (MQ=255)       tGTTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCg                                                                                                                                              >  2:385828/1‑150 (MQ=255)        gTTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGa                                                                                                                                             <  1:153729/150‑1 (MQ=255)             tGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCggg                                                                                                                                        >  1:248195/1‑150 (MQ=255)              gAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGa                                                                                                                                       <  2:557116/150‑1 (MQ=255)                           gACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCg                                                                                                                         <  1:200557/151‑1 (MQ=255)                            acggtaacggtaCCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCg                                                                                                                         >  2:392464/1‑150 (MQ=255)                             cggtaacggtaCCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGaa                                                                                                                       <  1:805398/151‑1 (MQ=255)                             cggtaacggtaCCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGaa                                                                                                                       <  1:715415/151‑1 (MQ=255)                             cggtaacggtaCCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGa                                                                                                                        >  1:285949/1‑150 (MQ=255)                                    ggtaCCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATtg                                                                                                                >  1:790105/1‑151 (MQ=255)                                       aCCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTgg                                                                                                             >  2:14844/1‑151 (MQ=255)                                        cccGACCATTTCTGACCAATTTTTTTGGTGATGATATCCCCTCCACCGCCGAGCGCATGGGAACCGTACAGCGACGACATCTGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTgg                                                                                                             <  1:77993/150‑1 (MQ=255)                                           gACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCgg                                                                                                          >  1:103384/1‑150 (MQ=255)                                            aCCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGa                                                                                                         >  2:780619/1‑150 (MQ=255)                                              cATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGa                                                                                                      <  2:434614/151‑1 (MQ=255)                                                 ttCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACg                                                                                                    <  1:470956/150‑1 (MQ=255)                                                     gACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCa                                                                                                 <  1:49962/149‑1 (MQ=255)                                                          aTTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCggg                                                                                          >  2:19793/1‑151 (MQ=255)                                                          aTTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCggg                                                                                          >  2:103366/1‑151 (MQ=255)                                                                  ggTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCAcg                                                                                  <  1:66158/151‑1 (MQ=255)                                                                           atTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCg                                                                         >  2:785585/1‑151 (MQ=255)                                                                              cACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTc                                                                       <  2:475350/150‑1 (MQ=255)                                                                                  acacCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTcgac                                                                    <  2:904815/149‑1 (MQ=255)                                                                                   cacCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTAACGTcgacg                                                                   >  1:269777/1‑149 (MQ=255)                                                                                      cgccgAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAAt                                                              <  1:486904/151‑1 (MQ=255)                                                                                       gccgAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAAt                                                              >  1:741302/1‑150 (MQ=255)                                                                                          gAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAgg                                                          <  2:468831/151‑1 (MQ=255)                                                                                            gcgcATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAggg                                                         >  1:830729/1‑150 (MQ=255)                                                                                               cATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGtata                                                     <  2:686376/151‑1 (MQ=255)                                                                                                 tCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGgata                                                     >  2:118120/1‑149 (MQ=255)                                                                                                 tCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAgc                                                   >  2:847657/1‑151 (MQ=255)                                                                                                   ggAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAgctg                                                 >  2:362871/1‑151 (MQ=255)                                                                                                     aaCCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAgctgct                                               >  1:291700/1‑151 (MQ=255)                                                                                                     aaCCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGGGCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAgctgc                                                <  1:300120/150‑1 (MQ=255)                                                                                                      aCCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGctgctg                                              <  1:322680/151‑1 (MQ=255)                                                                                                         gTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGCTGCTGTc                                            >  2:214548/1‑150 (MQ=255)                                                                                                           aCAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGCTGCTGTcc                                           >  2:572249/1‑149 (MQ=255)                                                                                                           aCAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGGATAGCTGCTGTcc                                           >  1:706725/1‑149 (MQ=255)                                                                                                             cgCGACGACATCGGCCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGGTCAGATCGAAATCATTTTGGCGGAAGACGGCATTGCGGGAGTTCACCGGTTTACCGTCGACGCAACTCAGGGTATAGCGGCTGTCCAGa                                        >  1:623368/2‑150 (MQ=255)                                                                                                               cgacgaCATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGCTGCTGTCCAGAcc                                      >  2:270106/1‑150 (MQ=255)                                                                                                                     cATCGGGCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGCTGCTGTCCAGACCACGAATa                               <  2:418095/151‑1 (MQ=255)                                                                                                                         gggCCACGGACCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGCTGCTGTCCAGACCACGAATACTaa                           <  1:381531/151‑1 (MQ=255)                                                                                                                                  aCCACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGCTGCTGTCCAGACCACGAATACTAACGCCCTTAc                  <  2:2270/151‑1 (MQ=255)                                                                                                                                    cACTTCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGCTGCTGTCCAGACCACGAATACTAACGCCCTTACgg                >  1:547524/1‑151 (MQ=255)                                                                                                                                        tCAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGCTGCTGTCCAGACCACGAATACTAACGCCCTTACGGtt              <  1:665235/149‑1 (MQ=255)                                                                                                                                         cAATACGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGCTGCTGTCCAGACCACGAATACTAACGCCCTTACGGTTAt            <  2:72727/150‑1 (MQ=255)                                                                                                                                             aCGTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGCTGCTGTCCAGACCACGAATACTAACGCCCTTACGGTTATCCCCt       <  1:226356/151‑1 (MQ=255)                                                                                                                                               gTTCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGCTGCTGTCCAGACCACGAATACTAACGCCCTTACGGTTATCCCCttc     >  2:586744/1‑151 (MQ=255)                                                                                                                                                 tCGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGCTGCTGTCCAGACCACGAATACTAACGCCCTTACGGTTATCCCCttcg    >  1:581105/1‑150 (MQ=255)                                                                                                                                                  cGAGGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGCTGCTGTCCAGACCACGAATACTAACGCCCTTACGGTTATCCCCttcgtt  <  2:248201/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TCGCGATGTTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGTTCGATGGAATCGACCGGGATCCAGTTCAGATCGAAATCATTGTGGCGGAAGACGGCATTGCGGGAGTTCACGCGTTTACCGTCGACGAGAATCAGGGTATAGCTGCTGTCCAGACCACGAATACTAACGCCCTTACGGTTATCCCCTTCGTT  >  NC_000913/2246234‑2246529

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