Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,246,914	T→C	intergenic (‑145/+149)	cirA ← / ← lysP	colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,246,914	0	T	C	100.0%	192.4 / NA	53	intergenic (‑145/+149)	cirA/lysP	colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/23);  total (30/23)

ATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAA  >  NC_000913/2246768‑2247060                                                                                                                                                   |                                                                                                                                                   aTTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGc                                                                                                                                                >  1:118667/1‑151 (MQ=255) aTTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGc                                                                                                                                                >  1:462881/1‑151 (MQ=255)     cATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCttt                                                                                                                                            <  1:812623/151‑1 (MQ=255)      atatCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGAC‑TATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCttttt                                                                                                                                          >  1:322445/1‑151 (MQ=255)        atCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCtttttt                                                                                                                                         >  2:254226/1‑151 (MQ=255)        atCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCtttttt                                                                                                                                         >  1:614355/1‑151 (MQ=255)           tcCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGAt                                                                                                                                      >  2:391609/1‑151 (MQ=255)             cATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTg                                                                                                                                    >  1:617522/1‑151 (MQ=255)              aTGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTgc                                                                                                                                   >  2:436191/1‑151 (MQ=255)              aTGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTgc                                                                                                                                   >  1:477865/1‑151 (MQ=255)                 aGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTgcggc                                                                                                                                <  1:840908/151‑1 (MQ=255)                     aaCTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTaa                                                                                                                            <  2:305162/151‑1 (MQ=255)                         aCGAAAATAAAATGGGTTATCGCTCACATCTTCTTCGCGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTAGGGAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGAt                                                                                                                        <  1:668108/151‑1 (MQ=255)                            aaaataaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAAt                                                                                                                     >  2:445961/1‑151 (MQ=255)                            aaaataaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAAt                                                                                                                     >  1:147267/1‑151 (MQ=255)                             aaataaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATg                                                                                                                    >  1:18003/1‑151 (MQ=255)                              aataaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGc                                                                                                                   <  1:69852/151‑1 (MQ=255)                              aataaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGc                                                                                                                   <  2:262641/151‑1 (MQ=255)                                 aaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCaa                                                                                                                 >  1:324592/1‑150 (MQ=255)                                      gggTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGAt                                                                                                            >  2:689671/1‑150 (MQ=255)                                           aTCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAAt                                                                                                      <  2:518552/151‑1 (MQ=255)                                              gCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTAt                                                                                                   >  1:132017/1‑151 (MQ=255)                                                tCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATc                                                                                                  >  2:120631/1‑150 (MQ=255)                                                          ttcACGTCACCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAGATGt                                                                                       >  1:713192/1‑151 (MQ=255)                                                             aCGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTa                                                                                     >  2:257334/1‑150 (MQ=255)                                                                 ccccTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAAttt                                                                                >  1:212222/1‑151 (MQ=255)                                                                  cccTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAAtttt                                                                               <  2:194658/151‑1 (MQ=255)                                                                   ccTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAAtttt                                                                               >  1:314563/1‑150 (MQ=255)                                                                       gCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTatat                                                                           <  1:138695/150‑1 (MQ=255)                                                                          tACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGc                                                                       >  1:210768/1‑151 (MQ=255)                                                                           aCGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGc                                                                       >  1:581141/1‑150 (MQ=255)                                                                            cGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCAt                                                                     <  1:557917/151‑1 (MQ=255)                                                                             ggCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATg                                                                    <  1:609178/151‑1 (MQ=255)                                                                             ggCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATg                                                                    >  1:872734/1‑151 (MQ=255)                                                                              gcagcaTCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGa                                                                   <  2:456970/151‑1 (MQ=255)                                                                                 gcaTCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGAttt                                                                <  1:380251/151‑1 (MQ=255)                                                                                       cgGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAgg                                                          <  2:669587/151‑1 (MQ=255)                                                                                             aCTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATg                                                    <  1:106000/151‑1 (MQ=255)                                                                                              cTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATgg                                                   >  1:408975/1‑151 (MQ=255)                                                                                                        aTGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGaaaa                                          <  2:694875/150‑1 (MQ=255)                                                                                                         tGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGaaaa                                          >  2:538788/1‑149 (MQ=255)                                                                                                         tGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGaaaa                                          >  2:510942/1‑149 (MQ=255)                                                                                                                 atcatcCCAGAAAGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCg                                <  2:189187/151‑1 (MQ=255)                                                                                                                           aaaGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAggg                      >  2:453033/1‑151 (MQ=255)                                                                                                                           aaaGCCGGAACAGCCTCTGATTTGCCAGCCTTTTTCGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAggg                      <  1:721066/151‑1 (MQ=255)                                                                                                                              gCCGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCtt                   <  1:510992/151‑1 (MQ=255)                                                                                                                               ccGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCtt                   >  2:304270/1‑150 (MQ=255)                                                                                                                               ccGGAACAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTa                  <  1:139490/151‑1 (MQ=255)                                                                                                                                     cAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaagg            <  2:572329/151‑1 (MQ=255)                                                                                                                                     cAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaagg            <  2:581091/151‑1 (MQ=255)                                                                                                                                     cAGCCTCTGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaagg            >  2:794752/1‑151 (MQ=255)                                                                                                                                            tGATTTGCCAGCCTTTTTTGATTGCGGCGTAACGCTAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAggg       >  2:739373/1‑149 (MQ=255)                                                                                                                                                ttGCCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGaa  <  1:16876/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   ATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAA  >  NC_000913/2246768‑2247060

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