Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,391,920	G→C	15.6%	G106A (GGT→GCT)	orn →	oligoribonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,391,920	0	G	C	15.6%	73.2 / 5.5	32	G106A (GGT→GCT)	orn	oligoribonuclease
Reads supporting (aligned to +/- strand):  ref base G (15/12);  new base C (3/2);  total (18/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 = 9.55e-01

GATGACTGGAACGTGCGCACCCATACCGCCAGCGGGCTGGTAGAGCGCGTGAAAGCGAGCACGATGGGCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCGGAGCTGGAAGCCTACTTCCACTACCGTTATCTCGATGTCAGCACCCTGAAAGAGCTGGCGCGCCGCTGGAAGCCGGAAATTCTGGATGGTTTTACC  >  NC_000913/4391778‑4392062                                                                                                                                               |                                                                                                                                               gatgACTGGAACGTGCGCACCCATACCGCCAGCGGGCTGGTAGAGCGCGTGAAAGCGAGCACGATGGGCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCgg                                                                                                                                                >  1:292938/1‑143 (MQ=255)       tGGAACGTGCGCACCCATACCGCCAGCGGGCTGGTAGAGCGCGTGAAAGCGAGCACGATGGGCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATGGCCGATTTGCTGTAagag                                                                                                                                          >  1:610047/1‑143 (MQ=255)         gAACGTGCGCACCCATACCGCCAGCGGGCTGGTAGAGCGCGTGAAAGCGAGCACGATGGGCGCTCGGGAAACTGAACTGGCAACGCTCGAATTTTTAAAACCGTGGGTGCCTGCGCGAAAACCGCCGAGGTGCCGTCACAGCa                                                                                                                                        >  2:174953/1‑143 (MQ=255)              cgctctcccATACCGACAGTGTGTTTGTTGAGCTCGTGAAAGCGACCACGATGTGCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGt                                                                                                                                   <  1:640305/137‑1 (MQ=255)               gcgcACCCATACCGCCAGCGGGCTGGTAGAGCGCGTGAAAGCGAGCACGATGGGCGATCGGGCAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGCTAACAGCATCGGTc                                                                                                                                  >  2:725994/1‑143 (MQ=255)                     cttttctgCCAGCTGTCTGGTAGAGCGCGTGAAAGCGAGCACGCTGGGCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGAcc                                                                                                                            <  1:575622/136‑1 (MQ=255)                          ccGCCAGCGGGCTGGTAGAGCGCGTGAAAGCGAGCACGATGGGCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGCTTTGCGCTCGCAACGACGGTGAGGACcgtcgt                                                                                                                       >  2:3238/1‑143 (MQ=255)                          ccGCCAGCGGGCTGGTAGAGCGCGTGAAAGCGAGCACGATGGGCGATCGGGAAGCTGAACTGGCAACGCTCGAAGTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGCTAACAGCACCGGTCAGGACcgtcgc                                                                                                                       >  2:578158/1‑142 (MQ=255)                           ggCCAGCGGGCTGGTAGAGCGCGTGAAAGCGAGCACGATGGGCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGtt                                                                                                                      <  1:61208/142‑1 (MQ=255)                             ccAGCGGGCTGGTAGAGCGCGTGAAAGCGAGCACGATGGGCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGCCCGTCGTTTc                                                                                                                    >  1:705878/1‑143 (MQ=255)                             ccAGCGGGCTGGTAGAGCGCGTGAAAGCGAGCACGATGGGCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGGCCGGGCCGTCGATTc                                                                                                                    >  2:660025/1‑143 (MQ=255)                             ccAGCGGGCAGGAAGAGCGCGGGAAAGCGAGCACGATGGGCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTc                                                                                                                    <  1:13376/143‑1 (MQ=255)                                      tGGTAGTGCGCTTGAAAGCGAGCACGATGGGCGACCTGTAAGCGGAACTGGCAACGCTCGAATTTTTATAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTaaa                                                                                                           <  2:47810/143‑1 (MQ=255)                                        gTAGAGCGAGTGAAAGCGAGCACGATGGGAGATCGGGAAGCTGAACTGGCAACGCTCGAACTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATa                                                                                                         <  1:379654/143‑1 (MQ=255)                                          agagCGCGTGAAAGCGAGCACGATGGGCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGGTTAAATACa                                                                                                       >  2:575832/1‑143 (MQ=255)                                                    aaaGCGAGCACGATGGGCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCGGAGCt                                                                                             >  2:255871/1‑143 (MQ=255)                                                     aaGCGAGCACGATGGGCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACCTGCCGGAGcgg                                                                                            >  2:604199/1‑141 (MQ=255)                                                       gCGAGCACGATGGGCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCGGAGCTGGa                                                                                          >  2:164737/1‑143 (MQ=255)                                                          aGCACGATGGGGGATCGGGAAGCTGAACGGGCAACGCGCGAATTTTGAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCGGAGCTGGAAGc                                                                                       <  1:491788/143‑1 (MQ=255)                                                                 tGGGCGATGGGGAAGTTGAATTGGGAGCGCTCGAATTGTTAAAACAGTGGGTGACTGCGGGAAAATCGGCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCGGAGCTGGAAGCCTACTTc                                                                                <  1:745946/143‑1 (MQ=255)                                                                    gCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATGTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACCCGCAGGAGCTGGAAGCCTACTTcccc                                                                             >  2:280638/1‑141 (MQ=255)                                                                                cTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTTCCTGCGGGATATTCGCCGATTTGCGCTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCGGAGCTGGAAGCCTACTTCCACTACCGTTAtctc                                                                 <  1:573295/143‑1 (MQ=255)                                                                                  gagcTGGGAACGCTCGAATTTTTTAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCGGAGCTGGAAGCCTACTTCCACTACCGTTATCTCGa                                                               <  2:484150/140‑1 (MQ=255)                                                                                       tgtaaagCTCGAATTTTGAAAGGCGGGTGTGACTGCGGGAAAATCGGCGAGTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCGGAGCTGGAAGCCTACTTCCACTACCGTTATCTCGATGTCa                                                          <  1:136485/137‑1 (MQ=255)                                                                                               tCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCTGAGCTGGAAGCCTACTCCCACTACCGTTATCTCGATGTCAGCACCCTg                                                  >  1:757695/1‑143 (MQ=255)                                                                                                 gAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCGGAGCTGGAAGCCTACTTCCACTACAGTTATCTCGATGTCAGCACCCTGaa                                                >  2:272691/1‑143 (MQ=255)                                                                                                  aaTTTTTAAAACAGTGGGTGCCTGCTGGATTATCGTCGATTTGCGGTAACAGCAACGGGCAGGACCGGCGTTTCCTGTTTAAATACATGCCGGAGCTGGAAGCCTACTTCCACTACCGTTATCTCGATGTCAGCACCCTGaaa                                               <  1:135653/143‑1 (MQ=255)                                                                                                      tttAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCGGAGCTGGAAGCCTACTTCCACTACCGTTATCTCGATGTCAGCACCCTGAAACAGc                                           >  2:244234/1‑143 (MQ=255)                                                                                                      tttAAAACAGTGGGTGCCTGCGAGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGCTTAAATACATGCCGGAGCTGGAAGCCTACTTCCCACACCGCTCACTCGATGTCAGCACCCAGAAAGAGc                                           >  2:62926/1‑143 (MQ=255)                                                                                                                          cGGGAGAATCGCCGATTTGCGGTAACGGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACAGGCCGGAGCTGGAAGCCTACTTCCACTACCGTTATCTCGATGTCAGCACCCTGAAAGAGCTGGCGCGCCGCTGGAAGCCg                       <  1:464051/143‑1 (MQ=255)                                                                                                                                     ccGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCGGAGCTGGAAGCCTACTTCCACTACCGTTATCTCGATGTCAGCACCCTGAAAGCGCTGGCGCGCCGCTGGAAGCCGGACATTCgggc            >  1:244465/1‑142 (MQ=255)                                                                                                                                     ccGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCGGAACAGCAACACGACTTCCGCTACCGTTATCTCGAAATCCACACCCAAAAAGAACAGAGGCGCAGCTGGAAGCCGGGAATTCTggt            >  2:260195/1‑142 (MQ=255)                                                                                                                                          gtGCGCTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCGGAGCTGGAAGCCTACTTCCACTACCGTTATCGCGATGTCAGCACCCTGAAAGAGCTGGCGCGCCGCTGGAAGCCGGAAATTCTGGATGGtt       <  1:752742/142‑1 (MQ=255)                                                                                                                                               gTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCGGAGCTGGAAGCCTACTTCCACTACCGTTATCTCGATGTCAGCACCCTGAAAGAGCTGGCGCGCCGCTGGAAGCCGGAAATTCTGGATGGTTTTAcc  <  1:575832/143‑1 (MQ=255)                                                                                                                                               |                                                                                                                                               GATGACTGGAACGTGCGCACCCATACCGCCAGCGGGCTGGTAGAGCGCGTGAAAGCGAGCACGATGGGCGATCGGGAAGCTGAACTGGCAACGCTCGAATTTTTAAAACAGTGGGTGCCTGCGGGAAAATCGCCGATTTGCGGTAACAGCATCGGTCAGGACCGTCGTTTCCTGTTTAAATACATGCCGGAGCTGGAAGCCTACTTCCACTACCGTTATCTCGATGTCAGCACCCTGAAAGAGCTGGCGCGCCGCTGGAAGCCGGAAATTCTGGATGGTTTTACC  >  NC_000913/4391778‑4392062

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑