Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,282,908	T→C	G1015G (GGT→GGC)	narG →	nitrate reductase A subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,282,908	0	T	C	97.7%	141.9 / ‑4.5	43	G1015G (GGT→GGC)	narG	nitrate reductase A subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (20/22);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAG  >  NC_000913/1282770‑1283036                                                                                                                                           |                                                                                                                                 gAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCtggtctggtc                                                                                                                                 >  1:235881/1‑140 (MQ=255)   aaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTggtccggtct                                                                                                                                >  2:656888/1‑139 (MQ=255)   aaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTggtccggtct                                                                                                                                <  1:656888/139‑1 (MQ=255)   aaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTggtccggtct                                                                                                                                >  2:334391/1‑139 (MQ=255)   aaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTggtccggtct                                                                                                                                <  1:334391/139‑1 (MQ=255)          ggCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATg                                                                                                                        >  1:90326/1‑140 (MQ=255)           gCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGa                                                                                                                       <  1:213583/140‑1 (MQ=255)               cccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAacac                                                                                                                   >  1:291199/1‑140 (MQ=255)                            tACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACg                                                                                                      <  2:244012/140‑1 (MQ=255)                                ggTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCgg                                                                                                  <  2:90326/140‑1 (MQ=255)                                    gTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTac                                                                                              >  1:54964/1‑140 (MQ=255)                                           tACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACg                                                                                       >  2:683173/1‑140 (MQ=255)                                            aCGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGt                                                                                      >  1:829703/1‑140 (MQ=255)                                            aCGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGt                                                                                      >  1:645400/1‑140 (MQ=255)                                               cATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCa                                                                                   >  2:630826/1‑140 (MQ=255)                                                                   agaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTa                                                               <  2:291199/140‑1 (MQ=255)                                                                     aCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGAACCCATGGCGTACg                                                             <  2:235881/140‑1 (MQ=255)                                                                     aCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGt                                                                <  2:35915/137‑1 (MQ=255)                                                                     aCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGt                                                                >  1:35915/1‑137 (MQ=255)                                                                     aCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACg                                                             >  1:604389/1‑140 (MQ=255)                                                                     aCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACg                                                             <  1:630826/140‑1 (MQ=255)                                                                            gaTCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                                                      >  1:285104/1‑140 (MQ=255)                                                                            gaTCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg                                                      >  2:404918/1‑140 (MQ=255)                                                                                    tCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAg                                              <  1:404918/140‑1 (MQ=255)                                                                                      cgcgATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCa                                            <  1:189708/140‑1 (MQ=255)                                                                                      cgcgATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCa                                            <  2:531545/140‑1 (MQ=255)                                                                                       gcgATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                           <  1:321947/140‑1 (MQ=255)                                                                                            atTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTa                                      <  1:133028/140‑1 (MQ=255)                                                                                              tCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATc                                    <  2:645400/140‑1 (MQ=255)                                                                                                 ggCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAgg                                 <  2:285104/140‑1 (MQ=255)                                                                                                  gCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                           <  1:572523/129‑1 (MQ=255)                                                                                                  gCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                           >  2:572523/1‑129 (MQ=255)                                                                                                  gCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGa                                <  2:54964/140‑1 (MQ=255)                                                                                                     caGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                           >  2:297430/1‑126 (MQ=255)                                                                                                     caGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                           <  1:297430/126‑1 (MQ=255)                                                                                                     caGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGa                                >  1:449948/1‑137 (MQ=255)                                                                                                     caGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGa                                <  2:449948/137‑1 (MQ=255)                                                                                                       gCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                           <  2:746490/124‑1 (MQ=255)                                                                                                       gCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                           >  1:746490/1‑124 (MQ=255)                                                                                                       gCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                           <  1:517228/124‑1 (MQ=255)                                                                                                       gCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                           >  2:517228/1‑124 (MQ=255)                                                                                                        ccGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                           >  2:97676/1‑123 (MQ=255)                                                                                                        ccGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCaa                                           <  1:97676/123‑1 (MQ=255)                                                                                                            gcAAAATTATCTCCAGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTgg                      <  2:653773/140‑1 (MQ=255)                                                                                                                          aGCCCGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCgg        >  2:809036/1‑140 (MQ=255)                                                                                                                            cccGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCg                 <  2:298262/129‑1 (MQ=255)                                                                                                                            cccGACCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCg                 >  1:298262/1‑129 (MQ=255)                                                                                                                                aCCTGGTCTGGCCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAg  >  1:632870/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                 GAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAG  >  NC_000913/1282770‑1283036

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