Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	3,761,329	A→C	42.6%	E96A (GAG→GCG)	CJ019_03574 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,761,329	0	A	C	42.6%	22.3 / 5.3	40	E96A (GAG→GCG)	CJ019_03574	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (8/13);  new base C (14/4);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 2.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCGG  >  ACB122/3761200‑3761462                                                                                                                                   |                                                                                                                                      gCAGGGGCCATGGTGAGCAAGACTGCAGCGACCCGACCCTGTCGGCCGGCGGGGACGGCGGGAGCGCCAGCCTGCTCCGCGCCACACTGCGGGCGGT‑CTTTTGCGGCCTGCTGACCCTGCGCGCCGGCGAGCCgg                                                                                                                                  <  4:143424/135‑1 (MQ=255)  cATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTGCTTTTTCGGCCTGCTGCCCCTGCACGCCGGCGCGCCgg                                                                                                                                  >  2:222878/1‑135 (MQ=255)  cATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTGCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCgg                                                                                                                                  >  2:218977/1‑135 (MQ=255)               gAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCACCAGCCTGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCCGGCAGCCCCTTAcc                                                                                                                     >  8:91530/1‑134 (MQ=255)                   aaGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCCGGCCGCCCCGTACcccacc                                                                                                                >  3:21218/1‑135 (MQ=255)                    aGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCGGGCAGCCCCGTAcccaccca                                                                                                               >  6:296418/1‑135 (MQ=255)                    aGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCGGCCAGCCCCTTACccaacca                                                                                                               >  4:115123/1‑135 (MQ=255)                                  cGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGT‑GTTTTTCGGCCTGCTGACCCTGCGCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGcc                                                                                                 <  8:113524/135‑1 (MQ=255)                                   gACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCCGCGCGCCCGGCAGCCCCGTTCCCCACCCACCCccgccccgcca                                                                                                >  2:151754/1‑134 (MQ=255)                                   gACCCGTTCTGCCGGCTGTGACTGCGTGATCGACAGCCTGTTCAGCGCAAAACTGCTGGCCGG‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCTAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACgac                                                                                                 <  6:83358/134‑3 (MQ=255)                                       cTGGCTGCCGGCGGGGACTGCGTGAGCGCCAGCCGGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGccc                                                                                            <  2:178077/135‑1 (MQ=255)                                            tGCCGGCGGTGACGGGGTGAGCGCCAGCCTGTGCCGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCGCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAagca                                                                                       <  5:100723/135‑1 (MQ=255)                                                   ggTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCGCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTgg                                                                                <  6:88394/135‑1 (MQ=255)                                                    gggACTGCGTGAGCGCCAGCCTGTTCCGCGCCACACTGCGGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCGCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGc                                                                               <  7:135009/133‑1 (MQ=255)                                                    gTGACTGGGGGCCCGCCAGCCTGTCCCGCCCCACACTGCGGGGCGG‑CTTTTTCGCCCGGCGGCCCCGGCGCGCCGGCGCGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGc                                                                               <  7:159707/135‑1 (MQ=255)                                                       aCTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCGGGCCGCCCCGTACCCACCCAACCTCCGCCCCGCCTGCCCCCGCACCGCGGGcccc                                                                            >  3:34465/1‑133 (MQ=255)                                                         tGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCGGGGAGCCCCGTACGCCACCAACCACCGCCCCGCCGGCCCAACCCGCGCGGGccccgg                                                                          >  7:293151/1‑131 (MQ=255)                                                         tGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCGGGCCGCCCCGTTCCCAACCACCCACCGCCCCGCCTGCCCCCCCAGCCCGGGCCAcag                                                                          >  3:148691/1‑133 (MQ=255)                                                              gATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGCCCCTGCTCGCCGGCGCGCCGGGCCGCCCCGTGCCCCCCCCCCCTCCGCCCCCCCGGCCCAGGCCGCGCGGGCCCCGGGTcca                                                                     >  5:706/1‑133 (MQ=255)                                                              gATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTAGGCCTGCTGACCCTGCTCGCCGGCGTGCCGGGCCGCGCCTTACCCACCCAACCTCCGCCCCGCCTGCACAAGCAGCCGAGGCCACTgggcca                                                                     >  6:117664/1‑131 (MQ=255)                                                                gcGCCAGCCTGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAAt                                                                    <  2:215029/133‑1 (MQ=255)                                                                        cTGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCCGCCCCGGACCCCACCAACCTCCGCCACGCCTGCCCAAGCCGCGCTGGCCCCTGGTCAATCCCCGCGcc                                                           >  3:296902/1‑135 (MQ=255)                                                                              aGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCCGGCCGCCCCGTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCGCTGGTCCATCCCCCCGCCGCGccc                                                     >  3:120091/1‑135 (MQ=255)                                                                              aGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCCGCGCGCCGGGCCCCCCCGGACCCACCCACCCACCCCCCCCCCGGCCCCAGCCCCGCTGGGCAGTGGGCAATCCCCCCGCCGCAccc                                                     >  1:97030/1‑135 (MQ=255)                                                                                   cacaCTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCCCTGGTCAATCCCCGCGCCTCGCCcgcttc                                                >  4:88529/1‑135 (MQ=255)                                                                                      cctggcggCGGT‑CTTTTGCGGCCGGGGCACCCGGCGCGCCGGCCCGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttcgct                                             <  1:151754/129‑1 (MQ=255)                                                                                       ctgcggGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttcgctt                                            <  2:217836/135‑1 (MQ=255)                                                                                             gCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCGGGCAGCCCCTTACCCAACCACCCTCCGCCACGCCTGCCCAAGCACAGCTGGCCACTCGTCAATCCCCGCGGCTCTCCCGCTGCGCTTCTggcc                                      >  4:159968/1‑134 (MQ=255)                                                                                             gCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCGCTCCCGCTTCGCGTCTcgcg                                      >  6:144325/1‑135 (MQ=255)                                                                                                   cTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCACCCTCCGCCACGCCTGCCCAAGCCGCGCTGGCCCCTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTcgg                                 >  4:119630/1‑133 (MQ=255)                                                                                                    ttttgCGGCCGGCTGACCCTGCGCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGa                                <  1:218977/135‑1 (MQ=255)                                                                                                       tgcGGCCTGCTGACCCTGCGCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGa                                <  1:222878/130‑1 (MQ=255)                                                                                                          ggCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCg                           <  6:67263/134‑1 (MQ=255)                                                                                                                         cTCGCCGGCGCGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTg                                 <  5:295809/113‑1 (MQ=255)                                                                                                                         cTCGCCGGCGCGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTg                                 <  5:291815/113‑1 (MQ=255)                                                                                                                         cTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTg                                 >  6:291815/1‑113 (MQ=255)                                                                                                                         cTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTg                                 >  6:295809/1‑113 (MQ=255)                                                                                                                             ccGGCGCGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCCCTGGTCAATCCCCGCGCCTCTCCCGCTTCGCGTCTCGCGGTCTGACATCGCCCCGTCCCCGTgccgcc         >  1:74465/1‑133 (MQ=255)                                                                                                                             ccGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCCCGCCTGCCCCAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTg       >  6:157393/1‑135 (MQ=255)                                                                                                                                  gAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCgg  <  5:156151/135‑1 (MQ=255)                                                                                                                                   |                                                                                                                                      GCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGT‑CTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCGG  >  ACB122/3761200‑3761462

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