Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,721,870	A→C	21.5%	V98G (GTC→GGC)	CJ019_01617 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	1,721,870	0	A	C	21.5%	41.0 / 2.5	36	V98G (GTC→GGC)	CJ019_01617	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/4);  new base C (4/4);  total (28/8)
Fisher's exact test for biased strand distribution p-value = 5.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCC‑GGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTAC  >  ACB122/1721737‑1721994                                                                                                                                      |                                                                                                                              gTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCCGCGGTCAGCATCTTCTTCCCCGGCTTCCTCGAc                                                                                                                              >  2:386036/1‑135 (MQ=255)    ttCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCa                                                                                                                           <  8:283023/135‑1 (MQ=255)     tCGCGGGTAAACCCATTCCTACAGAGACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAg                                                                                                                          >  4:164425/1‑135 (MQ=255)     tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGCCCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGCCCCAg                                                                                                                          >  3:346737/1‑135 (MQ=255)     tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCCGCGGTCAGCATCTTCGTCCCCGTCTTCCTCGCCCCAGc                                                                                                                         >  7:11484/1‑136 (MQ=255)     tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTGAGCATCTTCTTGCCCGTCTTCCTCGAccccg                                                                                                                          >  3:16118/1‑133 (MQ=255)     tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGcccccg                                                                                                                          >  1:319710/1‑133 (MQ=255)     tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAg                                                                                                                          >  3:307245/1‑135 (MQ=255)     tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGc                                                                                                                         >  7:73285/1‑136 (MQ=255)     tCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTGCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCCGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGCCCCAg                                                                                                                          >  3:373066/1‑135 (MQ=255)           gTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCAAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCCGCATCTTCTTCCCCGTCTTACTCGACCCCGCCAccc                                                                                                                    >  6:47003/1‑135 (MQ=255)           gTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGGTTCATGAAGTACGCCAGGTGATTTACTCATCACCGTCGGCGCTGAATCGGACCCAGCGGTCAGCATCTTAGTGCCCGGCATCCGAGACCCAACCAccc                                                                                                                    >  6:20542/1‑135 (MQ=255)             aaaCCCATTCCTACAGAAACGAGCAAGCCTGCCTTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCC‑Ggc                                                                                                                <  7:108327/136‑1 (MQ=255)              aaCCCATTCCTACATAAACGAGAAAGCCTGCGGTTTCATGATGTACGCCTGGTTATTTCCTCTTCACCGGCGGGGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCGGTCTTCCTCGACCCAGCCACCC‑Ggc                                                                                                                <  2:303898/135‑1 (MQ=255)                      ccTACAGAAAAGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCAACGTCGGCGCTGATGCGGACCCAGCCGTCACCATCTTCTTACCCCTCTTCCTCGAACCAGCCACCC‑GTCgcgcagcg                                                                                                        >  1:215584/1‑135 (MQ=255)                                  aGCACGCCTGCCGTTTCATGATGAACGCCCGGTGATTTCCTCTCCCCCGTCGGCGCTGCGGCGGCCCCAGCGGGCAGAATCTTCTCCCCGGTCTCCCTCCCCCCCCCCCCCC‑GGCGCGCCGCGCACTTCCACAccc                                                                                           <  8:394094/136‑1 (MQ=255)                                  aGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCCTCTTCTTCCCCGTCTTCCTCGACCCAGCCCCCC‑GGCGCGCCGCGCCCTTCCACAccc                                                                                           >  7:132183/1‑136 (MQ=255)                                      aGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCGTCTTCCCCGTCTTCCTCGGCCCAGCCACCC‑GGCGCGCAGCGCCCGTCCACACCCAgc                                                                                        >  4:192667/1‑135 (MQ=255)                                                    tgatgTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCCTCTTCTTCCCCGTCTTCCTCGCCCCAGCCACCC‑GGCGCGCCGCGCACGTCCACCCCCAGCGCGGCAAATgccgc                                                                          >  3:350665/1‑135 (MQ=255)                                                    tgatgTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCCGCCACCC‑GGCGCGCAGCGCACGTCCACACCCAGCGCGGCAAATgcggc                                                                          >  3:165174/1‑135 (MQ=255)                                                                   gATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCC‑GGCGCGCAGCGCACTTCCACACCCCGCGCGGCAAATGCGGCGCGGGCGCAGGCCAc                                                           >  1:209790/1‑135 (MQ=255)                                                                    aTTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCC‑GGCGCGCCGCGCCCTTCCACACCCCCCGCCGCAAATGCGGCGCGGGCGCAGGCAACg                                                          >  5:360197/1‑135 (MQ=255)                                                                    aTTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCC‑GGCGCGCAGCGCACTTACACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACg                                                          >  5:369268/1‑135 (MQ=255)                                                                    aTTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTAGACCCAGCCACCC‑GGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAAAGCGGCGCGGGCGCCGGCAACGt                                                         >  7:184095/1‑136 (MQ=255)                                                                    aTTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCATACCCGTCTTACTAGAACCAGCCACCC‑GGCGCGCAGACGAATTCCACCCCCACCGCGGCACATGCCGCACGGGCGCAGGCAACg                                                          >  2:161138/1‑135 (MQ=255)                                                                      ttCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCC‑CGCGCGCCGCGCACTTCCACACCCACCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTc                                                        >  4:23581/1‑135 (MQ=255)                                                                           ctTCACAGTCGGCGCTGAAGCGGACCCAGCGGTCAGCATCTTCTTCCACGTCTTCCTCGACCCAACCACCC‑GGCGCGCAGCGCACTTCAACACCCAGCGCGGCAAATGCGGCGCGGCCGCAGGCAACCTCATCCGc                                                  >  7:269973/1‑136 (MQ=255)                                                                              cACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCC‑GGCGCGCCGCGCACTTCCACACCCCGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGccc                                                >  3:109606/1‑135 (MQ=255)                                                                                 cGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCCGCCCCCC‑GGCGCGCAGCGCACTTCCACAACCAGCGCGGCTAAAGCGGCGCGGGCGCAGGGCACGGCATCCGGCCACg                                             >  5:293050/1‑135 (MQ=255)                                                                                      gcgcTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCC‑GGCGCGCAGCGCACTTCCACCCCCCGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCCCGGCGTc                                        >  5:88451/1‑135 (MQ=255)                                                                                          tggtGCGCACCCTGGGGCCTGCATGTTCCGCCCCGGCTGCC‑CGCCCCCCCCCCCCGGGCGCGCCCCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGt                                    <  5:305245/132‑1 (MQ=255)                                                                                             ggCGGACCCCGCGGTCAGCTTTTGCTGCCCCGCCTCCCGCGCCCCCCCCACCC‑GGGGCGCCCCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt                                <  8:132183/135‑1 (MQ=255)                                                                                              ggggACCCCGCGGACAGCATTTTCGTCCCCGTCTCCCCCCCCCCCGCCCCCC‑GGCGCGCCCCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCt                                <  6:71362/133‑1 (MQ=255)                                                                                                ggACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCC‑GGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCttt                              >  2:810/1‑135 (MQ=255)                                                                                                              ccATCTTCTGCCCCGGCAGCCTCCACCCCCCCCCCC‑GGCGCCCCGCCCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTg                <  6:153848/134‑1 (MQ=255)                                                                                                                           cGTCTTCCTCGACCCAGCCACCC‑GGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTa   >  5:153848/1‑135 (MQ=255)                                                                                                                            gTCTTCCTCGACCCAGCCACCC‑GGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTAc  >  6:134444/1‑135 (MQ=255)                                                                                                                                      |                                                                                                                              GTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCGTCTTCCTCGACCCAGCCACCC‑GGCGCGCAGCGCACTTCCACACCCAGCGCGGCAAATGCGGCGCGGGCGCAGGCAACGTCATCGGCCCACGGCGTCTGGTCGCTTTCCAGATACAGGCTGTTCCATTTCCCTAC  >  ACB122/1721737‑1721994

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