Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	485,798	G→A	R13H (CGC→CAC)	acrR →	DNA‑binding transcriptional repressor AcrR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	485,798	0	G	A	100.0%	178.9 / NA	52	R13H (CGC→CAC)	acrR	DNA‑binding transcriptional repressor AcrR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/39);  total (13/39)

CTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGTATCATCCACCTCG  >  NC_000913/485736‑485859                                                               |                                                              cTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCAcc                                                             <  1:1828155/65‑1 (MQ=255)    aCGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAac                                                          <  1:65757/65‑1 (MQ=255)     cGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACTAaca                                                         >  1:579040/1‑65 (MQ=255)     cGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAaca                                                         >  1:2636869/1‑65 (MQ=255)     cGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAaca                                                         <  1:881916/65‑1 (MQ=255)     cGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAaca                                                         >  1:658880/1‑65 (MQ=255)     cGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAaca                                                         >  1:2707526/1‑65 (MQ=255)     cGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAaca                                                         <  1:356830/65‑1 (MQ=255)      gCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAacac                                                        >  1:3230764/1‑65 (MQ=255)       ccTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAAcaca                                                       <  1:432104/65‑1 (MQ=255)       ccTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAAcaca                                                       >  1:468535/1‑65 (MQ=255)       ccTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAAcaca                                                       <  1:2171004/65‑1 (MQ=255)         tGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATc                                                     <  1:3175107/65‑1 (MQ=255)         tGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATc                                                     <  1:2980880/65‑1 (MQ=255)         tGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATc                                                     <  1:1442965/65‑1 (MQ=255)         tGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATc                                                     <  1:2704593/65‑1 (MQ=255)           tAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCt                                                   <  1:2966000/65‑1 (MQ=255)              aTTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGa                                                >  1:1708741/1‑65 (MQ=255)              aTTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGa                                                >  1:2046216/1‑65 (MQ=255)               ttCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGAt                                               <  1:1826467/65‑1 (MQ=255)                   cGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTgg                                           >  1:463381/1‑65 (MQ=255)                    gAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGc                                          <  1:2881422/65‑1 (MQ=255)                    gAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGc                                          <  1:2964558/65‑1 (MQ=255)                    gAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGc                                          <  1:2723397/65‑1 (MQ=255)                    gAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGc                                          <  1:99768/65‑1 (MQ=255)                    gAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGc                                          <  1:1285529/65‑1 (MQ=255)                     aaCATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGct                                         <  1:342297/65‑1 (MQ=255)                     aaCATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGct                                         <  1:1832861/65‑1 (MQ=255)                     aaCATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGct                                         <  1:1596924/65‑1 (MQ=255)                     aaCATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGct                                         <  1:756919/65‑1 (MQ=255)                     aaCATACGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGct                                         <  1:1488379/65‑1 (MQ=255)                      aCATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGctc                                        >  1:124977/1‑65 (MQ=255)                        atatGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTa                                      >  1:3246992/1‑65 (MQ=255)                          atGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACg                                    <  1:609964/65‑1 (MQ=255)                            ggCACGAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTc                                  >  1:2971241/1‑65 (MQ=255)                                 gAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTtc                             <  1:373242/65‑1 (MQ=255)                                 gAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTtc                             <  1:1816838/65‑1 (MQ=255)                                 gAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTtc                             <  1:844804/65‑1 (MQ=255)                                 gAAAAACCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTtc                             <  1:848656/65‑1 (MQ=255)                                      aCCAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTcaca                        <  1:1912069/65‑1 (MQ=255)                                       ccAAACAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCAcag                       >  1:2151816/1‑65 (MQ=255)                                          aaCAAGAAGCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCAcagcag                    <  1:1363590/65‑1 (MQ=255)                                             aagaagCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGGAgggg                 <  1:656915/65‑1 (MQ=255)                                             aagaagCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAgggg                 <  1:3140957/65‑1 (MQ=255)                                             aagaagCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAgggg                 <  1:2753919/65‑1 (MQ=255)                                             aagaagCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAgggg                 <  1:2056348/65‑1 (MQ=255)                                              agaagCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAggggg                <  1:463636/65‑1 (MQ=255)                                               gaagCGCAAGAAACGCACCAACGCATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGt               <  1:1203069/65‑1 (MQ=255)                                               gaagCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGt               <  1:3091372/65‑1 (MQ=255)                                               gaagCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGt               <  1:1906461/65‑1 (MQ=255)                                               gaagCGCAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGt               <  1:626576/65‑1 (MQ=255)                                                  gcgcAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGTatc            <  1:1600864/65‑1 (MQ=255)                                                  gcgcAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGTatc            <  1:1501535/65‑1 (MQ=255)                                                  gcgcAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGTatc            <  1:2922341/65‑1 (MQ=255)                                                   cgcAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGTatca           <  1:1056671/65‑1 (MQ=255)                                                   cgcAAGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGTatca           <  1:2142364/65‑1 (MQ=255)                                                       aGAAACGCACCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGTatcatcta       >  1:1957543/1‑63 (MQ=255)                                                            cgcaCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGTATCACCCACCTCg  <  1:1569889/65‑1 (MQ=255)                                                               |                                                              CTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCAGGGGGTATCATCCACCTCG  >  NC_000913/485736‑485859

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