Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	991,358	A→G	D246G (GAC→GGC)	pepN →	aminopeptidase N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	991,358	0	A	G	95.8%	74.3 / NA	24	D246G (GAC→GGC)	pepN	aminopeptidase N
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (6/17);  minor base C (0/1);  total (6/18)
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

CGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGACCTCGACATCTATATGATCGTCGCGGTGGATTTCTTCAATATGGGCGCAATGGAGAATAAGGGGCTGAATATCTTTAACTCCAAATATGTGCTGGCCCGCACCGA  >  NC_000913/991229‑991463                                                                                                                                  |                                                                                                          cgcgTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCt                                                                                                        <  1:322744/133‑1 (MQ=255)  gcgTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGAGGAAGAACCCTTAGGTCTGGAGGATGGCCTc                                                                                                       >  1:129970/1‑133 (MQ=255)   cgTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGATTTTGCGTGGAGTTTGGCCTCg                                                                                                      >  1:167110/1‑133 (MQ=255)    gTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGa                                                                                                     <  2:52485/133‑1 (MQ=255)     ttCTGGTCGCTAAGTCGCACTGGACCTGAACGCCGATCGCGCCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGAc                                                                                                    <  2:146133/133‑1 (MQ=255)     ttCTGGTCGCTAAGTAGCACTGGAACCGCACGACGATCGCGCCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGAc                                                                                                    <  2:238153/133‑1 (MQ=255)     ttCTGGTCGCGAAGTAGCGCTGGAGCTGTAGGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGAAGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGAc                                                                                                    <  1:233485/133‑1 (MQ=255)     ttCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAACTGGGATGAAGAACGCTTCGGCCTGGAGTATGGCCTCGAc                                                                                                    <  2:338464/133‑1 (MQ=255)     ttCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGAc                                                                                                    <  1:125077/133‑1 (MQ=255)     ttCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGAc                                                                                                    <  1:293941/133‑1 (MQ=255)                           gAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGAGATCGATGTGATCGACGCGGTgg                                                                              >  1:130079/1‑133 (MQ=255)                           gAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGACATCTATATGATCGTCGCGGTgg                                                                              >  2:21543/1‑133 (MQ=255)                                        aTCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGACATCTATATGATCGTCGCGGTGGATTTCTTCAatat                                                                 >  1:331097/1‑133 (MQ=255)                                          cgcgGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGACATCTATATGATCGTCGCGGTGGATTTCTTCAATATgg                                                               <  1:21543/133‑1 (MQ=255)                                             ggCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGACATCTATATGATCGTCGCGGTGGATTTCTTCAATATGGgcg                                                            <  1:443750/133‑1 (MQ=255)                                               cAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGACATCTATATGATCGTCGCGGTGGATTTCTTCAATATGGGCGCa                                                          <  2:8356/133‑1 (MQ=255)                                                  ccTTGATCGCCCGCCCTGGGCGATCACCTCGCCGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGACATCTATATGATCGTCGCGGTGGATTTCTTCAATATGGGCGCAATg                                                       <  2:130079/133‑1 (MQ=255)                                                   cTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGACATCTATATGATCGTCGCGGTGGATTTCTTCAATATGGGCGCAATgg                                                      <  2:31265/133‑1 (MQ=255)                                                    ttGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGACATCTATATGATCGTCGCGGTGGATTTCTTCAATATGGGCGCAATGga                                                     <  1:127126/133‑1 (MQ=255)                                                                          gACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGCCCTCGACATCTATATGATCGTCGCGGTGGATTTCTTCAATATGGGCGCAATGGAGAATAAGGGGCTGAATATCttt                               <  2:323958/133‑1 (MQ=255)                                                                               cGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGGCCTCGACATCTATATGATCGTCGCGGTGGATTTCTTCAATATGGGCGCAATGGAGAATAAGGGGCTGAATATCTTTAACTc                          <  1:378146/133‑1 (MQ=255)                                                                                 cTGAAAAACTCCATGCAATGGGATCAAGCACGCTTTGGCCTCGAGTATGGCCTCGACATCTATATGATCGTCGCGGTGGATTTCTTCAATATGGGCGCAATGGAGAATAAGGGGCTGAATATCTTTAACTCCa                        <  2:129970/133‑1 (MQ=255)                                                                                    aaaaaCTCCAGGAAATGGGATGAAGAACTCTTTGGCCTGGAATATAGCCTCGACATCTATATGATCGTCGCGGTGGATTTCTTCAATATGGGCGCAATGGAGAATAAGGGGCTGAATATCTTTAACTCCAAat                     <  2:167110/133‑1 (MQ=255)                                                                                                       aTGAAGAACGCTTTGGCCTGGAGTATGGCCTCGACATCTATATGATCGTCGCGGTGGATTTCTTCAATATGGGCGCAATGGAGAATAAGGGGCTGAATATCTTTAACTCCAAATATGTGCTGGCCCGCACCGa  >  2:406580/1‑133 (MQ=255)                                                                                                                                  |                                                                                                          CGCGTTCTGGTCGCGAAGTAGCACTGGAGCTGTACGTCGATCGCGGCAACCTTGATCGCGCGCCGTGGGCGATGACCTCGCTGAAAAACTCCATGAAATGGGATGAAGAACGCTTTGGCCTGGAGTATGACCTCGACATCTATATGATCGTCGCGGTGGATTTCTTCAATATGGGCGCAATGGAGAATAAGGGGCTGAATATCTTTAACTCCAAATATGTGCTGGCCCGCACCGA  >  NC_000913/991229‑991463

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