Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	1,486,333	A→T	17.4%	S497S (TCT→TCA)	U712_07630 ←	Ribonuclease J 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	1,486,333	0	A	T	17.4%	50.5 / 4.8	23	S497S (TCT→TCA)	U712_07630	Ribonuclease J 1
Reads supporting (aligned to +/- strand):  ref base A (5/13);  new base T (2/2);  total (8/15)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TTTTTTCGTACAAGAAAGGTGCAAGCGTGTCTGTAATTTCGTTTTTGATTTCAGACCATTGAGTCGTTTTTCGTTCCATTACTTTTTGTAAGTGATTTGAAATGAGCTCTTGAGCGTCGTTGATCAAGTCACCAGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAGTCGTCCATGTCAATGCTGACAACAACGATGACAAGTCCTTCTTCAGAGAGAATTCTGCGATCACGAAGTAC  >  CP006881/1486200‑1486464                                                                                                                                      |                                                                                                                                    ttttttCGTACAAGAAAGGTGCAAGCGTGTCTGTAATTTCGTTTTTGATTTCAGACCATTGAGTCGTTTTTCGTTCCATTACTTTTTGTAAGTGATTTGAAATGAGCTCTTGAGCGTCGTTGATCAAGTCACCAGAtt                                                                                                                                 <  1:182825/138‑1 (MQ=255) ttttttCGTACAAGAAAGGTGCAAGCGTGTCTGTAATTTCGTTTTTGATTTCAGACCATTGAGTCGTTTTTCGTTCCATTACTTTTTGTAAGTGATTTGAAATGAGCTCTTGAGCGTCGTTGATCAAGTCACCAGAtt                                                                                                                                 <  1:22089/138‑1 (MQ=255)           cAAGAAAGGTGCAAGCGTGTCTGTAATTTCGTTTTTGATTTCAGACCATTGAGTCGTTTTTCGTTCCATTACTTTTTGTAAGTGATTTGAAATGAGCTCTTGAGCGTCGTTGATCAAGTCACCAGATTCTCTCATGTa                                                                                                                       <  7:37834/138‑1 (MQ=255)                  ggTGCAAGCGTGTCTGTAATTTCGTTTTTGATTTCAGACCATTGAGTCGTTTTTCGTTCCATTACTTTTTGTAAGTGATTTGAAATGAGCTCTTGAGCGTCGTTGATCAAGTCACCAGATTCTCTCATGTACACAAATc                                                                                                               <  3:275397/139‑1 (MQ=255)                      cAAGCGTGTCTGTAATTTCGTTTTTGATTTCAGACCATTGAGTCGTGTTTGGGTACAATAATTTTTTTTATTGTTTTTTAATTATCTGTTGTGAGGCGTTGGTAAAATCACCCCATTATGTCATGTTAACAAAAcca                                                                                                             >  8:260373/1‑136 (MQ=255)                       aaGCGTGTCTGTAATTTCGTTTTTGATTTCAGACCATTGAGTCGTTTTTCGTTCCATTACTTTTTGTAAGTGATTTGAAATGAGCTCTTGAGCGTCGTTGATCAAGTCACCAGATTCTCTCATGTACACAAATCCTCTg                                                                                                          <  7:192198/139‑1 (MQ=255)                                      ttCGTTTTTGATTTCAGACCATTGAGTCGTTTTTCGTTCCATTACTTTTTGTAAGTGATTTGAAATGAGCTCTTGAGCGTCGTTGATCAAGTCACCAGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAgg                                                                                            <  4:169698/138‑1 (MQ=255)                                                              aGTCGTTTTTCGTTCCATTACTTTTTGTAAGTGATTTGAAATGAGCTCTTGAGCGTCGTTGATCAAGTCACCAGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAgtcgtc                                                                    <  4:27953/138‑1 (MQ=255)                                                                     tttCGTTCCATTACTTTTTGTAAGTGATTTGAAATGAGCTCTTGAGCGTCGTTGATCAAGTCACCAGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAGTCGTCCATCTCa                                                             <  8:230459/138‑1 (MQ=255)                                                                         gTTCCATTACTTTTTGTAAGTGATTTGAAATGAGCTCTTGAGCGTCGTTGATCAAGTCACCAGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAGTCGTCCATGTCAATGCt                                                        <  7:260373/139‑1 (MQ=255)                                                                                 aCTTTTTGTAAGTGATTTGAAATGAGCTCTTGAGCGTCGTTGATCAAGTCACCAGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAGTCGTCCATGTCAATGCTGacaaca                                                 <  1:188355/138‑1 (MQ=255)                                                                                       tGTAAGTGATTTGAAATGAGCTCTTGAGCGTCGTTGATCAAGTCACCTGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAGTCGTCCATGTCAATGCTGACAACAACGATGa                                          >  5:267363/1‑139 (MQ=255)                                                                                       tGTAAGTGATTTGAAATGAGCTCTTGAGCGTCGTTGATCAAGTCACCTGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAGTCGTCCATGTCAATGCTGACAACAACGATGa                                          >  5:271597/1‑139 (MQ=255)                                                                                         tAAGTGATTTGAAATGAGCTCTTGAGCGTCGTTGATCAAGTCACCAGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAGTCGTCCATGTCAATGCTGACAACAACGATGACa                                        <  7:152608/139‑1 (MQ=255)                                                                                                        gAGCTCTTGAGAGTCGTTGATCAAGTCACCAGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAGTCGTCCATGTCAATGCTGACAACAACGATGACAAGTCCTTCTTCaga                          <  6:243034/138‑1 (MQ=255)                                                                                                         aGCTCTTGAGCGTCGTTGATCAAGTCACCAGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAGTCGTCCATGTCAATGCTGACAACAACGATGACAAGTCCTTCTTCaga                          >  3:82981/1‑137 (MQ=255)                                                                                                              ttGAGCGTCGTTGATCAAGTCACCTGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAGTCGTCCATGTCAATGCTGACAACAACGATGACAAGTCCTTCTTCAGAGAGaa                     <  6:267363/137‑1 (MQ=255)                                                                                                              ttGAGCGTCGTTGATCAAGTCACCTGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAGTCGTCCATGTCAATGCTGACAACAACGATGACAAGTCCTTCTTCAGAGAGaa                     <  6:271597/137‑1 (MQ=255)                                                                                                                   cgtcgtTGATCAAGTCACCAGATTCTCTCCTGAACACAAATCCTCTGGAGATCAAATCAGGACCCGCTGAAATCTTTAAGTCGTCCATGTCAATGCTGCCAAAAACGATGACAAGTCCATCTTAAGAGACAAGTGTGCg              >  8:133645/1‑139 (MQ=255)                                                                                                                       gtTGATCAAGTCACCAGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAGTCGTCCATGTCAATGCTGACAACAACGATGACAAGTCCTTCTTCATAGAGAATTCTGCGATCa          >  5:216014/1‑139 (MQ=255)                                                                                                                        ttGATCAAGTCACCAGATTCTCTCATGTACACAAATCCTCTGGAAATCAACTCAGGACCCGCTGAAATCTTGAAGTCGTCCATGTCAATGCTGACAACAACGATGACAAGTCCTTCTTCAGAGAGAATTCTGCGATCAc         <  4:4601/139‑1 (MQ=255)                                                                                                                        ttGATCAAGTCACCAGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAGTCGTCCATGTCAATGCTGACAACAACGATGACAAGTCCTTCTTCAGAGAGAATTCTGCGATCAc         >  1:52729/1‑139 (MQ=255)                                                                                                                                gTCACCAGATTCTCTCATGTACACAAATCCGCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAGTCGTCCATGTCAATGCTGCCCACAACGATGACAAGTCCTTCGTCAGAGAGAATTCTGCGATCACGAAGTAc  >  4:99794/1‑138 (MQ=255)                                                                                                                                      |                                                                                                                                    TTTTTTCGTACAAGAAAGGTGCAAGCGTGTCTGTAATTTCGTTTTTGATTTCAGACCATTGAGTCGTTTTTCGTTCCATTACTTTTTGTAAGTGATTTGAAATGAGCTCTTGAGCGTCGTTGATCAAGTCACCAGATTCTCTCATGTACACAAATCCTCTGGAAATCAAATCAGGACCCGCTGAAATCTTGAAGTCGTCCATGTCAATGCTGACAACAACGATGACAAGTCCTTCTTCAGAGAGAATTCTGCGATCACGAAGTAC  >  CP006881/1486200‑1486464

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