Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,973,043	A→G	14.8%	N201S (AAC→AGC)	rpsA →	30S ribosomal protein S1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	1,973,043	0	A	G	14.8%	62.8 / 5.7	27	N201S (AAC→AGC)	rpsA	30S ribosomal protein S1
Reads supporting (aligned to +/- strand):  ref base A (11/12);  new base G (2/2);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-01

GACCAGAAGCGCAACAACGTTGTCGTTTCCCGTCGCAGCGTGCTGGAAGCCGAGAACAGCGCCGAGCGCGAAGCCCTGCTGGAAACCCTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCACATCACCGACATGGCCTGGAAGCGCATCAAGCACCCATCGGAAATCGT  >  ACB122/1972919‑1973145                                                                                                                             |                                                                                                       gACCAGAAGCGCAACAACGTTGTCGTTTCCCGTCGCAGCGTGCTGGAAGCCGAGAACAGCGCCGAGCGCGAAGCCCTGCTGGAAACCCTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGa                                                                                               <  6:158972/134‑1 (MQ=255) gACCAGAAGCGCAACAACGTTGTCGTTTCCCGTCGCAGCGTGCTGGAAGCCGAGAACAGCGCCGAGCGCGAAGCCCTGCTGGAAACCCTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGa                                                                                               <  8:82838/134‑1 (MQ=255)    cAGAAGCGCAACAACGTTGTCGTTTCCCGTCGCAGCGTGCTGGAAGCCGAGAACAGCGCCGAGCGCGAAGCCCTGCTGGAAACCCTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGACTACg                                                                                          <  4:201101/136‑1 (MQ=255)         gcgcAACAACGTTGTCGTTTCCCGTCGCAGCGTGCTGGAAGCCGAGAACAGCGCCGAGCGCGAAGCCCTGCTGGAAACCCTCCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCCCCGACTACGgcgc                                                                                      >  3:119862/1‑135 (MQ=255)          cgaaacaacGTTGTCGTTTCCCGTCGCAGCGTGCTGGATGCTGAGAACAGCGCCGAGCGCGAAGCCCTGCTGTAAACCCTGCAGGAAGGCCTGGAAGTCAAAGGTATCGTCAAGAACCTCACGGACTACGGCGCAt                                                                                    <  1:98889/133‑1 (MQ=255)                       tCGTTTCCCGTCGCAGCGTGCTGGAAGCCGAGAACAGCGCCGAGCGCGAAGCCCTGCTGGAAACCCTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGACTACGGCGCTTTCGTGGACCTggg                                                                       >  6:4656/1‑136 (MQ=255)                       tCGTTTCCCGTCGCAGCGTGCTGGAAGCCGAGAACAGCGCCGAGCGCGAAGCCCTGCTGGAAACCCTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGACTACGGCGCATTCGTGGACCTgg                                                                        <  8:195936/135‑1 (MQ=255)                          tttCCCGTCGCAGCGTGCTGGAAGCCGAGAACAGCGCCGAGCGCGAAGCCCTGCTGGACACCCTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGACTACGGCGCATTCGTGGACCTGGGCgg                                                                    >  4:135913/1‑136 (MQ=255)                               cGTCGCAGCGTGCTGGAAGCCGAGAACAGCGCCGAGCGCGAAGCCCTGCTGGAAACCCTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCg                                                               >  2:96523/1‑136 (MQ=255)                                         tgatGGAAGCCGAGAACAGCGCCGAGCGCGAAGCCCTGCTTGAAACCCTGCAGGATGGCCAGCAAGTCAATGGTATCGTCAAGAACCTCACCGACTACGGCGCTTTCGTGGACCTGGGCGGTATCGACGGCctgct                                                     <  4:119862/133‑1 (MQ=255)                                         tGCTGGAAGCCGAGAACAGCGCCGAGCGCGAAGCCCTGCTGGAAACCCTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAGCCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCctgct                                                     <  3:45489/136‑1 (MQ=255)                                                    gagaACATCGCCGAGCGCGATGCCCTGCTGGAAACCCTGCTGGAAGGCCAGCAAGTCAAAGGGATCGTCAAGAACCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCACATCACCg                                          <  4:145409/136‑1 (MQ=255)                                                               cGAGCGCGAAGCCCTGCTGGAAACCCTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCGCCTCACCGACATGGCCTgg                               >  5:193187/1‑136 (MQ=255)                                                                    gcgAAGCCCTGCTGGAAACCCTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAGCCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGc                                                          <  8:181222/104‑1 (MQ=255)                                                                    gcgAAGCCCTGCTGGAAACCCTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAGCCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGc                                                          >  7:181222/1‑104 (MQ=255)                                                                           ccTGCTGGAAACCCTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCACATCACCGACATGGCCTGGAAGCGCATCAAg                   <  3:184959/136‑1 (MQ=255)                                                                                   aaaCCCTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCACATCACCGACATGGCCTGGAAGCGCATCAAGCACCCATc           <  7:52849/136‑1 (MQ=255)                                                                                    aaCCCTGCAGGAAGGCCAGCAAGTCTAAGGTAGCGTCAAGAACCTCACAGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCACATCACCGACATGGCCTGGAAGCGCATCAAGCACCCATCg          <  2:154574/136‑1 (MQ=255)                                                                                       ccTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCACATCACCGACATGGCCTGGAAGCGCATCAAGCACCCATCGGa        >  4:38585/1‑135 (MQ=255)                                                                                          gCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCACATCACCGACATGGCCTGGAAGCGCATCAAGCACCCATCGGAAATc    >  8:198151/1‑136 (MQ=255)                                                                                            aGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAGCCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCACATCACCGCCATGGCCTGGAAGCGCATCAAGCACCCATCGGAAATCGt  >  5:153397/1‑136 (MQ=255)                                                                                            aGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCACATCACCGACATGGCCTGGAAGCGCATCAAGCACCCATCGGAAATCGt  >  5:22326/1‑136 (MQ=255)                                                                                            aGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCACATCACCGACATGGCCTGGAAGCGCATCAAGCACCCATCGGAAATCGt  >  5:123114/1‑136 (MQ=255)                                                                                                    cAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCACATCACCGACATGGCCTGGaa                             <  4:9597/101‑1 (MQ=255)                                                                                                    cAGCAAGTCAAAGCTATCGTCAAGAACCTCACCGACTAAGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCACATCACCGACATGGCCTGGaa                             >  3:9597/1‑101 (MQ=255)                                                                                                                        cAAGAACCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCACATCACCGACATGGCCTGGAAGCGCATCAAGCACCCATCg          <  4:171483/100‑1 (MQ=255)                                                                                                                        cAAGAACCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCACATCACCGACATGGCCTGGAAGCGCATCAAGCACCCATCg          >  3:171483/1‑100 (MQ=255)                                                                                                                             |                                                                                                       GACCAGAAGCGCAACAACGTTGTCGTTTCCCGTCGCAGCGTGCTGGAAGCCGAGAACAGCGCCGAGCGCGAAGCCCTGCTGGAAACCCTGCAGGAAGGCCAGCAAGTCAAAGGTATCGTCAAGAACCTCACCGACTACGGCGCATTCGTGGACCTGGGCGGTATCGACGGCCTGCTGCACATCACCGACATGGCCTGGAAGCGCATCAAGCACCCATCGGAAATCGT  >  ACB122/1972919‑1973145

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