Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,275,815	T→C	E5E (GAA→GAG)	narL ←	response regulator in two‑component regulatory system with NarX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,275,815	0	T	C	100.0%	129.3 / NA	37	E5E (GAA→GAG)	narL	response regulator in two‑component regulatory system with NarX
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/20);  total (17/20)

GACTCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGTTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATCGCG  >  NC_000913/1275681‑1275950                                                                                                                                       |                                                                                                                                        gaCTCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCt                                                                                                                                     <  1:949969/139‑1 (MQ=255)      cGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTa                                                                                                                                <  2:369529/139‑1 (MQ=255)               aaTACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGt                                                                                                                       <  2:871495/139‑1 (MQ=255)                 tACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTAt                                                                                                                     <  2:955608/139‑1 (MQ=255)                 tACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTAt                                                                                                                     <  2:431231/139‑1 (MQ=255)                 tACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTAt                                                                                                                     >  1:431231/1‑139 (MQ=255)                      tGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTAt                                                                                                                     <  1:1048775/134‑1 (MQ=255)                      tGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTAt                                                                                                                     >  2:1048775/1‑134 (MQ=255)                      tGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTAt                                                                                                                     <  1:674908/134‑1 (MQ=255)                      tGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTAt                                                                                                                     >  2:674908/1‑134 (MQ=255)                             cATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGt                                                                                                         >  1:357337/1‑139 (MQ=255)                                             aaCCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGTAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTTTCTGAAAGTTTTTTc                                                                                         <  2:758627/139‑1 (MQ=255)                                             aaCCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTc                                                                                         <  2:463419/139‑1 (MQ=255)                                                 aCGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGAtt                                                                                                                                 <  1:453864/95‑1 (MQ=255)                                                 aCGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGAtt                                                                                                                                 >  2:453864/1‑95 (MQ=255)                                                                   aTACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTc                                                                                         >  1:254837/1‑117 (MQ=255)                                                                   aTACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTc                                                                                         <  2:254837/117‑1 (MQ=255)                                                                        gATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGt                                                              <  2:249644/139‑1 (MQ=255)                                                                         aTAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCa                                                                                                                            >  2:131029/1‑76 (MQ=255)                                                                         aTAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCa                                                                                                                            <  1:131029/76‑1 (MQ=255)                                                                           aaGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGcc                                                           >  1:244466/1‑139 (MQ=255)                                                                              cTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCAcc                                                        >  1:807284/1‑139 (MQ=255)                                                                              cTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCAcc                                                        >  1:266348/1‑139 (MQ=255)                                                                                 ttttACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCa                                                                                                                            <  2:834917/68‑1 (MQ=255)                                                                                 ttttACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCa                                                                                                                            >  1:834917/1‑68 (MQ=255)                                                                                           gTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGAcgg                                           >  2:299338/1‑139 (MQ=255)                                                                                                     tCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGc                                 >  2:509076/1‑139 (MQ=255)                                                                                                        ggTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAAt                              >  2:608716/1‑139 (MQ=255)                                                                                                                tCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGaccac                                                                       >  1:740518/1‑90 (MQ=255)                                                                                                                tCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGaccac                                                                       <  2:740518/90‑1 (MQ=255)                                                                                                                  aaTCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTaa                    >  1:1083782/1‑139 (MQ=255)                                                                                                                     cagcagGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACt                 <  2:357337/139‑1 (MQ=255)                                                                                                                            aTAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGTGCa          <  2:244466/139‑1 (MQ=255)                                                                                                                                  gCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcg    >  1:861658/1‑139 (MQ=255)                                                                                                                                   ccGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcgc   <  2:861658/139‑1 (MQ=255)                                                                                                                                   ccGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcgc   <  2:266348/139‑1 (MQ=255)                                                                                                                                    cGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcgcg  <  2:807284/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        GACTCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGTTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATCGCG  >  NC_000913/1275681‑1275950

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