Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,278,990	G→A	G345D (GGC→GAC)	narK →	nitrate/nitrite transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,278,990	0	G	A	98.0%	172.4 / NA	51	G345D (GGC→GAC)	narK	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (21/29);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAA  >  NC_000913/1278856‑1279120                                                                                                                                       |                                                                                                                                   tGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGa                                                                                                                                <  1:61597/139‑1 (MQ=255)      cgcTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCtt                                                                                                                           >  1:120940/1‑139 (MQ=255)       gcTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTc                                                                                                                          <  1:1088206/139‑1 (MQ=255)        cTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCa                                                                                                                         >  2:331526/1‑139 (MQ=255)          ggCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATg                                                                                                                       <  2:353251/139‑1 (MQ=255)                   tGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGttcttc                                                                                                              <  1:876423/139‑1 (MQ=255)                             gCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTcc                                                                                                    <  1:232693/139‑1 (MQ=255)                                  aTCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGgcg                                                                                               <  1:331526/139‑1 (MQ=255)                                           tCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTg                                                                                      >  1:116818/1‑139 (MQ=255)                                           tCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTg                                                                                      >  1:608086/1‑139 (MQ=255)                                               ggCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAg                                                                                  >  2:492683/1‑139 (MQ=255)                                                    aaCTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAgctggg                                                                             <  1:913467/139‑1 (MQ=255)                                                    aaCTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAgctggg                                                                             <  2:120940/139‑1 (MQ=255)                                                    aaCTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAgctggg                                                                             <  2:1033537/139‑1 (MQ=255)                                                          tgtCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGt                                                                       >  1:1176072/1‑139 (MQ=255)                                                            tCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAg                                                                     <  1:862868/139‑1 (MQ=255)                                                                    gTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCa                                                             <  1:288395/139‑1 (MQ=255)                                                                    gTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCa                                                             >  1:83405/1‑139 (MQ=255)                                                                            tATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACgg                                                                                                                                 >  1:887954/1‑63 (MQ=255)                                                                            tATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACgg                                                                                                                                 <  2:887954/63‑1 (MQ=255)                                                                             aTTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGgcg                                                                                               <  2:1170338/96‑1 (MQ=255)                                                                             aTTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGgcg                                                                                               >  1:1170338/1‑96 (MQ=255)                                                                              ttCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAgat                                                   <  2:116818/139‑1 (MQ=255)                                                                              ttCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAgat                                                   >  2:1168099/1‑139 (MQ=255)                                                                              ttCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAgat                                                   >  1:668154/1‑139 (MQ=255)                                                                                   aTGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTc                                                               <  2:874245/122‑1 (MQ=255)                                                                                   aTGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTc                                                               >  1:874245/1‑122 (MQ=255)                                                                                   aTGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGAttt                                              >  2:160417/1‑139 (MQ=255)                                                                                   aTGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGAttt                                              >  1:370918/1‑139 (MQ=255)                                                                                   aTGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGAttt                                              <  2:608086/139‑1 (MQ=255)                                                                                           tttCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATc                                      >  1:40986/1‑139 (MQ=255)                                                                                                     ctgctgTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTTATCTTCCGTAAAc                            <  2:208439/139‑1 (MQ=255)                                                                                                     ctgctgTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAAc                            <  1:854791/139‑1 (MQ=255)                                                                                                     ctgctgTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAAc                            <  1:76027/139‑1 (MQ=255)                                                                                                     ctgctgTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAAc                            <  2:25789/139‑1 (MQ=255)                                                                                                     ctgctgTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAAc                            <  2:370918/139‑1 (MQ=255)                                                                                                     ctgctgTTCCTGACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAAc                            <  2:1176072/139‑1 (MQ=255)                                                                                                                gACCTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGAt                 >  2:978236/1‑139 (MQ=255)                                                                                                                  ccTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCa                                                             <  2:1164594/93‑1 (MQ=255)                                                                                                                  ccTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCa                                                             >  1:1164594/1‑93 (MQ=255)                                                                                                                  ccTTACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATcg               >  1:791778/1‑139 (MQ=255)                                                                                                                     tACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGAt                 >  1:818364/1‑134 (MQ=255)                                                                                                                     tACCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGt            <  1:492683/139‑1 (MQ=255)                                                                                                                     tACCGACTAACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGCGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGAt                 <  2:818364/134‑1 (MQ=255)                                                                                                                      aCCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTgg                                                                  <  1:692324/84‑1 (MQ=255)                                                                                                                      aCCGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTgg                                                                  >  2:692324/1‑84 (MQ=255)                                                                                                                       ccGACTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGa          <  2:83405/139‑1 (MQ=255)                                                                                                                          aCTGACGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAg       <  1:1021504/139‑1 (MQ=255)                                                                                                                               cGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAg                                                                                  <  2:1044839/59‑1 (MQ=255)                                                                                                                               cGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAg                                                                                  >  1:1044839/1‑59 (MQ=255)                                                                                                                               cGGGCAGGACGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGaa  <  1:1168099/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                   TGGTGCGCTGGCGCGTTCTGCAGGTGGTGCATTATCTGACCGTCTGGGCGGAACTCGTGTCACGCTGGTGAACTTTATTCTGATGGCGATTTTCAGCGGCCTGCTGTTCCTGACCTTACCGACTGACGGGCAGGGCGGAAGCTTCATGGCGTTCTTCGCGGTCTTCCTGGCGCTGTTCCTGACAGCTGGGCTGGGTAGTGGTTCCACTTTCCAGATGATTTCAGTGATCTTCCGTAAACTGACAATGGATCGCGTGAAAGCAGAA  >  NC_000913/1278856‑1279120

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