Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	1,005,859	G→A	7.3%	K407K (AAG→AAA)	prfC →	Peptide chain release factor RF3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	1,005,859	0	G	A	7.3%	180.6 / 3.7	55	K407K (AAG→AAA)	prfC	Peptide chain release factor RF3
Reads supporting (aligned to +/- strand):  ref base G (15/36);  new base A (2/2);  total (17/38)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

GGCCTGCACAACCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGG  >  ACB143/1005725‑1005992                                                                                                                                       |                                                                                                                                      ggCCTGCACAACCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAg                                                                                                                                       <  7:66317/135‑1 (MQ=255)          aaCCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCg                                                                                                                               <  5:235850/134‑1 (MQ=255)          aaCCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGt                                                                                                                              >  6:30151/1‑135 (MQ=255)                     cATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTgcag                                                                                                                  <  3:259286/136‑1 (MQ=255)                     cATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTgcag                                                                                                                  >  3:198757/1‑136 (MQ=255)                               ggCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGcc                                                                                                         >  7:235965/1‑135 (MQ=255)                                gCGACCCCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGCAAGACCCGCTGAAATCCACGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCga                                                                                                       <  4:183258/136‑1 (MQ=255)                                gCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCga                                                                                                       <  5:185805/136‑1 (MQ=255)                                gCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCga                                                                                                       <  3:49104/136‑1 (MQ=255)                                   acacCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCgaaga                                                                                                    <  4:190198/136‑1 (MQ=255)                                   acacCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCgaaga                                                                                                    <  7:206571/136‑1 (MQ=255)                                   acacCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCgaaga                                                                                                    <  7:206572/136‑1 (MQ=255)                                    cacCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCgaag                                                                                                     >  2:10422/1‑134 (MQ=255)                                    cacCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCgaag                                                                                                     >  3:193996/1‑134 (MQ=255)                                    cacCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCgaag                                                                                                     >  3:55335/1‑134 (MQ=255)                                             cGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACt                                                                                          <  3:88304/136‑1 (MQ=255)                                             cGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACt                                                                                          <  3:92485/136‑1 (MQ=255)                                             cGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACt                                                                                          <  8:16716/136‑1 (MQ=255)                                               aaGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCa                                                                                        <  5:101049/136‑1 (MQ=255)                                                  gCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGgt                                                                                     <  6:138121/136‑1 (MQ=255)                                                    gAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGgtgt                                                                                   >  8:187350/1‑136 (MQ=255)                                                    gAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCCACTCAGgtgt                                                                                   >  8:140721/1‑136 (MQ=255)                                                               cTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTTAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAg                                                                        <  8:100429/136‑1 (MQ=255)                                                                     cGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCgcag                                                                   >  4:135222/1‑135 (MQ=255)                                                                        tATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAAc                                                               >  5:152252/1‑136 (MQ=255)                                                                          tCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGa                                                             >  1:108746/1‑136 (MQ=255)                                                                          tCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGa                                                             <  6:198645/136‑1 (MQ=255)                                                                              gCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGTcatc                                                         <  3:6751/136‑1 (MQ=255)                                                                              gCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGAcatc                                                         <  2:101489/136‑1 (MQ=255)                                                                                  ttCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATcc                                                     >  8:224993/1‑136 (MQ=255)                                                                                       gcCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATAATCCTGGGt                                                >  8:259494/1‑136 (MQ=255)                                                                                         cGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCTCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGc                                              <  4:158343/136‑1 (MQ=255)                                                                                         cGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGc                                              <  5:247946/136‑1 (MQ=255)                                                                                         cGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGc                                              <  1:207673/136‑1 (MQ=255)                                                                                         cGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGc                                              <  7:193442/136‑1 (MQ=255)                                                                                         cGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGc                                              <  1:181006/136‑1 (MQ=255)                                                                                           gAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGcgg                                            <  8:57532/136‑1 (MQ=255)                                                                                           gAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGcgg                                            <  6:103725/136‑1 (MQ=255)                                                                                           gAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGcgg                                            <  5:127600/136‑1 (MQ=255)                                                                                           gAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGcgg                                            <  7:84935/136‑1 (MQ=255)                                                                                           gAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGcgg                                            <  6:225097/136‑1 (MQ=255)                                                                                                   ccgccgCGTACGCCTGAAAGACCCGCTGAAATCCAAGCATCTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGgtgt                                     <  4:12046/135‑1 (MQ=255)                                                                                                         cgTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGc                                                                                                                                <  5:237448/38‑1 (MQ=255)                                                                                                         cgTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGc                                                                                                                                >  6:237448/1‑38 (MQ=255)                                                                                                         cgTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTc                                                                                                                             <  8:74434/41‑1 (MQ=255)                                                                                                         cgTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTc                                                                                                                             >  7:74434/1‑41 (MQ=255)                                                                                                         cgTACGCCTGAAAGACCCGCTGAAATCCAAACAACTGCGTc                                                                                                                             <  2:235309/41‑1 (MQ=255)                                                                                                         cgTACGCCTGAAAGACCCGCTGAAATCCAAACAACTGCGTc                                                                                                                             <  2:231082/41‑1 (MQ=255)                                                                                                         cgTACGCCTGAAAGACCCGCTGAAATCCAAACAACTGCGTc                                                                                                                             >  1:235309/1‑41 (MQ=255)                                                                                                         cgTACGCCTGAAAGACCCGCTGAAATCCAAACAACTGCGTc                                                                                                                             >  1:231082/1‑41 (MQ=255)                                                                                                                        cccGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGcc               <  5:283052/136‑1 (MQ=255)                                                                                                                               aaaTCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCt        <  7:19056/136‑1 (MQ=255)                                                                                                                                  tCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTga     <  6:152252/136‑1 (MQ=255)                                                                                                                                    cAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTgaag   <  3:269028/136‑1 (MQ=255)                                                                                                                                    cAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTgaag   <  1:16931/136‑1 (MQ=255)                                                                                                                                     aaGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTgaagg  <  2:237673/136‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      GGCCTGCACAACCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGG  >  ACB143/1005725‑1005992

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