Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,536,304	C→T	G105D (GGC→GAC)	narV ←	nitrate reductase Z subunit gamma

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,536,304	0	C	T	96.2%	84.0 / NA	26	G105D (GGC→GAC)	narV	nitrate reductase Z subunit gamma
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/11);  minor base G (0/1);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ACGCCGTTGAGCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGCCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTC  >  NC_000913/1536076‑1536530                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                   aCGCCGTTGAGCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGc                                                                                                                                                                                                               <  2:83574/250‑1 (MQ=255)                                ggAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAgttgt                                                                                                                                                                               >  1:234246/1‑251 (MQ=255)                                      tGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCt                                                                                                                                                                         >  1:287992/1‑251 (MQ=255)                                                   cgcCCAACCCACCAGCTTCATCATTGCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACGTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATGTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCaaaa                                                                                                                                                            <  1:295886/251‑1 (MQ=255)                                                         aCCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGTCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGGCAATCAACGTCAACACACCGCAAATGCCGCCGAGCACCAGCGCCATGAGTTGTTTCGCTTCCACTGGCAAATACCAcg                                                                                                                                                      <  1:266807/251‑1 (MQ=255)                                                            caccaGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGt                                                                                                                                                   <  2:230481/251‑1 (MQ=255)                                                             accaGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAACAACCACGCGt                                                                                                                                                   >  2:240817/1‑250 (MQ=255)                                                                               cTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTa                                                                                                                                 >  1:17120/1‑250 (MQ=255)                                                                                  ccATCAGGATACAGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGCCACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACAt                                                                                                                             >  2:306166/1‑251 (MQ=255)                                                                                   cATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATg                                                                                                                            <  2:278274/251‑1 (MQ=255)                                                                                        ggATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATTCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCTCGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACTCCGCAAATACCGCCTAGCACCATCGCCATGAGTTGTTTCTCTGCCACTGGCAAAATCCTCGCGGACATCCAGTGCGTCGTGAACATGCCGaa                                                                                                                       <  2:300422/251‑1 (MQ=255)                                                                                               gtgCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGt                                                                                                                  <  1:277002/249‑1 (MQ=255)                                                                                               gtgCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGt                                                                                                                  <  1:104403/249‑1 (MQ=255)                                                                                                     aaaaCGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCgaag                                                                                                          <  2:54502/251‑1 (MQ=255)                                                                                                        aCGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCgaagaa                                                                                                        >  1:27553/1‑250 (MQ=255)                                                                                                          ggTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAAt                                                                                                     >  1:228245/1‑251 (MQ=255)                                                                                                                                              aaaaTGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCatat                                                                 >  2:307064/1‑251 (MQ=255)                                                                                                                                                           tgatgatATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGt                                                    >  2:181330/1‑251 (MQ=255)                                                                                                                                                                  atCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCg                                              >  2:236044/1‑250 (MQ=255)                                                                                                                                                                                  gTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAATCCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTgagg                              >  2:211303/1‑250 (MQ=255)                                                                                                                                                                                            cgcgCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGcc                     >  2:111167/1‑249 (MQ=255)                                                                                                                                                                                             gcgcTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAgg                  <  1:273106/251‑1 (MQ=255)                                                                                                                                                                                               gcTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGgtg                <  1:240817/251‑1 (MQ=255)                                                                                                                                                                                               gcTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGgtg                >  2:46515/1‑251 (MQ=255)                                                                                                                                                                                                        tCAGCCTGCGCCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCg       >  1:88332/1‑251 (MQ=255)                                                                                                                                                                                                               gcgcCACAGCAGTCCTGCGCCGTCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTc  <  1:263593/249‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                   ACGCCGTTGAGCATCTCTGACGATCCACCACGGAAAGTGACAATGCTTTGCGCCCAACCCACCAGCTTCATCATTTCGCTACCATCAGGATACTGTGCTGAAAACGGTATGGTACTTAGCCCCAACAGACACTGGATCAGCAAAATGCTCATGATGATGATATCCGGCGTAGTGGAAGTGGCACGCACGCGCTGATTTGTCAGCCTGCGCCACAGCAGTCCTGCGCCGCCAATCAACGTCAACACACCGCAAATACCGCCGAGCACCATCGCCATGAGTTGTTTCGCTGCCACTGGCAAAAACCACGCGTACATCCAGTGCGGCGTTAACATGCCGAACAAGTGCCCGAAGAAAATCCCCAAAATGCCGATATGGAACAGATTCGACCATATCACCATCCCGCGTTTATCGAGCATTTGACTTGAGGAGGCGCGCCAGGTGTACTGCCCGTAGTC  >  NC_000913/1536076‑1536530

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