Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,063,986	T→G	32.3%	D151A (GAC→GCC)	CJ019_01001 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	1,063,986	0	T	G	32.3%	32.6 / 2.4	40	D151A (GAC→GCC)	CJ019_01001	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/23);  new base G (5/8);  total (9/31)
Fisher's exact test for biased strand distribution p-value = 8.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCCAGGTCGGTGGCCACCTCGATGGCAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCTGCCCAGGTGGGCATTGGCGGACAGCGCATACACCAGC  >  ACB122/1063855‑1064120                                                                                                                                     |                                                                                                                                       ggCCCCGGTCGGGGGCCACCTCGATGGCAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCACTATCTGCAGCCGCGTACCTGACCAGAAGCGCCCGGCCCCGGACCTGGCG‑GCGTCCg                                                                                                                                     <  7:24547/135‑1 (MQ=255) ggCCCCGGGCGGTGGCCACCTCGAGGGCAGCTTCCAGAAGCACGTGGCCCGCCCGGCTCAGCCGGCCAAGGCCGACCAGTCTCTGCCGCCGCGTACCTGACCAGCCGCGCCCGGCCCCGGACCGGGCG‑GCGTCCg                                                                                                                                     <  6:155636/135‑1 (MQ=255) ggCCCAGGTCGGTGGCCACCTCGATGACAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGACGCGCCCGGCCCCGGACCTGGCG‑GCGGCCg                                                                                                                                     <  3:264625/135‑1 (MQ=255) ggCCCAGGGCGGTGGCCACCTGGATGGCAGCTGCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTCCCTGACCAGAAGCGCCCGGCCCCGGACCTGGCG‑GCGTCCg                                                                                                                                     <  3:131726/135‑1 (MQ=255) ggCCCAGGGCGGTGGCCACCTCGATGGCAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGACGCGCCCGGCCCCGGACCTGGCG‑GCGTCCg                                                                                                                                     <  8:82732/135‑1 (MQ=255) ggCCCAGGGCGGTGGCCACCTCGATGGCAGCTTCCAGAAGCACGTGGCCCGCCAGGCGCAGCAGGCCAAGCCCGACCAGTATCGGCAGCAGCGAACCTCACCAGCCGCGCCCGGCCCCGGACCTGGCG‑GCGGCCg                                                                                                                                     <  1:224759/135‑1 (MQ=255)       ggTCGGGGGCCACCTCGATGGCAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGACGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAg                                                                                                                               <  4:237510/135‑1 (MQ=255)           gttGGCCACCTCGATGGCAGCTTCAAGTAGCTAGGGCCCCGCCAGTCGCAGCAGGCCAAGCCCTACCCGAATCTGCAGCAGCGTCCGTCACCAGACGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTa                                                                                                                           <  6:282209/133‑1 (MQ=255)           ggTGGCCCCCGCGAGGGCCGCTGCCAGTACCCTGTGGCCCGCCAGGCTCAGAACGCCAGGGCCGACAAGTAGCGGCAGCACCGTACCTGACCAGACGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTa                                                                                                                           <  2:157704/135‑1 (MQ=255)           ggTGGCCACCTCGATGGCAGCTTCCAGTAGCACGGGGCCCGCCCGGCGCAGCAGGCCAAGGCCGACCAGTAGCTGCAGCAGCGTACCTGACCAGAAGCGCCCGGCCCCGGACCTGGCG‑GCGGCCGGGGCAGGGTa                                                                                                                           <  8:170585/135‑1 (MQ=255)             tGGCCACCTCGATGGCAGCTTCCATTCGAATGTGCCCCCCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCACAAGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAgg                                                                                                                         <  3:95217/135‑1 (MQ=255)             gggCCACCGCGATGGCAGCGCCCAGTAGCAAGGGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCATAAGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAgg                                                                                                                         <  2:317870/134‑1 (MQ=255)               gCCACCTCGATGGCAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGGCCCCCGGCCGGGCG‑GCGGGCGGGGGAggggtgggg                                                                                                                       >  8:209712/1‑129 (MQ=255)               gCCACCTCGATGGCAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGGCCCCCGGCCAGGCGCGCGCCCTGGGCgggtgtggg                                                                                                                        >  4:68306/1‑129 (MQ=255)               gCCACCTCGATGGCAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGGCCCCCGCCCGGGCG‑GCGGGCGGGGCCGGGGAgggg                                                                                                                       >  3:97779/1‑133 (MQ=255)               gCCACCTCGATGGCAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGGCCCCCGCCCAGGCC‑GCGGCCCGGggaggggagggg                                                                                                                       >  5:108326/1‑133 (MQ=255)                           gCAGCTGCCAGTCGCAAGGGGCCCCCCCGGCGCCCCCGGCCACGCCCGACCAGTATCTGCCGCAGCGAACCTGACCAGACGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATg                                                                                                            <  6:223076/134‑1 (MQ=255)                               cTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGAATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGcgcg                                                                                                       <  4:176018/135‑1 (MQ=255)                                   cAGTCGCACGTGGCCCGCCCGGCGCAGCAGCCCAAGGCCGACCAGTAGCTGCAGCAGCGTCCCTCACCACACGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGAc                                                                                                     <  3:164056/133‑1 (MQ=255)                                    aGTCGCACGTGGCCCGCCAGGCTCAGCAGCCCACGCCCGACCAGTATCTGCAGCCGCGTACCTGACCAGACGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGa                                                                                                  <  3:140156/135‑1 (MQ=255)                                       cccAAGGGGCCCGCCCGGCCCAGCACGCAAACCCCCACCAGAACCCCCCCCCGCGAACCTCCCCAGACGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGAt                                                                                                 <  6:228430/131‑1 (MQ=255)                                                cccGCCAGGCTCAGCAGGCCAAGGCCGACCTGTATCTCCCGCAGCGTCCCTGACCTCACGCGCCCGGCCCCGGACCTGGCG‑GCGGCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCAtt                                                                                      <  1:137945/135‑1 (MQ=255)                                                              cagGCCAAGGCCGACCAGTATCTGCCGCAGCGTACCTGACCAGAAGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGc                                                                        <  5:50369/135‑1 (MQ=255)                                                                ggCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCa                                                                      <  1:210876/135‑1 (MQ=255)                                                                  ccAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCACACGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCAGc                                                                    <  5:136992/135‑1 (MQ=255)                                                                  ccAAGGCCGACCAGAATCTCCAGCAGCGTACCTGACCAGACGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCAGc                                                                    <  6:349304/135‑1 (MQ=255)                                                                    aaGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGGCCCCCGGCCTGGCG‑GCGGCCCGGGGAGGGGTGGTGCGTGGGGAATGGCGCGCCAGCTGCCAGAGCAATTCCAGGTAACCCGCCAAccc                                                                  >  8:83626/1‑135 (MQ=255)                                                                              aGTATCTGCAGCAGCGTACCTGACCAGCAGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCAGCCCCAGCAATAGc                                                        <  6:75587/135‑1 (MQ=255)                                                                              aGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGGCCCCGGCCCGGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGGGTAATGGCGCGCCGGATGCCAGAGCCTTGCCAGGTAACCGGCCCGCCCCCGCAATa                                                          >  4:255467/1‑133 (MQ=255)                                                                                tATCTGCTGCTGCGTAGCTGACCAGACGCGCCCGGCCGCGGACCTGGCG‑GCGGCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCAGCCCCAGCAATAGCAg                                                      <  6:91889/135‑1 (MQ=255)                                                                                tATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCAGCCCCAGCAATAGCAg                                                      <  7:316776/135‑1 (MQ=255)                                                                                tATCTCCAGCAGCGTACCTGCCCAGAAGCGCCCGGCCCCGGACCTGGCG‑GCGGCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCAGCCCCAGCAATAGCAg                                                      <  7:319736/135‑1 (MQ=255)                                                                                          cgcGGACCTCCCCACACGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCa                                            <  3:68306/134‑1 (MQ=255)                                                                                          agcGTACCTGACCAGAAGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCa                                            <  4:103831/135‑1 (MQ=255)                                                                                               aCCTGACCAGAAGCGCCCGGCCCCGGGCCGGGCG‑GCGGGCGGGGGAGGGGTGGGGGGTGCGGTATGGGGCGCCTGATGGCAGAGCCTTTGCCGGTAACCCGCCCCCCCCCCCCAAAGCAGTGCCCCGTCACCCCt                                       >  2:232172/1‑135 (MQ=255)                                                                                                                 ggCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCTGCCCAGGTGGGCATTGGc                     >  8:317231/1‑135 (MQ=255)                                                                                                                   ccccGGACCTGGCG‑GCGGCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCTGCCCAGGTGGGCATTGGCgg                   <  1:157040/135‑1 (MQ=255)                                                                                                                       ggACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTACTGTCGCGACTGATGCCAGAGCATTGCCAGGGAAGCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCTGCCCAGGTGGGCATTGGCGGACAg               <  8:374764/135‑1 (MQ=255)                                                                                                                        gACCTGGCG‑GCGGCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCTGCCCAGGTGGGCATTGGCGGACAgc              <  6:124735/135‑1 (MQ=255)                                                                                                                                g‑gcgTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCTGCCCAGGTGGGCATTGGCGGACAGCGCATacac      >  6:91877/1‑135 (MQ=255)                                                                                                                                     tCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCTGCCCAGGTTGGCATTGTCTGACAGCGCATACACCAGc  <  5:136640/135‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       GGCCCAGGTCGGTGGCCACCTCGATGGCAGCTTCCAGTAGCAAGTGGCCCGCCAGGCTCAGCAGGCCAAGGCCGACCAGTATCTGCAGCAGCGTACCTGACCAGAAGCGCCCGGCCCCGGACCTGGCG‑GCGTCCGGGGCAGGGTAGGTGCGTGCGTAATGGCGCGACTGATGCCAGAGCATTGCCAGGTAACCGGCCAGCCCCAGCAATAGCAGTGCCCCTTCAACCCTGCCCAGGTGGGCATTGGCGGACAGCGCATACACCAGC  >  ACB122/1063855‑1064120

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