Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	456,874	T→C	25.5%	M425V (ATG→GTG)	iaaM ←	Tryptophan 2‑monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	456,874	0	T	C	25.5%	74.7 / 30.4	47	M425V (ATG→GTG)	iaaM	Tryptophan 2‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base T (15/20);  new base C (10/2);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 2.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

TCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGCGGTCGGTGAGGGTCATGCTCAGCAGGTCACGGCCGGTT  >  ACB143/456750‑457004                                                                                                                             |                                                                                                                                   tCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGt                                                                                                                         <  5:283825/136‑1 (MQ=255)   ggCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCg                                                                                                                       <  7:302649/136‑1 (MQ=255)       tCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCa                                                                                                                   >  6:126954/1‑136 (MQ=255)       tCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGACCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCa                                                                                                                   >  8:240103/1‑136 (MQ=255)        cccAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCAt                                                                                                                  >  7:50810/1‑136 (MQ=255)        cccAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCACCTTCAGCGCGTCGCTCAt                                                                                                                  >  8:232764/1‑136 (MQ=255)         ccAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCa                                                                                                                   <  1:217987/134‑1 (MQ=255)                     gctCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGAc                                                                                                     <  4:254808/134‑1 (MQ=255)                      aTCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACa                                                                                                    <  1:18731/136‑1 (MQ=255)                              gCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGAt                                                                                             >  1:35086/1‑135 (MQ=255)                              gCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATc                                                                                            <  2:201116/136‑1 (MQ=255)                               ccGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCa                                                                                           <  1:57726/136‑1 (MQ=255)                                             ggCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCACCTTCAGCGCGTCGCTAATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCg                                                                             >  4:108492/1‑136 (MQ=255)                                                 atatCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGc                                                                         <  5:172/136‑1 (MQ=255)                                                        tCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTc                                                                  <  8:199096/136‑1 (MQ=255)                                                        tCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTc                                                                  <  1:74309/136‑1 (MQ=255)                                                        tCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTc                                                                  <  5:37368/136‑1 (MQ=255)                                                        tCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCACCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTc                                                                  <  3:108492/136‑1 (MQ=255)                                                        tCTTCGGGTAGATCTTCTTCAGCGAATCCAGGGCAAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTc                                                                  <  4:237411/136‑1 (MQ=255)                                                            cGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGaa                                                               >  5:73291/1‑135 (MQ=255)                                                            cGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGaa                                                               >  3:310379/1‑135 (MQ=255)                                                             gggTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACa                                                             <  8:53967/136‑1 (MQ=255)                                                              ggTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACac                                                            >  4:255860/1‑135 (MQ=255)                                                               gTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAgg                                                           >  7:47151/1‑136 (MQ=255)                                                                     cttcttCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTg                                                     >  3:170495/1‑136 (MQ=255)                                                                     cttcttCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTg                                                     >  7:132236/1‑136 (MQ=255)                                                                     cttcttCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCACGGGCTTGTCGTTGCCGTTGTCGAAAAGGTAtgtg                                                     >  2:166768/1‑136 (MQ=255)                                                                           cAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCACCTTCAGCGCGTCGCTCATCCACGCGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAGCAGGTAAGTGCCGCgg                                               >  2:239684/1‑136 (MQ=255)                                                                            aGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCACCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCTCggg                                              >  2:148558/1‑136 (MQ=255)                                                                            aGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCACCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCggg                                              >  5:47856/1‑136 (MQ=255)                                                                             gcgcATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCCCCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGt                                             <  6:67227/136‑1 (MQ=255)                                                                                cATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAg                                          <  2:35086/136‑1 (MQ=255)                                                                                      gggCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGTTAAGTGCCGCGGGTGAGGcggtc                                    <  1:229123/136‑1 (MQ=255)                                                                                      gggCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGcggtc                                    <  5:130373/136‑1 (MQ=255)                                                                                      gggCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGcggtc                                    <  7:171780/136‑1 (MQ=255)                                                                                             ttGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGCGGTCGGTGAgg                             <  7:183650/136‑1 (MQ=255)                                                                                              tGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGCGGTCGGTGAggg                            >  4:101643/1‑136 (MQ=255)                                                                                              tGTACGCGCTTTTCCACCGGGTGCGGCAGCACCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGCGGTCGGTGAggg                            >  5:130355/1‑136 (MQ=255)                                                                                               gTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGCGGTCGGTGAGGGt                           >  6:220994/1‑136 (MQ=255)                                                                                                tACGCGCTTTTCCACCGGGTGCGGCAGCACCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGGAGATCACCGCGGTCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGCGGTCGGTGAGGGTc                          <  7:196755/136‑1 (MQ=255)                                                                                                tACGCGCTTTTCCACCGGGTGCGGCAGCACCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGCGGTCGGTGAGGGt                           >  5:190259/1‑135 (MQ=255)                                                                                                     gcTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGCGGTCGGTGAGGGTCATGCt                     >  7:62827/1‑136 (MQ=255)                                                                                                     gcTTTTCCACCGGGTGCGGCAGCACCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGCGGTCGGTGAGGGTCATGCt                     >  5:294400/1‑136 (MQ=255)                                                                                                             aCCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGCGGTCGGTGAGGGTCATGCTCAGCAGGt             >  3:308485/1‑136 (MQ=255)                                                                                                                    gCGGCAGCACCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGCGGTCGGTGAGGGTCATGCTCAGCAGGTCACGGcc      >  1:198409/1‑136 (MQ=255)                                                                                                                      ggCAGCACCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGCGGTCGGTGAGGGTCATGCTCAGCAGGTCACGGCCgg    >  1:150480/1‑136 (MQ=255)                                                                                                                        cagcaTCTTCAGCGCGTCTCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGCGGGCGGTGAGGGTCATGCTCAGCAGGTCACGGCCGGtt  <  3:73682/136‑1 (MQ=255)                                                                                                                             |                                                                                                                                   TCGGCCTCCCAGGAAACCGTGATCGGGTCGCCGATGATGTGCCCGGCGATATCGGTCTTCGGGTAGATCTTCTTCAGCGCATCCAGGGCCAGTTGTACGCGCTTTTCCACCGGGTGCGGCAGCATCTTCAGCGCGTCGCTCATCCACGAGTACGACAGGCAGATCACCCCGGGCTTGTCGTTGCCGTTGTCGAACAGGTAAGTGCCGCGGGTGAGGCGGTCGGTGAGGGTCATGCTCAGCAGGTCACGGCCGGTT  >  ACB143/456750‑457004

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: ‑