Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,585,544	T→C	intergenic (‑19/‑187)	aegA ← / → narQ	putative oxidoreductase Fe‑S subunit AegA/sensory histidine kinase NarQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,585,544	0	T	C	100.0%	193.5 / NA	54	intergenic (‑19/‑187)	aegA/narQ	putative oxidoreductase Fe‑S subunit AegA/sensory histidine kinase NarQ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/29);  total (25/29)

TGGTGTTGGCTCAGGACATGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCA  >  NC_000913/2585416‑2585670                                                                                                                                 |                                                                                                                               tGGTGTTGGCTCAGGACATGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCt                                                                                                                     <  2:679461/140‑1 (MQ=255)   gtgtTGGCTCAGGACATGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTa                                                                                                                   <  1:871818/140‑1 (MQ=255)   gtgtTGGCTCAGGACATGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTa                                                                                                                   <  2:914710/140‑1 (MQ=255)   gtgtTGGCTCAGGACATGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTa                                                                                                                   <  2:731241/140‑1 (MQ=255)   gtgtTGGCTCAGGACATGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTa                                                                                                                   <  1:59298/140‑1 (MQ=255)   gtgtTGGCTCAGGACATGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTa                                                                                                                   <  1:873605/140‑1 (MQ=255)       tGGCTCAGGACATGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCt                                                                                                               >  2:407335/1‑140 (MQ=255)             aGGACATGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTatat                                                                                                         <  1:406756/140‑1 (MQ=255)               gACATGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGa                                                                                                       <  1:383635/140‑1 (MQ=255)                   tGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTa                                                                                                                   >  2:77733/1‑124 (MQ=255)                   tGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTa                                                                                                                   <  1:77733/124‑1 (MQ=255)                   tGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAgtgt                                                                                                   >  2:1016221/1‑140 (MQ=255)                   tGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAgtgt                                                                                                   >  2:1063289/1‑140 (MQ=255)                   tGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAgtgt                                                                                                   >  2:158328/1‑140 (MQ=255)                         tcatcaTTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGttat                                                                                             >  2:1019146/1‑140 (MQ=255)                            tcatTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTAc                                                                                          <  1:407335/140‑1 (MQ=255)                            tcatTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTAc                                                                                          <  2:923093/140‑1 (MQ=255)                               ttGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCg                                                                                       <  1:1016221/140‑1 (MQ=255)                                       cATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAg                                                                               >  1:56906/1‑140 (MQ=255)                                       cATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAg                                                                               >  2:383598/1‑140 (MQ=255)                                       cATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAg                                                                               >  1:1050152/1‑140 (MQ=255)                                                ggCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCgg                                                                      >  2:997201/1‑140 (MQ=255)                                                ggCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCgg                                                                      <  1:742891/140‑1 (MQ=255)                                                   gATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCaa                                                                   >  2:976806/1‑140 (MQ=255)                                                    aTTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGa                                                                           >  1:864891/1‑131 (MQ=255)                                                    aTTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGa                                                                           <  2:864891/131‑1 (MQ=255)                                                          caAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCAt                                                            <  2:56906/140‑1 (MQ=255)                                                              cATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGc                                                        <  2:1050152/140‑1 (MQ=255)                                                               aTGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCa                                                                    <  2:778210/127‑1 (MQ=255)                                                               aTGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCa                                                                    >  1:778210/1‑127 (MQ=255)                                                               aTGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGcc                                                       <  1:976806/140‑1 (MQ=255)                                                                           aCATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGcc                                                       >  2:890686/1‑128 (MQ=255)                                                                           aCATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGcc                                                       <  1:890686/128‑1 (MQ=255)                                                                             aTTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACt                                         >  2:1039075/1‑140 (MQ=255)                                                                              ttGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTa                                        >  2:575113/1‑140 (MQ=255)                                                                                gCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTAtt                                      >  2:413455/1‑140 (MQ=255)                                                                                   gACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCg                                   <  1:383598/140‑1 (MQ=255)                                                                                     cTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCt                                 <  1:413455/140‑1 (MQ=255)                                                                                         tGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATaa                             >  1:40201/1‑140 (MQ=255)                                                                                           gCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAAc                           <  2:1033657/140‑1 (MQ=255)                                                                                           gCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAAc                           >  1:761154/1‑140 (MQ=255)                                                                                                    aaaCGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGa                  <  2:40201/140‑1 (MQ=255)                                                                                                    aaaCGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGa                  <  1:71245/140‑1 (MQ=255)                                                                                                     aaCGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGAt                 <  1:575113/140‑1 (MQ=255)                                                                                                     aaCGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGAt                 >  2:1066609/1‑140 (MQ=255)                                                                                                      aCGATTCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATc                <  1:755246/140‑1 (MQ=255)                                                                                                          ttCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAaat            >  1:210846/1‑140 (MQ=255)                                                                                                          ttCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAaat            <  1:112659/140‑1 (MQ=255)                                                                                                          ttCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAaat            >  2:418221/1‑140 (MQ=255)                                                                                                          ttCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAaat            <  1:569262/140‑1 (MQ=255)                                                                                                          ttCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAaat            >  2:318911/1‑140 (MQ=255)                                                                                                          ttCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAaat            <  1:967931/140‑1 (MQ=255)                                                                                                          ttCATAACGACCTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAaat            >  1:778553/1‑140 (MQ=255)                                                                                                                    ccTTCTTTTTTGGCTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCa  >  1:259216/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                               TGGTGTTGGCTCAGGACATGTTGCTCATCATTGTGAGCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCA  >  NC_000913/2585416‑2585670

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