Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,246,242	T→C	E176E (GAA→GAG)	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,242	0	T	C	100.0%	225.9 / NA	63	E176E (GAA→GAG)	cirA	colicin IA outer membrane receptor and translocator, ferric iron‑catecholate transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/29);  total (34/29)

TTTGCGGGTCATCCTTTTCACGTTTTGCCAGGCTGCCGTAAGCTTTCATTCCCAGCACACCATCAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGTTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGT  >  NC_000913/2246121‑2246377                                                                                                                          |                                                                                                                                        tttGCGGGTCATCCTTTTCACGTTTTGCCAGGCTGCCGTAAGCTTTCATTCCCAGCACACCATCAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATc                                                                                                                        >  2:452470/1‑139 (MQ=255) tttGCGGGTCATCCTTTTCACGTTTTGCCAGGCTGCCGTAAGCTTTCATTCCCAGCACACCAGCAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATTGTGGTATc                                                                                                                        <  1:34733/139‑1 (MQ=255)     cGGGTCATCCTTTTCACGTTTTGCCAGGCTGCCGTAAGCTTTCATTCCCAGCACACCATCAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGacg                                                                                                                    <  1:1544534/139‑1 (MQ=255)               ttttCACGTTTTGCCAGGCTGCCGTAAGCTTTCATTCCCAGCACACCATCAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTAc                                                                                                          <  2:909755/139‑1 (MQ=255)                    aCGTTTTGCCAGGCTGCCGTAAGCTTTCATTCCCAGCACACCATCAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGAc                                                                                                     <  2:7955/139‑1 (MQ=255)                       ttttGCCAGGCTGCCGTAAGCTTTCATTCCCAGCACACCATCAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCAt                                                                                                  >  2:1409811/1‑139 (MQ=255)                                              tCATTCCCAGCACACCATCAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATc                                                                                                                        <  1:448068/94‑1 (MQ=255)                                              tCATTCCCAGCACACCATCAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATc                                                                                                                        >  2:448068/1‑94 (MQ=255)                                               cATTCCCAGCACACCATCCATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGtgatg                                                                          <  2:1084006/139‑1 (MQ=255)                                               cATTCCCAGCACACCATCAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGtgatg                                                                          <  1:1409811/139‑1 (MQ=255)                                               cATTCCCAGCACACCATCAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGtgatg                                                                          <  1:199852/139‑1 (MQ=255)                                                        cacaCCATCAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACt                                                                 >  2:1004106/1‑139 (MQ=255)                                                             cATCAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGAc                                                                                                     <  1:411809/98‑1 (MQ=255)                                                             cATCAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGAc                                                                                                     >  2:411809/1‑98 (MQ=255)                                                                cAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACAccgcc                                                         >  1:290394/1‑139 (MQ=255)                                                                   ttAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAg                                                      <  1:200395/139‑1 (MQ=255)                                                                       tGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGa                                                                                                                      >  1:1430780/1‑71 (MQ=255)                                                                       tGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGa                                                                                                                      <  2:1430780/71‑1 (MQ=255)                                                                       tGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGtg                                                                             <  2:385185/112‑1 (MQ=255)                                                                       tGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGtg                                                                             >  1:385185/1‑112 (MQ=255)                                                                        ggTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCAt                                                 <  1:1465197/139‑1 (MQ=255)                                                                         gTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTAc                                                                                                          <  1:1532579/81‑1 (MQ=255)                                                                         gTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTAc                                                                                                          >  2:1532579/1‑81 (MQ=255)                                                                         gTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATAtt                                                                     >  2:680488/1‑118 (MQ=255)                                                                         gTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATAtt                                                                     <  1:680488/118‑1 (MQ=255)                                                                          tCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGacg                                                                                                                    >  1:1238344/1‑70 (MQ=255)                                                                          tCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGacg                                                                                                                    <  2:1238344/70‑1 (MQ=255)                                                                           ccACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGacggtaa                                                                                                                >  1:209934/1‑73 (MQ=255)                                                                           ccACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGacggtaa                                                                                                                <  2:209934/73‑1 (MQ=255)                                                                           ccACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCCTTTCTGACCGATTTTTTTGGTGATGATATTCCCTACACCGGCTAGCGCATCgg                                              >  1:1294758/1‑139 (MQ=255)                                                                           ccACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTg                                                                                >  1:1290392/1‑105 (MQ=255)                                                                           ccACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTg                                                                                <  2:1290392/105‑1 (MQ=255)                                                                           ccACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGtga                                                                            <  2:328866/109‑1 (MQ=255)                                                                           ccACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGtga                                                                            >  1:328866/1‑109 (MQ=255)                                                                                 gTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGt                                        >  2:1325322/1‑139 (MQ=255)                                                                                      gAACTGACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGc                                   >  2:821499/1‑139 (MQ=255)                                                                                           gACCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGcgacga                              >  2:1331736/1‑139 (MQ=255)                                                                                            aCCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCgacgac                             >  1:1330125/1‑139 (MQ=255)                                                                                            aCCGTTATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCgacgac                             >  1:485321/1‑139 (MQ=255)                                                                                                ttATAGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCg                         >  2:22057/1‑139 (MQ=255)                                                                                                 tataGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTa                                                                                                                          <  2:1546359/41‑1 (MQ=255)                                                                                                 tataGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTa                                                                                                                          >  1:1546359/1‑41 (MQ=255)                                                                                                 tataGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGa                                                                                            >  1:1157459/1‑71 (MQ=255)                                                                                                 tataGGTGTCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGa                                                                                            <  2:1157459/71‑1 (MQ=255)                                                                                                       tgtCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACg                  <  2:1358441/139‑1 (MQ=255)                                                                                                        gtCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGacg                                                                                                                    >  2:1114497/1‑40 (MQ=38)                                                                                                        gtCACCGCGATCGCGATGCTCCTGAATGGTGGTATCGacg                                                                                                                    <  1:1114497/40‑1 (MQ=38)                                                                                                           aCCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGa                                                                                                      >  2:1513327/1‑51 (MQ=255)                                                                                                           aCCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGa                                                                                                      <  1:1513327/51‑1 (MQ=255)                                                                                                           aCCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGtg                                                                             >  2:923387/1‑76 (MQ=255)                                                                                                           aCCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGtg                                                                             <  1:923387/76‑1 (MQ=255)                                                                                                           aCCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACAccgccg                                                        >  2:1563662/1‑97 (MQ=255)                                                                                                           aCCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACAccgccg                                                        <  1:1563662/97‑1 (MQ=255)                                                                                                           aCCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCgg                                              <  2:355746/107‑1 (MQ=255)                                                                                                           aCCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCgg                                              >  1:355746/1‑107 (MQ=255)                                                                                                           aCCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAAc                                           >  1:1466379/1‑110 (MQ=255)                                                                                                           aCCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAAc                                           <  2:1466379/110‑1 (MQ=255)                                                                                                           aCCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGAcc              >  1:316580/1‑139 (MQ=255)                                                                                                           aCCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGAcc              >  1:465093/1‑139 (MQ=255)                                                                                                           aCCGCGATCGCGATGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGAcc              >  1:1509568/1‑139 (MQ=255)                                                                                                                 atcgcgatGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGa                                                                                                      >  1:248604/1‑45 (MQ=255)                                                                                                                 atcgcgatGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGa                                                                                                      <  2:248604/45‑1 (MQ=255)                                                                                                                       atGCTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCAAGGACCACTTCAATACGt  >  1:535828/1‑139 (MQ=255)                                                                                                                          |                                                                                                                                        TTTGCGGGTCATCCTTTTCACGTTTTGCCAGGCTGCCGTAAGCTTTCATTCCCAGCACACCATCAATTAATGGTCCACTGGTAAAGAACTGACCGTTATAGGTGTCACCGCGATCGCGATGTTCCTGAATGGTGGTATCGACGGTAACGGTACCCGACCATTTCTGACCGATTTTTTTGGTGATGATATTCACTACACCGCCGAGCGCATCGGAACCGTACAGCGACGACATCGGGCCACGGACCACTTCAATACGT  >  NC_000913/2246121‑2246377

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