Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	3,986,622	T→G	35.3%	D31A (GAT→GCT)	CJ019_03763 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,986,622	0	T	G	35.3%	20.8 / 4.2	41	D31A (GAT→GCT)	CJ019_03763	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/8);  new base G (3/9);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 3.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAGGTTGACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGG  >  ACB122/3986488‑3986750                                                                                                                                       |                                                                                                                                 ggTCGGTTGACGCCATCAGCCCGTGCGCCCGCAACTGATCCGGTTACACGGCAAGGGGGGCCAGCATCTGGTCCAGTGGCTCCTGGGAGAACACCGGGGCCTTGGCTACGAGCGTGGCAGCGGACGGCGCCGGAt                                                                                                                                  <  8:122831/135‑1 (MQ=255)         gACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACCCCGGGTCCTTGGCTACGAGCGTGGCCGCGGACGGCGCCCGCTCCGCcgcg                                                                                                                          >  6:256163/1‑135 (MQ=255)         gACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCGGGGGGCGGcgcg                                                                                                                          >  5:335800/1‑135 (MQ=255)         gACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGGTCCGGcgcg                                                                                                                          >  7:277223/1‑135 (MQ=255)                   gCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGGCGGCGCCGGCGCGGGCCCggcgg                                                                                                                >  3:5683/1‑131 (MQ=255)                    cACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCACCCGCGCCGGCCCGTCttc                                                                                                               >  5:178411/1‑133 (MQ=255)                           gCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGACGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTg                                                                                                         <  4:27625/134‑1 (MQ=255)                                 aaCTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATACGGCGCGGGCCCGTCTGCGGCCTGCAccagg                                                                                                  >  4:249096/1‑135 (MQ=255)                                   cTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCCCGTCTGCGGCCTGCAccacgcc                                                                                                >  8:101429/1‑135 (MQ=255)                                           ggTACAAGGCAAGGGGGGCCAGCATCTGGTCCAGTTGCGCCTGGGACAACACCGGGTCCTTGGCTACGAGCGGGGCAGCGGACGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAg                                                                                        <  2:19323/135‑1 (MQ=255)                                                     aaTGGGGGCCTGCAGCTGGTCCAGTTGCTCCTGGGAGACCCCGGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGGGCGGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCa                                                                              <  8:320220/135‑1 (MQ=255)                                                        gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGACGGCGCCGGATCCGGCGCGGGCCCGGCTGCGGCCGGCCCCAGGCCAGCGGCCGCGCGCCGGCcaa                                                                            >  4:290339/1‑134 (MQ=255)                                                        gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCCGCGCGGGCACGTCTGCGGCCGCCCCCACGCCAGGGGCCGCGCGCAGGCCAAg                                                                           >  6:312343/1‑135 (MQ=255)                                                          gggCCAGCATCTGGTCCAGTTGCTCCAGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCGCGGGCCCGGCGGCGGCCGGCACCAGGCCCCGGGCAGCGCTCCGGCAGAGcc                                                                         >  5:292217/1‑134 (MQ=255)                                                               aGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGGCCGCTCCGCCGCGGGCCGGTCGGCGGCCTGCCCCCCCCCCCGGGCGCCGCGCAGGCACAGCAGTggc                                                                    >  7:49747/1‑133 (MQ=255)                                                               aGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCGCGGGCCCGGCTGCGGCCGGCCCCCCGCCCGGGGCCGCGCGCCCGCAAACCacagcc                                                                    >  5:198960/1‑130 (MQ=255)                                                               aGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGAACCGGCGCGGGCCCGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGcc                                                                    >  2:337447/1‑135 (MQ=255)                                                                gCATCTGGGCCTGTTGCGCCTGGGATAACCCCGGGTCCTTGGCTACGAGGGGGGCAGCGGAAGGCGCCGGACCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCg                                                                   <  7:40789/135‑1 (MQ=255)                                                                 cATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCCGCCCGGGCCCCGCGGCGGCCGGCCCCACGCCCCCGGCCGCGCGCCGGCCAACCACTGCCgc                                                                  >  6:342175/1‑134 (MQ=255)                                                                 cATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTACCGCGCGGGCCCGGCTGCGGCCGGCACCACGCCAGCGGCCGCGCGCCGGCCCACCAGTggcc                                                                   >  6:347939/1‑131 (MQ=255)                                                                        gccAGTTGCTCCTGGGAGAACACGGGGTCGTTGGCTCCGAGCGTGGTCGCGGAAGGCGCCGGACCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGaaa                                                           <  5:156478/134‑1 (MQ=255)                                                                         ccAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCGCCCGCGCGGGCCCCGCTGCGGCCGGCACCAGCCCCCCGGCAGCGCTCACGCCCAGCACAGCCCACAGCaaa                                                           >  7:57345/1‑134 (MQ=255)                                                                               gCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCCCCCGGCCAGGGGCAGCGCGCAGGCCAAGCAGTGCCGCCAGGAAAGCAAAc                                                     >  4:232817/1‑134 (MQ=255)                                                                                 tCCTGGGAGAACACGGGGTCCTTGGGCACGCGGGGGGCCGCGGACGGCGCCGGACCCGGCTCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAAt                                                  <  1:299358/135‑1 (MQ=255)                                                                                        agaACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCCGCTCCGGCGCGGGGCCGTCGGCGGCCGGCCCCACGCCCCGGGCCGCGCGCCGGCCAAGCAGTGCCGCCAGGAAAGCAAACAATGGCATgc                                           >  2:140376/1‑135 (MQ=255)                                                                                             acCGGGGCGTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGAGCGGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATc                                      <  4:5683/135‑1 (MQ=255)                                                                                               cGGGTCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGc                                    <  8:71610/135‑1 (MQ=255)                                                                                               cGGGGCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGc                                    <  2:50191/135‑1 (MQ=255)                                                                                                 gggccTTGGCTACGCGCGTGGCAGCGGACGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCtt                                  <  7:101429/132‑1 (MQ=255)                                                                                                 gctgcTTGGCTCCGCGCGTGGCAGCGGCAGGCGCGGGACCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCtt                                  <  7:43272/131‑1 (MQ=255)                                                                                                    ccTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCGCGGGCCCGTCTGCGGCCTGCCCCAGGCCCGGGGCAGCGCTCAGGCAAAGCAGTGCCGCCAGGAAAGCAAACAAAGGCCTGCGCATCGCGTATc                               >  2:144051/1‑135 (MQ=255)                                                                                                     ctgggCTCCGCGCGGGGCATCGGACGGCGCCGGAGCGGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTGTcc                              <  7:77586/132‑1 (MQ=255)                                                                                                      ttGGCTCCGCGGGTGGCTGGGGACGGCGCCGGAGGCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATccc                             <  5:312343/135‑1 (MQ=255)                                                                                                      ttGGCTACGAGCGGGGCAGCGTCAGGGGCGGCATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATccc                             <  6:133111/135‑1 (MQ=255)                                                                                                                   tgggcgCGTAAGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAc                <  1:286861/130‑1 (MQ=255)                                                                                                                    gccAGCGGACGGCGCGGGAGCCGGCGCGGGCACGGCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc               <  5:204377/133‑1 (MQ=255)                                                                                                                    gccAGCGGAAGGCGCCGGCTCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc               <  3:24761/133‑1 (MQ=255)                                                                                                                         cGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTg          <  6:172359/135‑1 (MQ=255)                                                                                                                          ggACGGCGCCGGAGCCGGGGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGCCAGACCGCGTgg         <  1:144051/135‑1 (MQ=255)                                                                                                                          ggAAGGCGCCGGAGCCGGCGCGGGCACGGCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTgg         <  8:186400/135‑1 (MQ=255)                                                                                                                          ggAAGGCGACGGATCCGGCGCGGGCACGGCTGCGGCCTGCACCAGGCCAGTTGCAGCGCTCAGGCAAAGAAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTgg         <  6:15183/135‑1 (MQ=255)                                                                                                                                 gcCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCCGGCCAGGGGCCGCGCTCAGGCAAAGCAGTGCCGCCAGGAAAGCAAACAATGGCAAGCGCATCGCTTTTCCCCATGAGGTCAGACCGCGTGGGTTTAgg  >  2:96759/1‑135 (MQ=255)                                                                                                                                       |                                                                                                                                 GGTAGGTTGACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGG  >  ACB122/3986488‑3986750

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