Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	318,084	G→C	19.0%	P107P (CCG→CCC)	oprD_3 →	Porin D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	318,084	0	G	C	19.0%	88.7 / 6.1	40	P107P (CCG→CCC)	oprD_3	Porin D
Reads supporting (aligned to +/- strand):  ref base G (22/10);  new base C (2/6);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 4.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCCAACTTCTCCTCCGGTTTCACCCAAGGCACCGTTGGCGTCGGTGTAGACGCATTCGGCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAACCACC‑AATGCCGAAGG‑CAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTG  >  ACB122/317954‑318213                                                                                                                                    |                                                                                                                                    tCGCCAACTTCTCCTCCGGTTTCACCCAAGGCACCGTTGGCGTCGGTGTAGACGCATTCGGCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCAc                                                                                                                                  <  1:224944/134‑1 (MQ=255) tCGCCAACTTCTCCTCCGGTTTCACCCAAGGCACCGTTGGCGTCGGTGTAGACGCATTCGGCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCAc                                                                                                                                  <  3:51842/134‑1 (MQ=255)  cGCCAACTTCTCCTCCGGTTTCACCCAAGGCACCGTTGGCGTCGGTGTAGACGCATTCGGCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGa                                                                                                                                <  2:408609/135‑1 (MQ=255)     cAACTTCTCCTCCGGTTTCACCCAAGGCACCGTTGGCGTCGGTGTAGACGCATTCGGCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAac                                                                                                                             <  1:63636/135‑1 (MQ=255)               tccGGTTTCACCCAAGGCACCGTTGGCGTCGGTGTTGACGCATTCGGCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACATATT‑AAACACCCAACAAACCACC‑aaaacc                                                                                                                  >  1:348661/1‑131 (MQ=255)               tccGGTTTCACCCAAGGCACCGTTGGCGTCGGTGTAGACGCATTCGGCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAACCACC‑AATg                                                                                                                    <  8:378936/133‑1 (MQ=255)                           cAAGGCACCGTTGGCGTCGGTGTAGACGCATTCGGCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAACCCCC‑AATGCCGAAGG‑CAACGcc                                                                                                     >  7:264807/1‑135 (MQ=255)                                      tGGCGTCGGTGTAGACGCATTCGGCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGGCATATTACACCCCCCCCGAAACCCCC‑AA‑GCCCCACG‑CCACGCCCGCAcgccgcc                                                                                          >  1:202736/1‑134 (MQ=255)                                      tGGCGTCGGTGTAGACGCATTCGGCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAACCACC‑AATGCCGAAGG‑CAACGCCAGCTCGCCGca                                                                                          >  4:344472/1‑135 (MQ=255)                                                            ggCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGCCTTCTT‑CAAGCCGCACGAAACCCCC‑AATGCCGAAGG‑CAACGCCAGCTCGCCGCCCAACCTGGCGCGTGGTGGtgct                                                                    >  2:126708/1‑135 (MQ=255)                                                            ggCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGCCTTCTT‑CAAGCCGCACGAAACCACC‑AATGCCGCAGG‑CAACGCCAGCTCGCCGCCCAACCTGGCGCggggggggggt                                                                    >  3:102856/1‑131 (MQ=255)                                                            ggCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGCCTTCTT‑CAAGCCGCACGAAACCACC‑AATGCCGAAGG‑CAACGCCCGCTCGCCGCCCCACCTGGCGCggggggggggg                                                                    >  5:118856/1‑131 (MQ=255)                                                            ggCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGCCGTCTT‑GAAGCCGCCCGCAACCCCC‑CATTCCGCAGG‑CAACGCCAGCGCGCCGCCCCACCGGGGGCgggggggtggg                                                                    >  1:315658/1‑133 (MQ=255)                                                            ggCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTGTT‑CAAGCCGCACGAAACCCCC‑AAAGCCCCAGG‑CAACCCCACCTCGCCGCACAACCCGGCGCgtgggggtggg                                                                    >  1:137680/1‑133 (MQ=255)                                                            ggCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAACCCCC‑AATGCCGAAGG‑CAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGtgct                                                                    >  8:261586/1‑135 (MQ=255)                                                            ggCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAACCACC‑AATGCCGAAGG‑CCACGCCAGCTCGCCACACAACCTGGCGCggggtggtggt                                                                    >  8:308255/1‑133 (MQ=255)                                                            ggCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAACCACC‑AATGCCGAAGG‑CAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGtgct                                                                    >  7:122817/1‑135 (MQ=255)                                                              ccTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAACCACC‑AATGCCGAAGG‑CAATGCCAGCTCGCCGCCCAACCGGGCGCGTGGGGGtgctgc                                                                  >  4:287116/1‑135 (MQ=255)                                                                    cGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCCCGAAACCACC‑AATGCCGAAGG‑CAACGCCAGCTCGCCGCCCAACCTGGCGCGGGGGGGTGCTGCTGTGaa                                                            >  7:135551/1‑135 (MQ=255)                                                                     gCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGCCTTCTT‑CAAGCCGCACGAAACCACC‑AATGCCGAAGG‑CAACGCCAGCTCGCCGCACACCCTGGCGCGTGGTGGTGCGGCTGGGAAg                                                           >  7:10864/1‑135 (MQ=255)                                                                     gCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAACCACC‑AATGCCGAAGG‑CAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGCAg                                                           >  7:325311/1‑135 (MQ=255)                                                                     gCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAACCACC‑AATGCCGAAGG‑CAACGCCACCACGCCGCACAACCTGGCGCGTGGTGGTGCTGGTGTGAAg                                                           >  4:103359/1‑135 (MQ=255)                                                                          gcgcCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAACCACC‑AATGCCGAAGG‑CCACGCCAGCTCGCCGCACAACCTGGCGCGTGGGGGTGCTGGTGTGAAGGTCCg                                                      >  6:245902/1‑135 (MQ=255)                                                                          gcgcCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAACCACC‑AATGCCGAAGG‑CAACGCCAGCTCGCCGGACAACCTGGCGCGGGGTGGTGCTGCTGTGAACTTCCg                                                      >  7:270022/1‑135 (MQ=255)                                                                          gcgcCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCCACTTCTT‑CAAGCCGCACGAAACCACC‑AATGCCGAAGG‑CAACGCCAGCTCGCCGCACAAACTGGCGCGTGGTGGTGCGGCTGTGAAGGTCCg                                                      >  4:115145/1‑135 (MQ=255)                                                                                       ggaaacggCCACACCGGGGGTGGGGGCGTCGACTTCTA‑CACCCGGCACGAACCCACA‑ATGGCCGACGC‑CAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACAcc                                         <  7:265956/129‑1 (MQ=255)                                                                                         cAAAGGCCACAACGGTGGTGGCGGCGTCGCCTTGTT‑CAAGCCGCACGAAACCACC‑AATGCCGCAGG‑CAACGCCAGCACGCCGCCCAACCTGGCGCGTGGTGGTGGTGCGGGGGAGGTCCGCATTTGCAACACCgg                                       >  6:148015/1‑134 (MQ=255)                                                                                          aaaGCCCACATGTGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAACCACG‑AATGCCGAAGG‑CAATGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTg                                      <  3:287116/135‑1 (MQ=255)                                                                                           aaGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAACCACC‑AATGCCGAAGG‑CAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACAccct                                       >  3:226882/1‑131 (MQ=255)                                                                                                           tggCGGCGTCGACTTATT‑AAAGCCGCACGAAACCACC‑AATGCCGAAGG‑CAACGCCAGCTAGCCGCACAACCTGGCGCGTGGTGGTGCTGGTGTGAAGTTACGCATTTCCAAAACCGTGCTTAATTACGGTgc                        >  4:103902/1‑131 (MQ=255)                                                                                                            gggggcggcGACTTATA‑CACCCCCCACCAACCCACA‑AAGCCCCACGC‑CAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAg                    <  6:118856/132‑1 (MQ=255)                                                                                                             gcggcggcgACTTC‑T‑CAACCCCCCCCTAACCCCCAAATCCCGACGC‑CAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGa                   <  3:344472/135‑1 (MQ=255)                                                                                                             gcggcgTCGACGTCTT‑CAAGCCGCACGAAACCACC‑AATTCCGCAGG‑CAACCCCAGCTCGCCGCACAACCTGGCGCGGGGTGGTGGGGGTGTGGAGGTGGGGATTTCCAACACCGGGGGGAAGTACgggggccc                     >  4:25289/1‑127 (MQ=255)                                                                                                                gcggcgACTTCTT‑CACGCCGCACTAAACCACC‑AATGCCGAAGG‑CAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGc                <  8:146285/135‑1 (MQ=255)                                                                                                                 cggctaCTGCTT‑CACGCCGAACGAACCCACA‑AAGCCCCACGC‑CCACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGcc               <  7:245286/130‑1 (MQ=255)                                                                                                                 cggccaCTTCTA‑CAAGCCCAACTAAA‑CACC‑AAATCCGAAGGCCAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGcc               <  7:113934/130‑1 (MQ=255)                                                                                                                 cgTCGACTTCTC‑AACCCCCCACTAACCCACA‑AATCCCCACGC‑CAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGcc               <  4:380620/135‑1 (MQ=255)                                                                                                                      aCTGCTA‑CACGCCCCACGAACCCCCC‑AAGGCCGCCGC‑CAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCAc          <  4:33866/135‑1 (MQ=255)                                                                                                                            tacacgCCGCAACAACCCACC‑AATCCCGACGC‑ACACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCa     <  7:305215/130‑1 (MQ=255)                                                                                                                               atgCCCCACGAACCCCCC‑ATTGCGGACGC‑CAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCGGCTGTGAAGGTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGtg  <  4:52157/133‑1 (MQ=255)                                                                                                                                    |                                                                                                                                    TCGCCAACTTCTCCTCCGGTTTCACCCAAGGCACCGTTGGCGTCGGTGTAGACGCATTCGGCCTGTACGCCGTGCGCCTGGATGGCGGCAAAGGCCACAACGGTGGTGGCGGCGTCGACTTCTT‑CAAGCCGCACGAAACCACC‑AATGCCGAAGG‑CAACGCCAGCTCGCCGCACAACCTGGCGCGTGGTGGTGCTGCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTG  >  ACB122/317954‑318213

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