Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,246,731	A→G	C13C (TGT→TGC)	cirA ←	ferric dihyroxybenzoylserine outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,246,731	0	A	G	100.0%	120.2 / NA	36	C13C (TGT→TGC)	cirA	ferric dihyroxybenzoylserine outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/20);  total (16/20)

AGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGG  >  NC_000913/2246594‑2246856                                                                                                                                          |                                                                                                                              aGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAgca                                                                                                                             <  1:335448/140‑1 (MQ=255)  ggTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAgcac                                                                                                                            <  1:58444/140‑1 (MQ=255)                 tGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACg                                                                                                             >  1:36161/1‑140 (MQ=255)                    cGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGaa                                                                                                           >  2:481182/1‑139 (MQ=255)                    cGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGaa                                                                                                           <  1:481182/139‑1 (MQ=255)                      cTGATACTGGCAGGCGCATTTTTAAGATTTTGTTCCACCGAAGATGCAGTGACAACCATCGTTTTGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAgg                                                                                                        >  1:764365/1‑140 (MQ=255)                        gATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGt                                                                                                      <  1:99767/140‑1 (MQ=255)                         aTACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGtt                                                                                                     <  2:165409/140‑1 (MQ=255)                               gCAGGTGCCTCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGaaa                                                                                                          >  1:207592/1‑129 (MQ=255)                               gCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGCCAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCt                                                                                               <  2:764365/140‑1 (MQ=255)                               gCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGaaa                                                                                                          <  2:207592/129‑1 (MQ=255)                                   gTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGTCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAAc                                                                                           <  2:514798/140‑1 (MQ=255)                                    tGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACa                                                                                          <  2:36161/140‑1 (MQ=255)                                        tCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTcc                                                                                      <  2:305211/140‑1 (MQ=255)                                        tCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTcc                                                                                      <  1:705498/140‑1 (MQ=255)                                             aaGATTTTGTTCCACGGAAGATGCAATGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTaaa                                                                                            <  2:146621/129‑1 (MQ=255)                                             aaGATTTTGTTCCACGGAAGATGCAATGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTaaa                                                                                            >  1:146621/1‑129 (MQ=255)                                                   ttGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATg                                                                           >  1:670795/1‑140 (MQ=255)                                                   ttGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATg                                                                           <  2:656391/140‑1 (MQ=255)                                                    tGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGa                                                                          >  2:189085/1‑140 (MQ=255)                                                           cGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACt                                                                   >  2:38679/1‑140 (MQ=255)                                                             gAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTAc                                                                 <  2:670795/140‑1 (MQ=255)                                                                 aTGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaa                                                             <  1:38679/140‑1 (MQ=255)                                                                       tGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaataaaa                                                       >  2:196790/1‑140 (MQ=255)                                                                            aCCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGt                                                  >  1:136813/1‑140 (MQ=255)                                                                            aCCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGt                                                  >  1:492567/1‑140 (MQ=255)                                                                                  gTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGc                                            <  2:168805/140‑1 (MQ=255)                                                                                       gCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTcaca                                       <  1:196790/140‑1 (MQ=255)                                                                                          atcatcatcGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACAtct                                    <  2:492567/140‑1 (MQ=255)                                                                                          atcatcatcGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACAtct                                    <  2:136813/140‑1 (MQ=255)                                                                                                    aCCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTc                          <  1:406871/140‑1 (MQ=255)                                                                                                       gCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTcccc                       >  1:600700/1‑140 (MQ=255)                                                                                                       gCTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTcccc                       >  2:97378/1‑140 (MQ=255)                                                                                                        cTAACACAGGCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCt                      >  1:371116/1‑140 (MQ=255)                                                                                                                 gCCATGCACAAGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACgg             <  2:242898/140‑1 (MQ=255)                                                                                                                           aGAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATcgcg   >  1:518866/1‑140 (MQ=255)                                                                                                                            gAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCgg  <  1:83920/140‑1 (MQ=255)                                                                                                                            gAAATAGCGGACAAGCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCgg  >  2:83920/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                              AGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGG  >  NC_000913/2246594‑2246856

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