Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,787,133	A→G	21.0%	D347G (GAC→GGC)	gabT →	4‑aminobutyrate aminotransferase, PLP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,787,133	0	A	G	21.0%	33.0 / 2.6	19	D347G (GAC→GGC)	gabT	4‑aminobutyrate aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (10/5);  new base G (2/2);  total (12/7)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GTAACCCGATTGCCTGCGTGGCTGCGCTGGAAGTGTTGAAGGTGTTTGAGCAGGAAAATCTGCTGCAAAAAGCCAACGATCTGGGGCAGAAGTTGAAAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGTACGCGGGCTGGGGGCGATGATCGCCATTGAGCTGTTTGAAGACGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTATTACA  >  CP009273/2786992‑2787276                                                                                                                                              |                                                                                                                                                gTAACCCGATTGCCCGCGTGGCTGCGCTGGAAGTGTTGAAGGTGTTTGAGCAGGAAAATCTGCTGCAAAAAGCCAACGATCTGGGGCAGAAGTTGAAAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGTAcgc                                                                                                                                          <  1:168932/149‑1 (MQ=255)         aTTGCCTGCGTGGCTGCGCTGGAAGTGTTGAAGGTGTTTGAGCAGGAAAATCTGCTGCAAAAAGCCAACGATCTGGGGCAGAAGTTGAAAGACTGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGTACGCgg                                                                                                                                        <  2:255873/143‑1 (MQ=255)           tGCCTGCGTGGCTGCGCTGGAAGTGTTGAAGGTGTTTGAGCAGGAAAATCTGCTGCAAAAAGCCAACGATCTGGGGCAGAAGTTGAAAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGTACGCGGGc                                                                                                                                      >  2:101109/1‑143 (MQ=255)                           cTGGAAGTGTTGAAGGTGTTTGAGCAGGAAAATCTGCTGCAAAAAGCCAACGATCTGGGGCAGAAATGGATAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGTACGCGGGCTGGGGGCGATGATCGCCATTgg                                                                                                                >  1:145578/1‑148 (MQ=255)                              gAAGTGTTGAAGGTGTTTGAGCAGGAAAATCTGCTGCAAAAAGCCAACGATCTGGGGCAGAAGTTGAAAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGGCGTACGCGGGCTGGGGGCGATGATCGCCATTGAGCt                                                                                                             >  1:121267/1‑149 (MQ=255)                              gAAGTGTTGAAGGTGTTTGAGCAGGAAAATCTGCTGCAAAAAGCCAACGATCTGGGGCAGAACTTGAAAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGGCGTATGCGGGCTGGGGGCGATGATCGCCATTGAGCt                                                                                                             >  1:222730/1‑149 (MQ=255)                               aaGTGTTGAAGGTGTTTGAGCAGGAAAATCTGCTGCAAAAAGCCAACGATCTGGGGCAGAAGTTGAAAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGTACGCGGGCTGGGGGCGATGATCGCCATTGAGCTg                                                                                                            >  1:46481/1‑149 (MQ=255)                                       aaGGTGTTTGAGCAGGAAAATCTGCTGCAAAAAGCCAACGATCTGGGGCAGAAGTTGAAAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGTACGCGGGCTGGGGGCGATGATCGCCATTGAGCTGTTTGAAGa                                                                                                    >  1:227599/1‑149 (MQ=255)                                            gtTTGAGCAGGAAAATCTGCTGCAAAAAGCCAACGATCTGGGGCAGAAGTTGAAAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGGACGCGGGCTGGGGGCGATGATCGCCATTGAGCTGTTTGacg                                                                                                     >  2:187120/1‑141 (MQ=255)                                                 aGCAGGAAAATCTGCTGCAAAAAGCCAACGATCTGGGGCAGAAGTTGAAAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGTACGCGGGCTGGGGGCGATGATCGCCATTGAGCTGTTTGAAGACGGCGATcac                                                                                          >  1:290845/1‑149 (MQ=255)                                                   cAGGAAACGCTGCTGCAACAAGCCAACGTTCTGGGGCAGATGTTGGTCTTCGGAGGGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGTACGCGGGCTGGGGGCGATGATCGCCATTGAGCTGTTTGAAGATGGCGATCacaa                                                                                        <  1:76020/149‑1 (MQ=255)                                                      gAAAATCTGCTGCAAAAAGCCAACGATCTGGGGCAGAAGTTGAAAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGTACGCGGGCTGGGGGCGATGATCGCCATTGAGCTGTTTGAAGACGGCGATCAcaacaa                                                                                     >  1:66982/1‑149 (MQ=255)                                                          aTCTGCTGCAAAAAGCCAACGATCGGGGGCAGAAGTTGAAAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGTACGCGGGCTGGGGGCGATGATCGCCATTGAGCTGTTTGAAGACGGCGATCACAACAAGCCg                                                                                 <  1:97285/149‑1 (MQ=255)                                                           tCTGCTGCAAAAAGCCAACGATCTGGGGCAGAAGTTGAAAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGTACGCGGGCTGGGGGCGATGATCGCCATTGAGCTGTTTGAAGACGGCGATCACAACAAGCCgg                                                                                <  1:311642/149‑1 (MQ=255)                                                                    aaaaGCCAACGATCTGGGGCAGAAGTTGAAAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGTACGCGGGCTGGGGGCGATGATCGCCATTGAGCTGTTTGAAGACGGCGATCACAACAAGCCGGACg                                                                             >  2:220297/1‑143 (MQ=255)                                                                                    gggCAGAAGTTGAAAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGTACGCGGGCTGGGGGCGATGATCGCCATTGAGCTGTTTGAAGACGGCGATCGCAACAAGCCGGACGCCAAACTCACCGCCGAGATCGt                                                       >  1:288522/1‑149 (MQ=255)                                                                                                                      ccGAAAAACAGCCGGAGATCGGCGGCGTACGCGGGCTGGGGGCGATGATCGCCATTGAGCTGTTTGAAGACGGCAATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGAtt                           <  2:197730/143‑1 (MQ=255)                                                                                                                      ccGAAAAACACCCGGAGATCGGCGGCGTACGCGGGCTGGGGGCGATGATCGCCATTGAGCTGTTTGAAGACGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGAtt                           <  2:106702/143‑1 (MQ=255)                                                                                                                                         cGGCGACGTACGCGGGCTGGGGGCGATGATCGCCATTGAGCTGTTTGAAGACGGCGATCCCAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTATTACa  >  1:294423/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                GTAACCCGATTGCCTGCGTGGCTGCGCTGGAAGTGTTGAAGGTGTTTGAGCAGGAAAATCTGCTGCAAAAAGCCAACGATCTGGGGCAGAAGTTGAAAGACGGATTGCTGGCGATAGCCGAAAAACACCCGGAGATCGGCGACGTACGCGGGCTGGGGGCGATGATCGCCATTGAGCTGTTTGAAGACGGCGATCACAACAAGCCGGACGCCAAACTCACCGCCGAGATCGTGGCTCGCGCCCGCGATAAAGGCCTGATTCTTCTCTCCTGCGGCCCGTATTACA  >  CP009273/2786992‑2787276

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