Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	4,259,370:1	+A	15.6%	coding (85/1254 nt)	nicP_10 ←	Porin‑like protein NicP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	4,259,370	1	.	A	15.6%	91.7 / 20.1	32	coding (85/1254 nt)	nicP_10	Porin‑like protein NicP
Reads supporting (aligned to +/- strand):  ref base . (9/18);  new base A (3/2);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

CTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGA‑AATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTGGCATCTTCCAAGAAGCCCCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCGCCCCAGGCCA  >  ACB143/4259237‑4259502                                                                                                                                        |                                                                                                                                     cTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGA‑AATCCCGGTTGAAATAGTAGTAGCGCAAAGTCAGCCCCCCCTTGGCATAGTCCAAGCAGCCCCCctcat                                                                                                                                     >  8:217150/1‑134 (MQ=255)  tGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGA‑AATCGCGGGTGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTGGCATCATCCAAGAAGCCCCCCTC‑TTg                                                                                                                                   >  6:157932/1‑136 (MQ=255)    aaCTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGA‑AATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTGGCATCTTCCAAGAAGCCCCCCTC‑TTgc                                                                                                                                  <  2:242483/135‑1 (MQ=255)             aTGAAACCCTGGGCCCACTCTTCAACCTTGTCCTTGGCGGCCCCCGGGTCGCGGATAAACG‑GGGTGAACTAGTAGTTGCGGAATGTCAGCCCCCCCATTGCATATAACAACCAGCACCCCAC‑TTGCGCCCAGCCCg                                                                                                                        >  8:303915/1‑136 (MQ=255)                 aaCCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGA‑AATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTGGCATCTTCCAAGAAGCCCCCCTC‑TTGCGCCCAGGCCGGCaa                                                                                                                    <  7:371180/136‑1 (MQ=255)                 aaCCCTGCGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGA‑AATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTGGCATCTTCCAAGAAGCCCCCCTC‑TTGCGCCCAGGCCGGCaa                                                                                                                    <  2:185267/136‑1 (MQ=255)                                     aCCTTGCCCTTGGCGGCCCCCGGGTCGCGGA‑AATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCCCCCATGGCATCTTCCAAGAAGCCCCCCTC‑CTGCGCCCAGGCCGGCAAACCCAAGCCAAACGCCCGc                                                                                                 >  2:212027/1‑135 (MQ=255)                                           cccTGGGCGGCCCCCGGGTCGCGGA‑AATCGCGTTTGATGTAGTAGTTGCGCAATTTCAGCCCACCCTTGGCTTGTTCCAAGAAGCCCCCCTCATTGCGCCCAGGCCGGCAAACCAAGGCAAAACGCCAGCGAAAc                                                                                            <  4:18/135‑1 (MQ=255)                                            ccTTGGCGGCCCCCGGGTCGCGGA‑AATCGCGGTTGAAGTAGTAGTTGCGAAATGTCAGCCCAACCGTGGCATCTTCCAAAAAGCCCCCCTC‑TTTCGCCAACGCCGGCACAACAAGGCCAACCGCCAGCGCAACAGc                                                                                         >  7:25009/1‑136 (MQ=255)                                            ccTTGGCGGCCCCCGGGTCGCGGA‑AATCGCGGGTGAAGTAGTAGTTGCGCAATGTCAGCCCCCCCTTGGCATCTTCCAAGCAGCCCCCCTCATTGCGCCCACGCCCGCCAACCAAGGCCCACCGCCACCGAAACAg                                                                                          >  7:125699/1‑136 (MQ=255)                                             cTTGGGGGCCCCCGGGGCGCGGA‑AATCGCGGTTGCAGTAGTAGTTGCGCAATGTCAGCCCACCCTTGGCATCTTCCAAGAAG‑CCCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCa                                                                                       <  7:286199/136‑1 (MQ=255)                                              ttGGCGGCCCCCGGGTCGCGGG‑AATCGCGGGTGTAGTAGTAGTTGCGCAATGTGAGCCCCCCCTTTGGATCTTCCAAGAAG‑CCCCCGC‑ATGCGGCCAGGCCGGCCAACCAAGGCCCAACGCCCGCGAAACAGCCac                                                                                      >  2:35083/1‑135 (MQ=255)                                                  cggcCCCCGGGTCGCGGA‑AATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTGGCATCTTCCAAGAAGCCCCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGAc                                                                                   <  5:364603/136‑1 (MQ=255)                                                        ccGGGTCGCGGA‑AATCGCGGGTGTAGTAGTATTTGGGCAAAGTCAGCCCACCCGTGGCCTCTTACAAGAAG‑CCCCCAC‑ATGCGCCCCCGCCCGCAAACCAAGCCCCAAACCCCGCCAAACCACCACGAcccacccg                                                                            >  6:3712/1‑130 (MQ=255)                                                                           ggTTTAAGTAGTAGTTGCGCAATGTCCGCCCACCCTTGGCTTCTTCCAAGAAGCCCCCCTA‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGTCCGATCGGGCCAGCGTGCAAGTtgt                                                           <  5:215344/136‑1 (MQ=255)                                                                           ggTTGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTGGCATCTTCCAAGAAGCCCCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTtgt                                                           <  1:379712/136‑1 (MQ=255)                                                                             ttGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTGGCATCTTCCAAGAAG‑CCCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTc                                                        <  8:318609/136‑1 (MQ=255)                                                                             ttGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTGGCATCTTCCAAGAAG‑CCCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTc                                                        <  4:89188/136‑1 (MQ=255)                                                                                                 tGTCAGCCCACCCTTGGCATCTTCCAAGAAG‑CCCCCTC‑TTGCGCCCAGGCCGGCAAACCAATGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTcc                                    <  1:41062/136‑1 (MQ=255)                                                                                                    cAGCCCACCCTTGTCATCTTCAAAGAAGCCCAACTCATTGCGCCCAGGCCGGCAAACCAAGGCCAAACTCCAGAGAAACAGACAGGACCGAACGGGCCAGCGTGAAAGTTTTGTCAATCATTTTTATTGTtatc                                     >  6:98379/1‑134 (MQ=255)                                                                                                    cAGCCCACCCTTGGCATCTTCCAAGAAGCCCCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGCACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCAc                                  >  2:218325/1‑136 (MQ=255)                                                                                                         cACCCTTGGCATCTTCCAAGAAGCCCCCCTCATTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGtt                               >  3:53819/1‑135 (MQ=255)                                                                                                              ttGGCATCTTCCAAGAAGCCCCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAgg                        >  2:287453/1‑136 (MQ=255)                                                                                                                      ttCCAAGAAGCCCCCCTA‑TTGCGCCCAGGCCGGAAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCgcagca                 <  6:37442/135‑1 (MQ=255)                                                                                                                        ccAAGAAGCCGCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGGGCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCgcagca                 <  8:23547/133‑1 (MQ=255)                                                                                                                        ccAAGAAGCCCCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCgcagca                 <  2:41295/133‑1 (MQ=255)                                                                                                                        ccAAGAAG‑CCCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCgcagca                 <  1:35083/132‑1 (MQ=255)                                                                                                                         cAAGAAGCCCCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTc                                                        >  7:138570/1‑93 (MQ=255)                                                                                                                         cAAGAAGCCCCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTc                                                        <  8:138570/93‑1 (MQ=255)                                                                                                                          aagaagCCCCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCg            <  6:185568/136‑1 (MQ=255)                                                                                                                          aagaagCCCCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGCCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCg            <  8:31549/136‑1 (MQ=255)                                                                                                                           agaagCCCCCCTCATTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCg            <  4:49772/136‑1 (MQ=255)                                                                                                                                    ccTC‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCTAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCGCCCCAGGCCa  <  1:148063/136‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     CTGAACTTGAAGATGAAACCCTGGGCCCACTCTTCAACCTTGCCCTTGGCGGCCCCCGGGTCGCGGA‑AATCGCGGTTGAAGTAGTAGTTGCGCAATGTCAGCCCACCCTTGGCATCTTCCAAGAAGCCCCCCTC‑TTGCGCCCAGGCCGGCAAACCAAGGCCAAACGCCAGCGAAACAGCCAGGACCGAACGGGCCAGCGTGCAAGTTGTGTCCATCATTATTATTGTTCTCCACGTTCAATAGGCGCAGCAGACCGCCCCAGGCCA  >  ACB143/4259237‑4259502

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