Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,247,638	G→T	10.0%	E223D (GAG→GAT)	argG →	argininosuccinate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,247,638	0	G	T	10.0%	113.8 / 4.2	40	E223D (GAG→GAT)	argG	argininosuccinate synthase
Reads supporting (aligned to +/- strand):  ref base G (19/17);  new base T (2/2);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

GAAGTCGCCGTACTCGATGGACGCCAACCTGCTGCACATCTCCTACGAGGGCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTATCGGTGGTGCCAACGGCATCGGCCGTCTGGACATCGTCGAAAACCGTTACGTCGGCATG  >  NC_002947/1247491‑1247779                                                                                                                                                    |                                                                                                                                              gaagTCGCCGTACTCGATGGACGCCAACCTGCTGCACATCTCCTACGAGGGCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGc                                                                                                                                              >  2:129419/1‑149 (MQ=255)       gCCGTACTCGATGGACGCCAACCTGCTGCACATCTCCTACGAGGGCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGAcc                                                                                                                                         >  3:279593/1‑148 (MQ=255)                gATGGACGCCAACCTGCTGCACATCTCCTACGAGGGCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAAt                                                                                                                                <  4:279593/148‑1 (MQ=255)                gATGGACGCCAACCTGCTGCACATCTCCTACGAGGGCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATg                                                                                                                               >  1:241510/1‑149 (MQ=255)                           aCCTGCTGCACATCTCCTACGAGGGCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGtt                                                                                                                    <  4:321768/149‑1 (MQ=255)                                ctgcACATCTCCTACGAGGGCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCAt                                                                                                               <  3:366742/149‑1 (MQ=255)                                    acaTCTCCTACGAGGGCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGAc                                                                                                           >  5:125058/1‑149 (MQ=255)                                      aTCTCCTACGAGGGCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGa                                                                                                                                           <  2:312774/115‑1 (MQ=255)                                           cTACGAGGGCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTc                                                                                                     <  4:323653/148‑1 (MQ=255)                                            tACGAGGGCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTc                                                                                                     >  6:52543/1‑147 (MQ=255)                                              cGAGGGCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAg                                                                                                                                               <  8:66516/103‑1 (MQ=255)                                              cGAGGGCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAg                                                                                                                                               >  7:66516/1‑103 (MQ=255)                                              cGAGGGCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCgaga                                                                                                 <  3:183122/149‑1 (MQ=255)                                                  ggCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAAt                                                                                              <  1:129419/148‑1 (MQ=255)                                                     ggTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATcccc                                                                                          >  7:228082/1‑149 (MQ=255)                                                                  aTACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTg                                                                             <  6:125058/149‑1 (MQ=255)                                                                   tACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTgg                                                                            <  2:17047/149‑1 (MQ=255)                                                                      cTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGt                                                                                                                             <  7:209502/97‑1 (MQ=255)                                                                      cTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGt                                                                                                                             >  8:209502/1‑97 (MQ=255)                                                                        ggACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGATAATGCCACGGAACAGGCTACCTACATCCAGCTGACATACCGAAATGTTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCCCCGTCCTGGAAGa                                                                        >  7:273364/1‑148 (MQ=255)                                                                        ggACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGAc                                                                       >  6:100731/1‑149 (MQ=255)                                                                        ggACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGAc                                                                       >  1:181156/1‑149 (MQ=255)                                                                            cGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTg                                                                    >  3:355337/1‑148 (MQ=255)                                                                               gCACGAAGAAGACATCTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAAcc                                                                <  8:228082/149‑1 (MQ=255)                                                                                   gaagaaGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGATCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTAt                                                            >  5:52287/1‑149 (MQ=255)                                                                                   gaagaaGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGATCTGACCTACCGCAAAGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGCCAGACCTGAACCGTAt                                                            >  5:46707/1‑149 (MQ=255)                                                                                    aagaagACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTAt                                                            >  1:228394/1‑148 (MQ=255)                                                                                    aagaagACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTAt                                                            <  2:228394/148‑1 (MQ=255)                                                                                       aagaCATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGAc                                                                       <  4:139114/134‑1 (MQ=255)                                                                                       aagaCATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATACCCGGCCACCGTCCTGGCAGAc                                                                       >  3:139114/1‑134 (MQ=255)                                                                                                     tGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTATCGGTGGTGCCAACGGCAt                                          >  8:7681/1‑149 (MQ=255)                                                                                                                  cTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTATCGGTGGTGCCAACGGCATCGGCCGTCTGGa                              >  7:37815/1‑148 (MQ=255)                                                                                                                     agaATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTATCGGTGGTGCCAACGGCATCGGCCGTCTGGACAtc                          >  2:192316/1‑149 (MQ=255)                                                                                                                           ccccGGACCAGGCTACCTACATCGATCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTATCGGTGGTGCCAACGGCATCGGGCGTCTGGACATCGTCGaa                    <  6:46707/149‑1 (MQ=255)                                                                                                                           ccccGGACCAGGCTACCTACATCGATCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTATCGGTGGTGCCAACGGCATCGGCCGTCTGGACATCGTCGaa                    <  6:52287/149‑1 (MQ=255)                                                                                                                               ggACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTATCGGTGGTGCCAACGGCATCGGCCGTCTGGACATCGTCGAAAAcc                <  1:192316/149‑1 (MQ=255)                                                                                                                                  ccAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTATCGGTGGTGCCAACGGCATCGGCCGTCTGGACATCGTCGAAAACCGtt             <  2:181156/149‑1 (MQ=255)                                                                                                                                  ccAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTATCGGTGGTGCCAACGGCATCGGCCGTCTGGACATCGTCGAAAACCGtt             <  5:52543/149‑1 (MQ=255)                                                                                                                                     ggCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTATCGGTGGTGCCAACGGCATCGGCCGTCTGGACATCGTCGAAAACCGTTa            >  3:196421/1‑147 (MQ=255)                                                                                                                                     ggCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTATCGGTGGTGCCAACGGCATCGGCCGTCTGGACATCGTCGAAAACCGTTa            <  4:196421/147‑1 (MQ=255)                                                                                                                                              cATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTATCGGTGGTGCCAACGGCATCGGCCGTCTGGACATCGTCGAAAACCGTTACGTCGGCATg  >  6:161327/1‑148 (MQ=255)                                                                                                                                               aTCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCATGGCAGACCTGAAAAGTATCGGTGGTGCCAACGGCATCGGCCGTCTGGACATCGTCGAAAACCGTTACGTCGGCATg  >  3:365570/1‑147 (MQ=255)                                                                                                                                                    |                                                                                                                                              GAAGTCGCCGTACTCGATGGACGCCAACCTGCTGCACATCTCCTACGAGGGCGGTGTCCTGGAAGATACCTGGACCGAGCACGAAGAAGACATGTGGCGCTGGAGTGTCTCGCCTGAGAATGCCCCGGACCAGGCTACCTACATCGAGCTGACCTACCGCAATGGTGACATCGTTGCCATCGACGGCGTCGAGAAATCCCCGGCCACCGTCCTGGCAGACCTGAACCGTATCGGTGGTGCCAACGGCATCGGCCGTCTGGACATCGTCGAAAACCGTTACGTCGGCATG  >  NC_002947/1247491‑1247779

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