Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	3,598,388	T→G	28.1%	L184R (CTT→CGT)	CJ019_03412 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,598,388	0	T	G	28.1%	49.4 / 4.1	39	L184R (CTT→CGT)	CJ019_03412	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (4/23);  new base G (10/2);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 7.91e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCCAGCATTCGGCATCGGTCCATCCGCGCTCGACTTACCGCTTCAACGATGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATGGGCTTGCCCACCGGCCTTGCCTGGTCGCACGACGGCCATACGGTGTACTTCGTCGATTCGGCGGCGCGTGCGATCTACCGC  >  ACB122/3598263‑3598518                                                                                                                              |                                                                                                                                   gAGCCAGCATTCGGCATCGGTCCATCCGCGCTCGACTTACCGCTTCAACGATGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCGTTTTAGcgt                                                                                                                           >  7:302596/1‑133 (MQ=255)   gCCAGCATTCGGCATCGGTCCATCCGCGCTCGACTTACCGCTTCAACGATGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTg                                                                                                                         <  3:373314/135‑1 (MQ=255)                     tCCATCCGCGCTCGACTTACCGCTTCAACGATGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTGACTGGCAGTTGAGCGAcc                                                                                                       <  6:139687/135‑1 (MQ=255)                       cATCCGCGCTCGACTTACCGCTTCAACGATGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCGTTTTCGCGTTGACGGGCAGTTGAGCGCCCAg                                                                                                     >  7:27679/1‑135 (MQ=255)                        aTCCGCGCTCGACTTACCGCTGCAACGATGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTGACTGGCAGTTGAGCGACCAgg                                                                                                    <  4:194923/135‑1 (MQ=255)                              gcTCGACTTACCGCTTCAACGATGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTGACTGGCAGTTGAGCGACCAGGTGATTc                                                                                              <  7:59087/135‑1 (MQ=255)                                  gACTTACCGCTTCAACGATGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGGTTTGGCTGGCAGTTGAGCGGCCACGTGAttaccga                                                                                          >  7:293738/1‑130 (MQ=255)                                  gACTTACCGCTTCAACGATGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCGTTTTAGCGTTGTCAGGCAGTTGAGCGACCAGGGGATTAACGa                                                                                          >  4:277794/1‑135 (MQ=255)                                    cTTACCGCTTCAACGATGGCGCCTGCGCCCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGGTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATg                                                                                        <  3:318914/135‑1 (MQ=255)                                    cTNACCGCTTACACGATGGCGCCTGCGACCCACAGGGACGCGTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTGTTTCTCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATg                                                                                        <  6:238007/135‑1 (MQ=255)                                         cGCTTCAACGATGGCGCCTGCGACCCACAGGGACGCTGCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTGACTGGCAGTTTAGCGACCAGGTGATTCACGATGACATg                                                                                   <  6:167267/135‑1 (MQ=255)                                          gCTTCAACGCTGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCGTTTTCGCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATgg                                                                                  <  7:349092/135‑1 (MQ=255)                                          gCTGCACCGATGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATgg                                                                                  <  7:353615/135‑1 (MQ=255)                                           cTTCAACGAGGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACGGGGGCGCTTTTTCGCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATggg                                                                                 <  1:269169/135‑1 (MQ=255)                                                ccGATGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGAGGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATGGGCTTGc                                                                            <  1:128773/134‑1 (MQ=255)                                                 cGATGGCGCCTGCGCCCCACAGGGACGCTTACTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGTGGCGCGATTTCGCTTTGACTAGCACTTTACCGCACACGTCAATCAAGATCACATTGGcattcc                                                                           >  8:295407/1‑130 (MQ=255)                                                 cGATGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATGGGCTTGcc                                                                           <  3:67913/135‑1 (MQ=255)                                                 cGATGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCGTTTTGACTTTTACTGGCATTTTAGCGCCCAGGTGATTGACGATTCCATTGGCGTGcc                                                                           >  3:94892/1‑135 (MQ=255)                                                    tGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTAGCTTTTACTGGCAGTTTAGCGACCATGTGATTCACGATGACATGGGCTTTCCCAc                                                                        >  4:103138/1‑135 (MQ=255)                                                    tGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCGTTTTCGCTTTGTCTGGCAGTTGCGCGACCAGGTGATTCACCATGACATGGGCTTGGCCAc                                                                        >  4:97649/1‑135 (MQ=255)                                                                    ctgggACGCTTCGTCACCGGCCTGAGGGACGCGGCCCGCAGTGCCAATACCGGGGCGGTTTTTCGCTTTAACTGGCTGTTGAGCGACCAGGTGATTCACGATGACATGGGCTTGCCCACCGGCCTTGCCTGGTCg                                                        <  4:94892/133‑1 (MQ=255)                                                                    caGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATGGGCTTGCCCACCGGCCTTGCCTGGTCg                                                        <  6:63308/135‑1 (MQ=255)                                                                    caGGGACGCGTCGGCCCCGGCCTGCTGGACGAGGCCCTCAGGGACCATCCCGGGGCGCTGTTTCGCTTTGACTGGCAGTTGCGCGACCAGGTGATTCACGATGACATGGGCTTGCCCACCGGCCTTGCCTGGTCg                                                        <  2:57934/135‑1 (MQ=255)                                                                      gggACGCTTCGTCACCGGCCTGCTGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATGGGCTTGCCCACCGGCCTTGCCTGGTCGCa                                                      <  5:121339/135‑1 (MQ=255)                                                                      gggACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATGGGCTTGCCCACCGGCCTTGCCTGGTCGCa                                                      <  3:97649/135‑1 (MQ=255)                                                                      gggACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATGGGCTTGCCCACCGGCCTTGCCTGGTCGCa                                                      <  3:103138/135‑1 (MQ=255)                                                                           gCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATGGGCTTGCCCACCGGCCTTGCCTGGTCGCacgacg                                                 <  7:261790/135‑1 (MQ=255)                                                                             ttCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTTTCTGGCAGTTGAGCGCCCACGGGATTCACGATGTCATGGGCTTGCCCACCGGCCGTGCCTGGTCGCACGACGGc                                               >  5:197442/1‑135 (MQ=255)                                                                             ttCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCGTTTTCGCGTTTACGGGCAGGTGAGCGACCAGGGGGTTCGCGATGACATGGGCTTGCCCCCCCGCCATGCCGGGTCGCCCGCCGGc                                               >  5:56059/1‑135 (MQ=255)                                                                                gTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCGTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATGGGCTTGCCCACCGGCCTTGCCTGGTCGCACGACGGCCAt                                            <  3:36423/135‑1 (MQ=255)                                                                                 tCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTGACTGGCAGGTCAGCGACCAGGTGATTCACGAAGACATGGGCTTGCCCACCGGCCTTGCCTGGTCGCACGACGGCCAt                                            <  3:277794/134‑1 (MQ=255)                                                                                  cACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCGTTTTCGGGTTGACTGGGAGTTGAGCGCCCAGGGGATTGACGATGTCATGGGCTTGCCCCCCCGCCGTGCCGGGTCGCACGAAGGCCATAc                                          >  6:368417/1‑135 (MQ=255)                                                                                             tgtccGAGGCGGTCCGTGACATTACGGGGGGGCTTTTTCGGTTTGTGTTGATGTTTAGCGTCCTGGTGATTCACGATGACCTGGTCTTGCCCACCGGCCTTGCCTGGTCGCACGTCGGCCATACGGTGTACtt                                 <  5:254607/129‑1 (MQ=255)                                                                                                             tGACAATACCGGGGCGCGTTTTGGCGTTGACGGGCAGTTGAGCGCCCCGGTGATTAACGATGACATGGGCTTGCCCACCGGCCATGCCTGGTCGCCCGCCGGCCATACGGTGGACTTTGTCGATTCGGCGGcgcg               >  8:355659/1‑135 (MQ=255)                                                                                                              gACAATACCGGGGCGCGTTTTCGCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATGGGCTTGCCCACCGGCCTTGCCTGGTCGCACGACGGCCATACGGTGTACTTCGTCGATTCGGCGGCGCGt              <  3:231788/135‑1 (MQ=255)                                                                                                                        gggCGCTTTTTCGCTTGGACGGGCAGTTTAGCGACCAGGTGATTCACGATGACATGGGCTTGCCCACCGGCCTTGCCTGGTCGCACGACGGCCATACGGTGTACTTCGTCGATTCGGCGGCGCGTGCGATCTAcc    <  5:8137/135‑1 (MQ=255)                                                                                                                          gcgcTTTTTCGCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATGGGCTTGCCCACCGGCCTTGCCTGGTCGCACGACGGCCATACGGTGTACTTCGTCGATTCGGCGGCGCGTGCGATCTACcgc  <  1:292298/135‑1 (MQ=255)                                                                                                                          gcgcTTTTTAGCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATGGGCTTGCCCACCCGCCTTGCCTGGTCGCACGACGGCCATACCGTGTACATCGTCGATTAGGCGGCGCGTGGGATCTAcacc  >  6:263988/1‑132 (MQ=255)                                                                                                                          gcgcGTTTTCGCTTTGACTGGCAGTTGAGCGCCCAGGTGATTCACGATGACATGGGCTTGCCCACCGGCCGTGGCTGGTCGCACGCCGGCCATACGGTGTACTTCGTCGATTCGGCGGCGCGTGGGGTCTAcccc  >  7:49715/1‑133 (MQ=255)                                                                                                                              |                                                                                                                                   GAGCCAGCATTCGGCATCGGTCCATCCGCGCTCGACTTACCGCTTCAACGATGGCGCCTGCGACCCACAGGGACGCTTCGTCACCGGCCTGATGGACGAGGCCCTCAGTGACAATACCGGGGCGCTTTTTCGCTTTGACTGGCAGTTGAGCGACCAGGTGATTCACGATGACATGGGCTTGCCCACCGGCCTTGCCTGGTCGCACGACGGCCATACGGTGTACTTCGTCGATTCGGCGGCGCGTGCGATCTACCGC  >  ACB122/3598263‑3598518

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