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%	101.1 / NA	30	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 (15/15);  total (15/15)

ACTCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGTTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATCGCG  >  NC_000913/1275682‑1275950                                                                                                                                      |                                                                                                                                        aCTCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTg                                                                                                                                    >  2:722320/1‑139 (MQ=255)   tCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGAt                                                                                                                                  <  2:341810/139‑1 (MQ=255)   tCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGAt                                                                                                                                  <  2:575402/139‑1 (MQ=255)      gCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTAc                                                                                                                               <  1:723767/139‑1 (MQ=255)      gCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTAc                                                                                                                               <  2:118842/139‑1 (MQ=255)              aaTACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGTCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGt                                                                                                                       <  2:651598/139‑1 (MQ=255)              aaTACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGt                                                                                                                       <  1:657312/139‑1 (MQ=255)                tACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTAt                                                                                                                     <  2:372547/139‑1 (MQ=255)                     tGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCt                                                                                                                >  1:683675/1‑139 (MQ=255)                     tGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCt                                                                                                                >  2:190306/1‑139 (MQ=255)                       ttCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTAt                                                                                                                     >  1:546471/1‑132 (MQ=255)                       ttCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCATCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTAt                                                                                                                     <  2:546471/132‑1 (MQ=255)                          gCCATTACTCGCTTCGCCAACCACGGTGATAACTGGTGCCATACTCATAAGCAGTTTTACGCCAGTTCGCAGCATCGGGTGCTCGTCCATCACCAGTATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGAc                                                                                                           >  2:712720/1‑139 (MQ=255)                                                    tGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATaa                                                                                 <  2:516798/139‑1 (MQ=255)                                                       tatCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAgg                                                                              <  2:683675/139‑1 (MQ=255)                                                              tGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGaccac                                                                       >  1:369358/1‑139 (MQ=255)                                                                   tACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTAtct                                                                                                                   <  2:28665/90‑1 (MQ=255)                                                                   tACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTAtct                                                                                                                   >  1:28665/1‑90 (MQ=255)                                                                   tACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACtt                                                                  >  1:72031/1‑139 (MQ=255)                                                                    aCTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTc                                                                 >  1:478497/1‑139 (MQ=255)                                                                       gATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGt                                                              >  1:487260/1‑139 (MQ=255)                                                                             cTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCAcc                                                        >  2:565840/1‑139 (MQ=255)                                                                                ttttACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGa                                                     >  2:683217/1‑139 (MQ=255)                                                                                     cGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCac                                                >  2:514333/1‑139 (MQ=255)                                                                                                                  aTCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATc                             >  1:663820/1‑129 (MQ=255)                                                                                                                  aTCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATc                             <  2:663820/129‑1 (MQ=255)                                                                                                                               tagCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGAt      >  1:287971/1‑139 (MQ=255)                                                                                                                                   cGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcgcg  <  1:271230/139‑1 (MQ=255)                                                                                                                                   cGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcgcg  <  2:478497/139‑1 (MQ=255)                                                                                                                                   cGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcgcg  <  1:514333/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                        ACTCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGTTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATCGCG  >  NC_000913/1275682‑1275950

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: ‑