Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	3,772,505	A→G	F313F (TTT→TTC)	LY180_18030 ←	glutamate decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	3,772,505	0	A	G	94.7%	24.2 / ‑6.0	19	F313F (TTT→TTC)	LY180_18030	glutamate decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (5/13);  minor base C (0/1);  total (5/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

AGGAATTCATAGTACTGTGCAATTACCTGACCCGCCGGGCGGGAGAAGTTGATGGCAAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTA  >  CP006584/3772449‑3772799                                                         |                                                                                                                                                                                                                                                                                                       aGGAATTCCTGGTACTGTGAAATTACCTGACCCGCCGGGCGGGAGAAGTTGAGGGCGAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTAt                                                       <  2:544074/298‑1 (MQ=2)        cATAGTACTGTGCAATTACCTGACCCGCCGGGCGGGAGAAGTTGATGGCGAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTg                                                <  1:47757/298‑1 (MQ=14)          tAGTACTGTGCAATTACCTGACCCGCCGGGCGGGAGAAGTTGATGGCGAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGa                                              <  1:443416/298‑1 (MQ=14)                 gtgCAATTACCTGACCCGCCGGGCGGGAGAAGTTGATGGCGAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATc                                       >  1:142605/1‑298 (MQ=14)                    cAATTACCTGACCCGCCGGGCGGGAGAAGTTGATGGCGAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAg                                    <  1:565555/298‑1 (MQ=14)                      aTTACCTGACCCGCCGGGCGGGAGAAGTTGATGGCGAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAgcg                                  >  1:237165/1‑298 (MQ=14)                         aCCTGACCCGCCGGGCGGGAGAAGTTGATGGCGAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCAt                               <  1:751661/298‑1 (MQ=14)                            ggACCCGCCGGGCGGGGGAAGTTGGTGGCGAAAGTGCCAAGTTGACCGCCCAGGTAGTCAACTTGGACCACCAGTTCCTGCGGCAGGGCTTCTTCGGCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGt                            <  2:452359/297‑1 (MQ=11)                              aCCCGCCGGGCGGGAGAAGTTGATGGCGAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGc                          <  1:537941/298‑1 (MQ=14)                                     gggcgggAGAAGTTGATGGCGAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGtt                   >  1:494816/1‑298 (MQ=14)                                      ggcgggAGAAGTTGATGGCGAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtg                  >  2:318804/1‑298 (MQ=14)                                       gcaggAGAAGTTGATGGCGAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTtgt                 <  1:376115/295‑1 (MQ=14)                                          ggAGAAGTTGATGGCGAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTggg              <  1:651755/298‑1 (MQ=14)                                           gagaAGTTGATGGCGAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGGGCCAGGAAGCCACCCCTGGCAGCGCCGAGGGGCATGTCGGTTTCGATCCCGGTATCGGCCTGGAATTTATCCAGCGCACCGTTCCGCGGGTGTgggc             >  2:321757/1‑297 (MQ=11)                                             gaAGTTGATGGCGAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAAc           <  1:552466/298‑1 (MQ=14)                                               aGTTGATGGCGAAAGAACCAAGTTGACGACCCGGGTAGTCAACGTTGAACACCAGTGCCTGCGGCAGCGCTTCTTCGTCACGCCCGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGGCTGAAGCACTGGTCGATTGCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTc         <  2:159447/298‑1 (MQ=11)                                                 gtggcggCGAAAGGACCAATTGGACCCCCCAGGTGGTCAACGTTGACGCCCATCTGCTCCGGTCGCCCTGCTTGGTCCGGCAAGATAACCCCGCCGTCGCCCCGCTGAGACAGACGGAATTGAGGGACTGAAGTATGGAGTGATGTAACACGCGGAGGGGGGAGGTCCCAAGCGAGATCCGGGGCGACGAACGGTGCTAGGAAGCCACCGCTGGCAGCTTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCAt       <  2:634390/293‑1 (MQ=2)                                                    aTGGCGAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGt    <  1:137412/298‑1 (MQ=14)                                                      ggCCAAAGTACCAATGTAACCAACCAGGGTGTCAACGTGTAACACCAGTACCTGCGGCAGCGCCTCGTCGCCACGCCAGGTAACCCAGCCGCAGCCCAGCGGGGCCAGACCGGATTTATGGCCTGAAGCACGGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTa  <  2:374014/298‑1 (MQ=11)                                                         |                                                                                                                                                                                                                                                                                                       AGGAATTCATAGTACTGTGCAATTACCTGACCCGCCGGGCGGGAGAAGTTGATGGCAAAAGTACCAATTTGACCACCCAGGTAGTCAACGTTGAACACCAGTTCCTGCGGCAGCGCTTCTTCGTCACGCCAGATAACCCAGCCGCAGCCCAGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTATCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTA  >  CP006584/3772449‑3772799

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