Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	1,113,257	C→T	7.6%	L18L (CTG→CTA)	pepA ←	Cytosol aminopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	1,113,257	0	C	T	7.6%	170.7 / 2.4	52	L18L (CTG→CTA)	pepA	Cytosol aminopeptidase
Reads supporting (aligned to +/- strand):  ref base C (26/22);  new base T (2/2);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGGTCTGGCCCGGCTTACCGGCCAGGTCGCCGCGCTTGAGCACGGCGCTGATGGCGCCTTCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAGCGCTCTGGACGGGCCTGGCCGCGTCGTCGCGC  >  ACB143/1113131‑1113388                                                                                                                               |                                                                                                                                    cagGGTCTGGCCCGGCTTACCGGCCAGGTCGCCGCGCTTGAGCACGGCGCTGATGGCGCCTTCGCAGGCCTGGTCGACAGCCTGGGCACCAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCg                                                                                                                            >  2:129586/1‑136 (MQ=255)      tCTGGCCCGGCTTACCGGCCAGGTCGCCGCGCTTGAGCACGGCGCTGATGGCGCCTTCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGtttt                                                                                                                       >  6:66520/1‑136 (MQ=255)           cccGGCTTACCGGCCAGGTCGCCGCGCTTGAGCACGGCGCTGATGGCGCCTTCTCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGTAAGGACCAGTGTGTCGGTTTTTACgg                                                                                                                  >  7:276332/1‑136 (MQ=255)                  tACCGGCCAGGTCGCCGCGCTTGAGCACGGCGCTGATGGCGCCTTCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATgcag                                                                                                            >  1:173433/1‑135 (MQ=255)                  tACCGGCCAGGTCGCCGCGCTTGAGCACGGCGCTCATGGCGCCTTCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATgcat                                                                                                            >  7:101283/1‑134 (MQ=255)                    ccGGCCAGGTCGCCGCGCTTGAGCACTGCGCTGATGGCGCCTTCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATgcagca                                                                                                          <  5:18788/135‑1 (MQ=255)                       gCCAGGTCGCCGCGCTTGAGCACGGCGCTGATGGCGCCTTCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTa                                                                                                      >  8:40183/1‑136 (MQ=255)                            gTCGCCGCGCTTGAGCACGGCGCTGATGGCGCCTTCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCtt                                                                                                 <  6:170736/136‑1 (MQ=255)                                        gAGCAAGGCGCTGATGGCGCCTTCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCTGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAg                                                                                      <  2:303223/135‑1 (MQ=255)                                                 gcTGATGGCGCCTTCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGt                                                                              >  7:201020/1‑134 (MQ=255)                                                   tGATGGCGCCTTCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAg                                                                           >  6:171615/1‑135 (MQ=255)                                                         cgcCTTCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGGTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATcccc                                                                    <  8:212136/136‑1 (MQ=255)                                                            cTTCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGaa                                                                 <  1:122521/136‑1 (MQ=255)                                                             ttCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCtt                                                                                                 >  6:114215/1‑103 (MQ=255)                                                             ttCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCtt                                                                                                 <  5:114215/103‑1 (MQ=255)                                                             ttCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAAt                                                                >  8:202073/1‑136 (MQ=255)                                                             ttCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACATACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAAt                                                                >  4:220951/1‑136 (MQ=255)                                                               cGCAGGCCTGGTCGACAGCCTGGGCAAGAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAAt                                                                <  6:95132/134‑1 (MQ=255)                                                               cGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAAt                                                                <  3:227419/134‑1 (MQ=255)                                                                  aGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATc                                                            >  1:60624/1‑135 (MQ=255)                                                                  aGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATc                                                            >  2:207581/1‑135 (MQ=255)                                                                     ccTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATa                                                        <  4:307186/136‑1 (MQ=255)                                                                      cTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATAc                                                       <  2:12019/136‑1 (MQ=255)                                                                       tGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATAc                                                       >  8:9618/1‑135 (MQ=255)                                                                         gTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGaa                                                    >  1:178358/1‑136 (MQ=255)                                                                            gACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATgg                                                 >  3:25690/1‑136 (MQ=255)                                                                            gACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATg                                                  <  1:51033/135‑1 (MQ=255)                                                                            gACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATg                                                  <  2:203690/135‑1 (MQ=255)                                                                            gACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATg                                                  <  5:33768/135‑1 (MQ=255)                                                                            gACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATg                                                  <  3:212190/135‑1 (MQ=255)                                                                                gCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCt                                             <  6:102411/136‑1 (MQ=255)                                                                                 ccTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTGCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCtg                                            >  8:34835/1‑136 (MQ=255)                                                                                       cAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGtt                                      <  6:182680/136‑1 (MQ=255)                                                                                       cAACAGTGCCGAGCTTGCGGTTTACACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGtt                                      <  6:305550/136‑1 (MQ=255)                                                                                       cAACAGAGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGtt                                      <  7:290496/136‑1 (MQ=255)                                                                                        aaCAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCTGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGtt                                      <  6:269994/135‑1 (MQ=255)                                                                                              gCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAgcgc                               >  2:124297/1‑136 (MQ=255)                                                                                                cGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAGCGctc                             >  2:186830/1‑136 (MQ=255)                                                                                                 gAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAGCGctct                            >  1:2394/1‑136 (MQ=255)                                                                                                  aGCTTGCGGTTTTCACCTACCGGAAGGACTAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAGCGctct                            <  6:241244/135‑1 (MQ=255)                                                                                                  aGCTTGCGGTTTTCACCTACCGGAAGGACTAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAGCGctct                            <  6:236703/135‑1 (MQ=255)                                                                                                  aGCTTGCGGTTTTCACCTACCGGAAGGACTAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAGCGctct                            >  5:241244/1‑135 (MQ=255)                                                                                                  aGCTTGCGGTTTTCACCTACCGGAAGGACTAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAGCGctct                            >  5:236703/1‑135 (MQ=255)                                                                                                  aGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAGCGCTCTg                           >  3:40351/1‑136 (MQ=255)                                                                                                       gCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAGCGCTCTGGACgg                      >  8:7423/1‑136 (MQ=255)                                                                                                          gTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGCTGCAGCAGCTACGCTTTTTACAACCAGTTACATGTCCGGATCACCGAATGATCGATACGAATGGCGCTGGGTGTTCCAGCGCTCTGGACGGGc                    >  7:111215/1‑135 (MQ=255)                                                                                                          gTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAGCGCTCTGGACGGGcc                   >  3:312418/1‑136 (MQ=255)                                                                                                                      cGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAAACAGTTCCATGTCAGGATACCCGAATGATAGATACTAATGGCGCTGTGTGTTCCATCGCTCTGGACGGGCCTGGCCGcgtcgt       >  2:276130/1‑136 (MQ=255)                                                                                                                       ggAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAGCGCTCTGGACGGGCCTGGCCGCgtcgtc      <  3:31509/136‑1 (MQ=255)                                                                                                                        gAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAGCGCTCTGGACGGGCCTGGCCGCGtcgtcg     >  3:282954/1‑136 (MQ=255)                                                                                                                          aGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAGCGCTCTGGACGGGCCTGGCCGCGTCGTcgcg   <  1:207581/136‑1 (MQ=255)                                                                                                                           ggACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAGCGCTCTGGACGGGCCTGGCCGCGTCGTCgcgc  <  4:40351/136‑1 (MQ=255)                                                                                                                               |                                                                                                                                    CAGGGTCTGGCCCGGCTTACCGGCCAGGTCGCCGCGCTTGAGCACGGCGCTGATGGCGCCTTCGCAGGCCTGGTCGACAGCCTGGGCAACAGTGCCGAGCTTGCGGTTTTCACCTACCGGAAGGACCAGTGTGGCGGTTTTTACGGATGCAGCAGCTACGCTTTTTACAACCAGTTCCATGTCAGGATCCCCGAATGATCGATACGAATGGCGCTGTGTGTTCCAGCGCTCTGGACGGGCCTGGCCGCGTCGTCGCGC  >  ACB143/1113131‑1113388

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