Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	4,214,479	A→G	13.0%	K384R (AAG→AGG)	CJ019_03954 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	4,214,479	0	A	G	13.0%	81.8 / 4.2	32	K384R (AAG→AGG)	CJ019_03954	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (15/12);  new base G (2/2);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CGTTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTAGACATCCCCGCACCTG  >  ACB122/4214348‑4214613                                                                                                                                    |                                                                                                                                       cgtTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGTGGATg                                                                                                                                    >  1:86239/1‑136 (MQ=255) cgtTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATAGCCGCAGCCAGTACCCTGGGAGGATg                                                                                                                                    >  1:88914/1‑136 (MQ=255)          cGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATATACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATgatc                                                                                                                           <  1:201056/136‑1 (MQ=255)           ggCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTTGAGTGGTAGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTTCCCTGGGAAGATGGCAATgatcg                                                                                                                          <  3:47286/136‑1 (MQ=255)                agagCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGc                                                                                                                     <  6:132130/136‑1 (MQ=255)                 gagCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGGCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCa                                                                                                                    <  8:68333/136‑1 (MQ=255)                    cgcATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAAct                                                                                                                 <  7:37134/136‑1 (MQ=255)                     gcATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAActg                                                                                                                <  6:102646/136‑1 (MQ=255)                               ccGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAActggggctgggg                                                                                                       >  6:151606/1‑135 (MQ=255)                                      ggTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACATGCTGATCGCCGCAGCAAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTTGGGAAGCCCgt                                                                                               >  1:157672/1‑135 (MQ=255)                                                  gCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCCGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTc                                                                                   >  5:161127/1‑136 (MQ=255)                                                  gCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTc                                                                                   >  6:113913/1‑136 (MQ=255)                                                                   tGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGcc                                                                   >  7:61638/1‑135 (MQ=255)                                                                   tGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGcc                                                                   >  6:135097/1‑135 (MQ=255)                                                                           tGGTCGTTAGCACGGTAGTCAATCTACCAGCTGATCGCCGAAGCCAGTACCCTGGGAAGATGGCAATGATCGATAGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTc                                                          >  3:104183/1‑136 (MQ=255)                                                                             gTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCcg                                                        >  5:72501/1‑136 (MQ=255)                                                                             gTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCCCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCcg                                                        >  2:160944/1‑136 (MQ=255)                                                                               cGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCcggc                                                      <  8:61638/136‑1 (MQ=255)                                                                                gTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCcggcg                                                     >  7:127561/1‑136 (MQ=255)                                                                                      aCGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAggtg                                               <  4:153104/136‑1 (MQ=255)                                                                                             tCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGAtt                                                                             <  7:103915/99‑1 (MQ=255)                                                                                             tCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGAtt                                                                             >  8:103915/1‑99 (MQ=255)                                                                                                       tgctgATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCAc                              >  5:107847/1‑136 (MQ=255)                                                                                                         ctgATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAggtg                                               <  8:84758/117‑1 (MQ=255)                                                                                                         ctgATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAggtg                                               >  7:84758/1‑117 (MQ=255)                                                                                                         ctgATCGCCGCAGCACGTACCCTGGGATGATGGCAATGAGCCATCGCAACTGGGGCTGGGGAAGCCCGGTTGATTAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACg                             <  5:113913/135‑1 (MQ=255)                                                                                                            aTCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTca                          >  7:96915/1‑135 (MQ=255)                                                                                                                ccgcAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAAtt                     <  8:96915/136‑1 (MQ=255)                                                                                                                   cAGCCTGTACCGTGGGTAGATGGCAATGATCGATCGCAACTGGGTCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTa                  <  5:151606/136‑1 (MQ=255)                                                                                                                              cTGGGAGGATGGCAATTATCGATCGCAACTGGGGCGGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTAGACATCCCCGc       <  2:86239/136‑1 (MQ=255)                                                                                                                              cTGGGAGGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTAGACATCCCCGc       <  2:88914/136‑1 (MQ=255)                                                                                                                                  gAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTAGACATCCCCGCCCCt   >  4:127298/1‑136 (MQ=255)                                                                                                                                   aaGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTAGACATCCCCGCACCTg  >  7:169608/1‑136 (MQ=255)                                                                                                                                    |                                                                                                                                       CGTTAAACACGGCCAAGAGCGCATTTATGACCGCAATGGTGAGCTTGAAGCTGTGTTGATGTGGAATGCTGGAGTGGTCGTTAGCACGGTAGTCAATCTACCTGCTGATCGCCGCAGCCAGTACCCTGGGAAGATGGCAATGATCGATCGCAACTGGGGCTGGGGAAGCCCGGTTGATGAGATTCCAGATTTCGGAGAGCCTAGCCGCTCCGGCGGAGGTGGTTACGCTATCGACCACGTCACAATTGTAGACATCCCCGCACCTG  >  ACB122/4214348‑4214613

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