Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,284,434	G→A	43.5%	T368T (ACC→ACT)	argH ←	argininosuccinate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,284,434	0	G	A	43.5%	21.5 / 59.6	78	T368T (ACC→ACT)	argH	argininosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base G (6/38);  new base A (13/21);  total (19/59)
Fisher's exact test for biased strand distribution p-value = 1.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GCTTCACCAACAATATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTGGACGTTTCACCTGAATGCC  >  minE/2284366‑2284496                                                                     |                                                               gCTTCACCAACAATATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCg                                                                >  1:2597367/1‑69 (MQ=255) gCTTCACCAACAATATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCg                                                                >  1:1646807/1‑69 (MQ=255) gCTTCACCAACAATATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCg                                                                >  1:2990521/1‑69 (MQ=255) gCTTCACCAACAATATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCg                                                                >  1:2375766/1‑69 (MQ=255)       ccAACAATATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGtggc                                                           <  1:555585/68‑1 (MQ=255)       ccAACAATATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGtggc                                                           <  1:1835767/68‑1 (MQ=255)       ccAACAATATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGtggc                                                           <  1:831208/68‑1 (MQ=255)       ccAACAATATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGtggc                                                           <  1:153779/68‑1 (MQ=255)         aacaaTATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGtggcgtt                                                        <  1:739970/69‑1 (MQ=255)         aacaaTATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGtggcgtt                                                        <  1:1845686/69‑1 (MQ=255)         aacaaTATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGtggcgtt                                                        <  1:1108832/69‑1 (MQ=255)         aacaaTATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGtggcgtt                                                        <  1:2087134/69‑1 (MQ=255)          acaaTATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGtggcgtt                                                        <  1:829380/68‑1 (MQ=255)          acaaTATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGtggc                                                           <  1:2419452/65‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCTCCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:1095723/68‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:2530461/68‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:732713/68‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:2323765/68‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:519858/68‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:635815/68‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:2594347/68‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:2146481/68‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:1640982/68‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:802705/68‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:1467721/68‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:494677/68‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:3049074/68‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:913335/68‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:386437/68‑1 (MQ=255)                    gcgcCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAAccc                                              <  1:3304457/68‑1 (MQ=255)                          cgcgGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAACCctgctgcg                                       >  1:2685519/1‑67 (MQ=255)                          cgcgGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAACCctgctgcg                                       >  1:3341926/1‑67 (MQ=255)                          cgcgGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAACCctgctgcg                                       >  1:3357985/1‑67 (MQ=255)                          cgcgGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAACCctgctgcg                                       >  1:921012/1‑67 (MQ=255)                          cgcgGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAACCctgctgcg                                       >  1:2092886/1‑67 (MQ=255)                          cgcgGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAACCctgctgcg                                       >  1:1967606/1‑67 (MQ=255)                          cgcgGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAACCctgctgcg                                       >  1:1925731/1‑67 (MQ=255)                          cgcgGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAACCctgctgcg                                       >  1:1802613/1‑67 (MQ=255)                          cgcgGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAACCctgctgc                                        >  1:2026436/1‑67 (MQ=255)                          cgcgGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAACCctgctgc                                        >  1:3091339/1‑67 (MQ=255)                          cgcgGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAACCctgctgc                                        >  1:2474568/1‑67 (MQ=255)                          cgcgGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAACCctgctgc                                        >  1:371832/1‑67 (MQ=255)                          cgcgGAACGGTACGCCTTTCGCCACAAGATAATCCGCCAGTTCAGTAGCGTTGGCGTAACCctgctgcg                                       >  1:1382473/1‑67 (MQ=255)                             ggAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCctg                                            >  1:1106065/1‑61 (MQ=255)                             ggAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGc                                      <  1:3235632/67‑1 (MQ=255)                             ggAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGc                                      <  1:1421024/67‑1 (MQ=255)                             ggAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGc                                      <  1:1566820/67‑1 (MQ=255)                             ggAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGc                                      <  1:2159411/67‑1 (MQ=255)                             ggAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGc                                      <  1:279036/67‑1 (MQ=255)                             ggAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGc                                      <  1:330337/67‑1 (MQ=255)                             ggAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGc                                      <  1:2401652/67‑1 (MQ=255)                             ggAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGc                                      <  1:60557/67‑1 (MQ=255)                             ggAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCg                                    >  1:1881631/1‑69 (MQ=255)                             ggAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCg                                    >  1:1885438/1‑69 (MQ=255)                                        ccTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCaa                         <  1:2824732/69‑1 (MQ=255)                                        ccTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCaa                         >  1:458066/1‑69 (MQ=255)                                        ccTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCaa                         <  1:508610/69‑1 (MQ=255)                                        ccTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCaa                         <  1:1727011/69‑1 (MQ=255)                                              gCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTg                     <  1:674541/67‑1 (MQ=255)                                              gCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTg                     <  1:1022429/67‑1 (MQ=255)                                              gCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTg                     <  1:968564/67‑1 (MQ=255)                                              gCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTg                     <  1:3261859/67‑1 (MQ=255)                                              gCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTg                     <  1:1051026/67‑1 (MQ=255)                                              gCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTg                     <  1:116901/67‑1 (MQ=255)                                              gCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTg                     <  1:1285351/67‑1 (MQ=255)                                              gCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTg                     <  1:2215210/67‑1 (MQ=255)                                              gCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTg                     <  1:2565270/67‑1 (MQ=255)                                              gCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTg                     <  1:2891844/67‑1 (MQ=255)                                              gCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTg                     <  1:3411299/67‑1 (MQ=255)                                                    aGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTGGACGTTTc             >  1:1518541/1‑69 (MQ=255)                                                    aGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTGGACGTTTc             >  1:912243/1‑69 (MQ=255)                                                         aTCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTGGACGTTTCACCTg        <  1:2063656/69‑1 (MQ=255)                                                         aTCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTGGACGTTTCACCTg        <  1:2034143/69‑1 (MQ=255)                                                         aTCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTGGACGTTTCACCTg        <  1:1393491/69‑1 (MQ=255)                                                         aTCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTGGACGTTTCACCTg        <  1:3165543/69‑1 (MQ=255)                                                         aTCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTGGACGTTTCACCTg        <  1:846079/69‑1 (MQ=255)                                                               cAGTTCAGTAGCGTTGGCGTAACCCTGCTGCGCCGCTTCCTGGCAACGCGGACGTTTCACCTGAATGcc  <  1:225554/69‑1 (MQ=255)                                                               cAGTTCAGTAGCGTTGGCGTAACCCTGCTGCGCCGCTTCCTGGCAACGCGGACGTTTCACCTGAATGcc  <  1:3036713/69‑1 (MQ=255)                                                               cAGTTCAGTAGCGTTGGCGTAACCCTGCTGCGCCGCTTCCTGGCAACGCGGACGTTTCACCTGAATGcc  <  1:1892476/69‑1 (MQ=255)                                                               cAGTTCAGTAGCGTTGGCGTAACCCTGCTGCGCCGCTTCCTGGCAACGCGGACGTTTCACCTGAATGcc  <  1:18671/69‑1 (MQ=255)                                                               cAGTTCAGTAGCGTTGGCGTAACCCTGCTGCGCCGCTTCCTGGCAACGCGGACGTTTCACCTGAATGcc  <  1:3248874/69‑1 (MQ=255)                                                                  ttCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTg                     <  1:12371/47‑1 (MQ=255)                                                                     |                                                               GCTTCACCAACAATATGGTGCGCCTCGCGGAACGGTACGCCTTTCGCCACCAGATAATCCGCCAGTTCGGTGGCGTTGGCGTAACCCTGCTGAGCCGCTTCCTGGCAACGTGGACGTTTCACCTGAATGCC  >  minE/2284366‑2284496

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