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%	122.1 / 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 (16/21);  total (16/21)

ACTCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGTTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATCGCG  >  NC_000913/1275682‑1275950                                                                                                                                      |                                                                                                                                        aCTCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTg                                                                                                                                    >  1:526659/1‑139 (MQ=255)  cTCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGa                                                                                                                                   <  1:239550/139‑1 (MQ=255)   tCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGAt                                                                                                                                  >  1:146850/1‑139 (MQ=255)       ccAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACt                                                                                                                              <  1:153619/139‑1 (MQ=255)       ccAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACt                                                                                                                              <  2:534782/139‑1 (MQ=255)         aGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCa                                                                                                                            <  2:526659/139‑1 (MQ=255)          gTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCAt                                                                                                                           <  2:331237/139‑1 (MQ=255)            tCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATgg                                                                                                                         <  1:1075350/139‑1 (MQ=255)              aaTACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGt                                                                                                                       <  2:200215/139‑1 (MQ=255)              aaTACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGt                                                                                                                       <  1:997912/139‑1 (MQ=255)              aaTACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGt                                                                                                                       <  1:866634/139‑1 (MQ=255)                     tGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCt                                                                                                                >  2:64904/1‑139 (MQ=255)                       ttCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATGCTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCa                                                                                                                            >  2:641808/1‑125 (MQ=255)                       ttCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCa                                                                                                                            <  1:641808/125‑1 (MQ=255)                                         gCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGtttt                                                                                            >  2:969283/1‑139 (MQ=255)                                            aaCCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTc                                                                                         <  1:969283/139‑1 (MQ=255)                                                         tCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTc                                                                                         >  2:559166/1‑126 (MQ=255)                                                         tCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTc                                                                                         <  1:559166/126‑1 (MQ=255)                                                                   tACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGt                                                                                                                       >  1:274017/1‑86 (MQ=255)                                                                   tACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGt                                                                                                                       <  2:274017/86‑1 (MQ=255)                                                                   tACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTc                                                                                         >  1:690528/1‑116 (MQ=255)                                                                   tACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTc                                                                                         <  2:690528/116‑1 (MQ=255)                                                                        aTAAGCTTTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAGTCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTg                                                             >  2:517525/1‑139 (MQ=255)                                                                         tAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGc                                                            >  1:1045877/1‑139 (MQ=255)                                                                             cTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCCCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCAcc                                                        >  2:672828/1‑139 (MQ=255)                                                                                          gTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTAt                                                                                                                     <  2:455228/65‑1 (MQ=255)                                                                                          gTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTAt                                                                                                                     >  1:455228/1‑65 (MQ=255)                                                                                                    tCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCa                                                                                                                            <  1:498021/48‑1 (MQ=255)                                                                                                    tCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCa                                                                                                                            >  2:498021/1‑48 (MQ=255)                                                                                                            tCGTCAATCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCt                         >  2:489107/1‑139 (MQ=255)                                                                                                                   tCAGCAGGATAGTAGCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAAc                  >  1:477377/1‑139 (MQ=255)                                                                                                                                 gCCGGCTCCTGATTAGTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGGCCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcg    <  1:517525/139‑1 (MQ=255)                                                                                                                                 gCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcg    <  1:690153/139‑1 (MQ=255)                                                                                                                                 gCCGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcg    >  1:470529/1‑139 (MQ=255)                                                                                                                                   cGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGCCGGCGGACGCGCCAAGCGCCTCGTAAACTTTGCGCACGATcgcg  <  2:1093199/139‑1 (MQ=255)                                                                                                                                   cGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcgcg  <  2:477377/139‑1 (MQ=255)                                                                                                                                   cGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcgcg  <  2:470529/139‑1 (MQ=255)                                                                                                                                   cGGCTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATcgcg  <  2:146850/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                        ACTCCGCCAGTTCAATACCCTGTTCGCCATTACTCGCTTCGCCAACCACGGTGATATCTGGTGCCATACTGATAAGCTGTTTTACGCCAGTTCGCAGCATCGGGTGATCGTCAATCAGCAGGATAGTAGCCGGTTCCTGATTACTCATGGGTATCTCCTTGGACGTCTGTGAAAGTTTTTTCGGGAATAAAGGTGACCACCACTTCGGTGCCACCTGATTCACGACGGCGGACGCGGCAATCGCCTCGTAAACTTTGCGCACGATCGCG  >  NC_000913/1275682‑1275950

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