Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	2,094,817	T→G	44.9%	Q183P (CAG→CCG)	CJ019_01995 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	2,094,817	0	T	G	44.9%	2.1 / 8.2	35	Q183P (CAG→CCG)	CJ019_01995	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (7/10);  new base G (5/9);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATG  >  ACB143/2094707‑2094947                                                                                                                |                                                                                                                                   cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTTCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCCCGCGCGCCAGGCCCGGCAGGCGGTTTTTGAgggg                                                                                                           >  7:55640/1‑135 (MQ=255) cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCACGCGCCCCAGGCCGGGCTGGCCATTTTGCAGGGc                                                                                                           >  6:281877/1‑136 (MQ=255) cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCACGCCCGCCCGGCCGGGCGGGCGGGTTTTCAgggg                                                                                                           >  3:83533/1‑135 (MQ=255) cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCAGCTCTTCATACACCTTCGGCGCCGGCAGCTTGCGCAACTAC‑TGCACTCCGCGCGCGAAGGCCCAGCGCGCCGGGACCGGCTGGCCGTTTTgtcatgc                                                                                                           >  8:270247/1‑130 (MQ=255)      gCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGTCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCCTGCGCGCGGGGCCCGGGAGGGCGTTTTGACGGGGCTggt                                                                                                      >  2:186010/1‑134 (MQ=255)      gCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCCCGCGCGCCGGGGCGGGCGGGGCGGTTTGGAGGGcggggg                                                                                                      >  3:132095/1‑131 (MQ=255)              cAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGGGGGGCATTTTGGAGGGCCTGGTGGGCGAgg                                                                                              >  6:22150/1‑136 (MQ=255)                   gTTTCACATGATAGGTAAGGTGCTTTTCATACACGTTCGGCTCCGGGTGGTTGCTCAAATAC‑TCCACGCCGCGCTGGAAGGCCAGGGGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt                                                                                         <  2:179444/136‑1 (MQ=255)                   gTTTCACAAGATACTTCAGCTGCTCTACATACACCTTCGGCGCCGGCTGCTTGCCCAACTAC‑TCCACGCCGCGCGCGAAGGCCTGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt                                                                                         <  2:175446/136‑1 (MQ=255)                   gTTGCTCAAGATACGTCTGCTGCTCTTCATACACCTTCGGCGCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCGCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt                                                                                         <  2:41191/136‑1 (MQ=255)                    tttCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCACGCGCGCCTGGCCCGGCTGGCCGTTTTGCAGGGCCTGGTGGCCGAGggtgtt                                                                                         >  3:117463/1‑135 (MQ=255)                     ttgtcAATATCCGTCAGCTGCTCTTCAAACAGATTGGGGTCGGGGTGCTTTC‑AAAACACATGCCCCCGGCGCTAGAAGGCCAGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGttgttg                                                                                        <  6:41456/132‑1 (MQ=255)                      tctcAAGATACGGCACCTGCTATCAAAACACCTTCGGCTCCGGCTGGTTGCCAAAATAC‑TCCCCGCCGCGCGCGACGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTga                                                                                       <  8:88898/135‑1 (MQ=255)                         cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCAGGCCCGCCTGGCCCGGCGGGCCATTTTGAAGGGCCAGTTGGCCGAGGTTGTTgaaga                                                                                    >  8:96483/1‑135 (MQ=255)                                    aGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGACGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgc                                                                         <  8:269084/135‑1 (MQ=255)                                     gctgctCTTCATACACCTTCGGCTCCGGCTGCTTGCGAAAATAC‑TGCCCGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcag                                                                       <  8:156388/136‑1 (MQ=255)                                        gatctTCATACACGTTCCGCTCGGGCTGCTTGCGAAAATAC‑GGCCCGCGGCGCTCGAAGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag                                                                    <  1:93640/134‑1 (MQ=255)                                         ctctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag                                                                    <  3:149941/135‑1 (MQ=255)                                         ctctTCAAACACCTTCGGCTCCGGCTGCTGCCGAAAATAC‑TGCACGCCGCGCCCGAAGGCCTGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag                                                                    <  6:218162/135‑1 (MQ=255)                                            ttCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCACGCGCGCCTGGCCGGGCTGGCCTTTTTTCAGGGCCTGTTGGGCGAGGTTGTTGAAGAACACTATTTGCAGCAGCAcc                                                                >  6:242942/1‑136 (MQ=255)                                                      ttCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCAGGCGCGCCCGGCCGGGCGGGCCGGTTTTGAGGGCCGGGTGGGCGCGGTTGGTTAAGAAATCTATATGAAGGAGCCCCCCGATATGGc                                                      >  6:29980/1‑136 (MQ=255)                                                        cGGCGCGGGCTGGTTGCCAAAATAC‑CGCCCGCCGCGCTCGCCGGCCCGGCGGGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGc                                                      <  4:149403/134‑1 (MQ=255)                                                          gCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGGCGGGGTGGGCGTTTTTCAGGGCCGGTTGGGCGAGGTTGTTGTAGAACTCGATATGCAGCAGCACCAGGATATGGCGGAt                                                  >  2:169278/1‑136 (MQ=255)                                                            cccGGCGGCTTGCCAAAATAC‑TCCACGCCGCGCTCGAAGCCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAttt                                                <  3:223859/135‑1 (MQ=255)                                                                       gcaaaaaTAC‑CCCCCGCCGCGCGCCACGGCCAGCCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCaa                                     <  2:156009/132‑1 (MQ=255)                                                                                     ccgccgcGCTCGACGGCCAGGGGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGa                        <  3:52750/135‑1 (MQ=255)                                                                                      cgccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCGTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGa                        >  7:212710/1‑135 (MQ=255)                                                                                       gccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGGTGGTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATTCTTTGCCTCGACa                      >  1:171542/1‑136 (MQ=255)                                                                                           cgcTCGAAGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCgg                  <  8:157014/136‑1 (MQ=255)                                                                                            gcTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCGTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATTGCGGATTTCCTTTCCCCAATGCTTGGCCTCGACAGCggg                 >  5:166443/1‑136 (MQ=255)                                                                                             cgcGACGGCCCGGCGCGCGGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGc                <  3:9120/134‑1 (MQ=255)                                                                                                gAAGGCCAGGCGGGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag             <  6:205886/136‑1 (MQ=255)                                                                                                     cccggCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGcc        <  8:55640/133‑1 (MQ=255)                                                                                                          cgcgcCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCAt   <  4:15583/136‑1 (MQ=255)                                                                                                           gcgcCTGGCCGGGCTGGCCTTTTTTCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCACTGCTTGGCCTCGACAGCGGGCAAGCAGCCCTCATg  >  3:174777/1‑136 (MQ=255)                                                                                                                |                                                                                                                                   CATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATAC‑TGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATG  >  ACB143/2094707‑2094947

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