Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	4,342,664	T→G	34.7%	Q17P (CAG→CCG)	tsdA ←	Thiosulfate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	4,342,664	0	T	G	34.7%	22.8 / 2.7	32	Q17P (CAG→CCG)	tsdA	Thiosulfate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (6/13);  new base G (9/3);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 2.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGCCTTCACGCACCATTTTGCCAAAGGCATTATCGGGCAACTCACGCTCAGCAGGCGGCACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGCCTGGGCGCTGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑CAGGGTTTTCATGGCGTGACCTCCTGGCTGGCGGGGCGGGCGAAGAACTCGGCAACGGCCTTGATCTCCTCATCGGTCATGGCCTTGGCCACGC  >  ACB122/4342531‑4342795                                                                                                                                      |                                                                                                                                     gTGGCCTTCACGCACCATTTTGCCAAAGGCATTATCGGGCAACTCACGCTCAGCAGGCGGCACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGCCTgg                                                                                                                                    <  8:364841/136‑1 (MQ=255)      cTTCACGCACCATTTTGCCAAAGGCATTATCGGGCAACTCACGCTCAGCAGGCGGCACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGCCAGGgcgc                                                                                                                                >  5:168324/1‑135 (MQ=255)           cGCACCATTTTGCCAAAGGCATTATCGGGCAACTCACGCGCAGCCGGCGGCACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCACCTGCGACTGGGCTTCCATGGCGATCGGCGGGGCCTGGGCGCTGccc                                                                                                                           <  6:178006/135‑1 (MQ=255)            gCACCATTTTGCCAAAGGCATTATCGGGCAACTCACGCTCAGCAGGCGGCACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGGCTGGGCGCTGccccc                                                                                                                         >  8:93887/1‑134 (MQ=255)             cACCATTTTGCCAAAGGCATTATCGGGCAACTCACGCTCAGCAGGCGGCACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGGCGGGGGGCTGGCCCCg                                                                                                                        >  7:97362/1‑136 (MQ=255)                   tttGCCAAAGGCATTATCGGGCAACTCACGCTCAGCAGGCGGCACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGCCTGGGCGCTGCCCCCCGCcac                                                                                                                   >  2:253606/1‑134 (MQ=255)                     tGCCAAAGGCATTATCGGGCATCTCACGCTCAGCAGGCGGCACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGCCTGGGCGCTGGCCACCGCcagcag                                                                                                                >  7:395522/1‑136 (MQ=255)                                      ggCAACGCCCGCTCCGCCGGGGGCACGAACTGGGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGGGTGGCCTGGGCGCTGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑Ca                                                                                              <  4:209694/136‑1 (MQ=255)                                            tCACGCGCCGCAGGCGGCACGAACTGCGGTGCGGCCGGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGCCTGGGCGCTGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑CAGGGtt                                                                                         <  1:246683/135‑1 (MQ=255)                                                gCTCAGCAGGCGGCACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGCCGGGGCGCGGGCCCCGGCCCGCCCCAACAAGGCCGGGAT‑CACGGTTTTGAt                                                                                    >  7:27999/1‑136 (MQ=255)                                                    agcagGCGGCACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGGCGGGGCGCGGCCCCCGGCCCGCAGCCACCACGCGGGGAT‑AAGGGGTTTGATGGc                                                                                 >  1:271272/1‑135 (MQ=255)                                                    agcagGCGGCACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGCCTGGGCGCTGCCCCCGGCCAGCCGCCACAAGGCGGGCAT‑CAGGGGTTTCATGGCg                                                                                >  4:161145/1‑136 (MQ=255)                                                    agcagGCGGCACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGCCTGGGCGCTGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑CAGGGTTTTCATGGc                                                                                 <  6:168324/135‑1 (MQ=255)                                                        ggcggcACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCAATGGCGATCGGCGTGCCCGGTGCCGTGCCACCCGCCCCCAGCAACAACGCGGGCAA‑AAGGGGTTACATTGccgta                                                                             >  6:103738/1‑131 (MQ=255)                                                         gcggcACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGCCGGGGGGCGGGCCCCGGCCCGCCGCCACAAAGCGGGGAT‑CCTGGTTTTGATGGCGTggc                                                                            >  6:345179/1‑133 (MQ=255)                                                             cACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGGCGGGGCGCTGCCCCCCGCCCGCAGCCACAATGCGGGCAT‑CTGGGTGTTCATGGGGTGcccg                                                                          >  6:229319/1‑132 (MQ=255)                                                                         tgctgcCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGGCGGGGCGCTGCCCCCCGGCAGCACCCACCAGGGGGGGTTGCAGGG‑TTTAATGGCGTGGCCACCTGGCgggggggg                                                            >  1:310080/1‑135 (MQ=255)                                                                              cctgcACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGCCTGGGCGCTGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑CAGGGTTTTCATGGCGTGACCTCCTGGCTGGCGgggcgggc                                                       >  1:33423/1‑135 (MQ=255)                                                                                        gggAACCTTCAACTGCGACTTGTCTTCCATGGCGTGGGGGGTGGCCTGGGCGCTGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑CAGGGTTTTCATGGCGTGACCTCCTGGCTGGCGGGGCTGGCGAAGAACTCgg                                            <  7:188386/136‑1 (MQ=255)                                                                                          gAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGCCTGGGCGCGGCCCCCCGCCACCAGCAACAACGCGGGAAT‑CACGGTTTTCATGGCGTGACCACCAGGCTGGCGGGGCGGGGGCAGAAcaaggc                                           >  7:260664/1‑130 (MQ=255)                                                                                           aaCCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGCCGGGGCGCTGCCCCCCGCCAGCCCCCACAAGGCGGGAAA‑AAGGGGGTTCATGGCGGGCACGCCAGGCTGGCCGGGCGGGCGCAGCACGCGGCaa                                         >  4:140953/1‑136 (MQ=255)                                                                                                   ccTGCGACTGGGCTTCCTTGGCGATCGGCGGGGCGGGGGCGCGGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑CAGGGTTTTCATGGCGTGACCTCCTGGCTGGCGGGGCGGGCGAAGAACTCGGCAACGGCCtt                                  <  5:292491/134‑1 (MQ=255)                                                                                                          cTGGTCTTCCATGGCGATCGGCGTGGGCGGGGCGCAGCCCCCCGCCACCACCAAAAAGGCGGGAAT‑AATGGGTTTGATTGCGTGACCGCCAGGCTGGCGGGGCGGGCGAAGAAAAAGGCAACGGCCGTGAtctac                           >  5:273050/1‑133 (MQ=255)                                                                                                          cTGGTCTTCCATGGCGATCGGCGTGGCCTGGGCGCTGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑CAGGGTCTTCATGGCGTGACCTCCTGGCTGGCGGGGCGGGCGAAGAACTCGGCAACGGCCTTGATctcct                          <  3:104649/136‑1 (MQ=255)                                                                                                          cTGGTCTTCCATGGCGATCGGCGGGGCCTGGGCGCTGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑CAGGGTTTTCATGGCGTGACCTCCTGGCTGGCGGGGCGGGCGAAGAACTCGGCAACGGCCTCGATctcct                          <  4:376437/136‑1 (MQ=255)                                                                                                           tGGTCTTCCATGGCGATCGGCGTGGCCTGGGCGCTGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑CAGGGTTTTCATGGCGTGACCTCCTGGCTGGCGGGGCGGGCGAAGAACTCGGCAACGGCCTTGATctcct                          <  2:33423/135‑1 (MQ=255)                                                                                                           tGGTCTTCCATGGCGATCGGCGTGGCCTGGGCGCTGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑CAGGGTTTTCATGGCGTGACCTCCTGGCTGGCGGGGCGGGCGAAGAACTCGGCAACGGCCTTGATctcct                          <  6:273050/135‑1 (MQ=255)                                                                                                           tGGTCTTCCATGGAGATCGGCGTGGCCTGGGCGCTGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑CAGGGTTTTCATGGCGTGACCTCCTGGCTGGCGGGGCGGGCGAAGAACTCGGCAACGGCCTTGATctcct                          <  1:211109/135‑1 (MQ=255)                                                                                                                  ccATGGCGATCGGCGTGGCCTGGGCGCTGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑CAGGGTTTTCATGGCGTGACCTCCTGGCTGGCGGGGCGGGCGAAGAACTCGGCAACGGCCTTGATCTCCTCATCgg                    <  3:16829/134‑1 (MQ=255)                                                                                                                    aTGGCGATCGGCGTGGCCTGGGCGCTGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑CAGGGTTTTCATGGCGTGACCTCCTGGCTGGCGGGGCGGGCGAAGAACTCGGCAACGGCCTTGATCTCCTCATCGGTCAt                <  4:233538/136‑1 (MQ=255)                                                                                                                               cGTGGCCGGGGCGCTGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑CAGGGTTTTCATGGCGTGACCTCCTGGCTGGCGGGGCGGGCGAAGAACTCGGCAACGGCCTTGATCTCCTCATCGGTCATGGCCTTGGCCa     <  3:302186/136‑1 (MQ=255)                                                                                                                                  ggCCGGGGCGCTGCCCACGGCCAGCAGCAACTAGGCGGGACT‑CAGGGTTTTCATGGCGTGACCTCCTGGCTGGCGGGGCGGGCGAAGAACTCGGCAACGGCCTTGATCTCCTCATCGGTCATGGCCTTGGCCACGc  <  7:370715/136‑1 (MQ=255)                                                                                                                                      |                                                                                                                                     GTGGCCTTCACGCACCATTTTGCCAAAGGCATTATCGGGCAACTCACGCTCAGCAGGCGGCACGAACTGCGGTGCTGCCTGCACACCGGGAACCTTCAACTGCGACTGGTCTTCCATGGCGATCGGCGTGGCCTGGGCGCTGCCCACGGCCAGCAGCAACAAGGCGGGAAT‑CAGGGTTTTCATGGCGTGACCTCCTGGCTGGCGGGGCGGGCGAAGAACTCGGCAACGGCCTTGATCTCCTCATCGGTCATGGCCTTGGCCACGC  >  ACB122/4342531‑4342795

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