Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	5,247,847	A→C	21.4%	L260R (CTG→CGG)	CJ019_04943 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	5,247,847	0	A	C	21.4%	60.5 / 3.9	46	L260R (CTG→CGG)	CJ019_04943	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/13);  new base C (9/2);  total (31/15)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGC  >  ACB122/5247729‑5247979                                                                                                                       |                                                                                                                                     cTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCt                                                                                                                      <  8:131303/135‑1 (MQ=255)             acTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCCGCGCTTGAAGTGGCTGGCGAGGTctgc                                                                                                          >  2:142490/1‑135 (MQ=255)                 cccTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGGTCTGCggtg                                                                                                      >  4:232204/1‑135 (MQ=255)                    tGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGCGTTCTGCTGTgggg                                                                                                   >  4:16322/1‑134 (MQ=255)                       gcgcCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCg                                                                                                  <  8:200529/133‑1 (MQ=255)                       gcgcCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCg                                                                                                  <  4:260707/133‑1 (MQ=255)                        cgcCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGGCt                                                                                               >  7:105033/1‑135 (MQ=255)                         gcCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGGTTGAAGTGGCTGGCGAGGTCTGCTGTGGGGGGCTg                                                                                              >  5:175095/1‑135 (MQ=255)                         gcCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCGTGCAGTGCCTGGGGAGTTGTGCTGTggggggc                                                                                                >  3:51652/1‑129 (MQ=255)                              ggcaggcTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCGGCGCGTGAAGTGCCTGGCGCGTTATTCTGTGGGCGACGGCTgg                                                                                          >  4:140091/1‑134 (MQ=255)                              ggcaggcTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGCAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATgg                                                                                          >  8:223629/1‑134 (MQ=255)                                          ggCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCGGCGCGGGCAGTGCCTGGCGCGTTGTGGGGTGGGCGGCGGCTGGGGAAACAGCCCCa                                                                             >  8:249610/1‑135 (MQ=255)                                            cTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCGTGAAGTGCCTGGCGAGTTGTGGGGTGGGGGACGGATGGGGAAACCGCACCAgg                                                                           >  7:2481/1‑135 (MQ=255)                                             tGCCTGCAAGGCTTTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCGGCGCTTGAAGTGCCGGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAggg                                                                          <  8:195616/135‑1 (MQ=255)                                             tGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAggg                                                                          <  4:79777/135‑1 (MQ=255)                                             tGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAggg                                                                          <  4:96140/135‑1 (MQ=255)                                                cagccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCGGCGCGTGCAGTGCCGGGCGAGTTCTGCGGTGGGCGCCTGCTGGGATAACAGCACCAGGGCGt                                                                       >  3:127401/1‑135 (MQ=255)                                                cagccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCGTGCAGTGCCGGGCGCGGTGTGCGGGGGGGGGCGGAGGGGGACACAGCCCCCGGGCgg                                                                       >  3:220766/1‑134 (MQ=255)                                                  gccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCGGCGCGTGCAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGCTGGGGAAACAGCCCCCGGGCGTcc                                                                     >  6:63607/1‑135 (MQ=255)                                                  gccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCGTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGTTGGGTAAACAGCACCAGGGCGTcc                                                                     >  2:15375/1‑135 (MQ=255)                                                       gCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCGGCGCTTGAAGTGCCTGGCGAGTTCTGGGGTGGGCGCCGGATGGGAAAACAGCACCCGGGCGTCCGGCtt                                                                >  8:225882/1‑135 (MQ=255)                                                         tGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCa                                                              <  2:199565/135‑1 (MQ=255)                                                         tGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCa                                                              <  4:264304/135‑1 (MQ=255)                                                         tGTCCTGTACCGCTTCCGCTGCTGCACCACCAAGTAACAACGGGCATTCCAGCAGCGACGCCCGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCa                                                              <  6:210000/135‑1 (MQ=255)                                                         tGTCCTGAACCGCTTCCGCTGCTTCACCCGCCATGAACAACGGCCATTCCAGCATCGACGCCAGGCGCTTGAAGTCCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCa                                                              <  5:91876/135‑1 (MQ=255)                                                              tGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGt                                                         >  8:295863/1‑135 (MQ=255)                                                                   cGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACCCCAGGCGCTTGAATTTGCAGGCGCGTTATGCTGTTGGCGTCCGAAGGGCAAACACCACCATGGCGTCCGGCATAAGCCggtagg                                                     >  2:115405/1‑129 (MQ=255)                                                                   cGCTTCCGCGGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGCCGCCAGGCGCTTGAAGTGCCGGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCa                                                    <  7:28338/135‑1 (MQ=255)                                                                        ccGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGCCTGATGGGCAAACAGCACCAGGGCGTCCGGCGTCAGCCGTTCGCCGAccc                                               >  8:372796/1‑134 (MQ=255)                                                                        ccGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCGGCGCGTGAAGTGCCTGGCGAGTTATGCTGTGGGCGACGGATGGGAAAACCGCCCCAGGGCGTGCGGCTTCCGCCGTTCGCAGCCCa                                               >  2:106970/1‑135 (MQ=255)                                                                        ccGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACGGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGGTCGCAGACCa                                               >  4:66860/1‑135 (MQ=255)                                                                        ccGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCGTGAAGTGCCTGGCGAGTTCTGCTGTGGGGGGCTGTTGGGGAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCa                                               >  6:271273/1‑135 (MQ=255)                                                                        ccGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCACGCGCCTGAAGTGCCGGGGGGCTTCTGCTGTGGGGGACTGAAGGGTTAACAGCACCACGGCGGCCGGAATCACCCCTTCGCCGCCCa                                               >  2:142689/1‑135 (MQ=255)                                                                          gctgctTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGGAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAgc                                             >  2:115216/1‑135 (MQ=255)                                                                          gctgctTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAgc                                             >  2:13800/1‑135 (MQ=255)                                                                          gctgctTCACCAGCCAGTAACAACGGGCATTCCACCATCGACGCCAGGCGCTTGAAGGGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAgc                                             <  1:115405/135‑1 (MQ=255)                                                                           ctgctTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCGTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGCCTGATGGGTAAACAGCACCAGGGCGTCCGGCGTCAGCCGTTCGCAGACCcgcg                                            >  2:279645/1‑135 (MQ=255)                                                                             gctTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCGTGAAGTGCCAGGCGAGTTCTGCTGTGGGCGCCGGATGGGAAAACAGCACCAGGGCGTCCCGCTTCCGCCGTTCGCCGACCAgcgcc                                          >  2:297491/1‑134 (MQ=255)                                                                                  aCCAGCCAGTAACAACGGGCATTCCAGCAGCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCt                                      <  7:372796/134‑1 (MQ=255)                                                                                           tAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTTCTGGTGGCGAAGGATGGGAAAACAGCACCAGGGCGTCCGGCTTAAGCCGTTCGCACACCATCGCAAGCACCagcagccg                            >  8:369605/1‑133 (MQ=255)                                                                                                   ggCATTCCAGCATCGACCCCCGGCGCGTGTAGTGCCTGGGGAGTTGTGGGGTGGGCGACTGATGGGGAAACCGCACCAGGGCGTCCCGCTTCAGCCGGTCGCAGACCAGCGCCAGCACCTTCAGCGGCTGcccct                    >  8:18595/1‑133 (MQ=255)                                                                                                    tcATTCCAGCAGCGACGCCAGGCGCGTGAAGTGCCTGGCGAGTTCTGCTGTGGGCTACTGATGGGAAAACAGCACCAGGGCGTCCGGCCTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGtt                   <  3:247921/134‑1 (MQ=255)                                                                                                      aTTCCAGCATCGACGCCAGGCGCGTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGCCGGATGGGAAAACCGCACCAGGGCGTCCGGCTTCAGCCGTTCGCCGACCAGCGCACGCTCCTGCAGCGGCTGCCCCTTCg                 >  6:97955/1‑135 (MQ=255)                                                                                                       ttCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGc                <  8:105033/135‑1 (MQ=255)                                                                                                                     ccAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTAAGCCGTTCGCAGACCAGCCCAAGCTCCTGCCTCGGCGGCCCCGTCGCCACCCGGCTCAcgc  >  2:309471/1‑135 (MQ=255)                                                                                                                     ccAGGCGCTTGAAGTGCCGGGCGAGTTCTGCTGTTGGCGACTGATGGGAAAACAGCCCCAGGGCGTCCGGCTTCAGCCGGTCGCAGACCAGCGCAAGCTCCTGCCGCGGCTGCCCGTTAGCCACCCGGATCCcgc  >  2:262049/1‑135 (MQ=255)                                                                                                                       |                                                                                                                                     CTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGC  >  ACB122/5247729‑5247979

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