Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,271,666	C→T	20.4%	intergenic (+50/‑50)	PP_t32 → / → PP_t33	tRNA‑Asp/tRNA‑Asp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,271,666	0	C	T	20.4%	91.2 / 15.0	44	intergenic (+50/‑50)	PP_t32/PP_t33	tRNA‑Asp/tRNA‑Asp
Reads supporting (aligned to +/- strand):  ref base C (16/19);  new base T (4/5);  total (20/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

GAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAA  >  NC_002947/2271529‑2271804                                                                                                                                          |                                                                                                                                             gAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCC‑CAAGGCATGTGCC‑TTgtggc                                                                                                                                 <  8:55925/149‑5 (MQ=21)   aCGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCC‑CAAGGCATGTGCC‑TTgt                                                                                                                                    >  1:208244/1‑143 (MQ=255)   aCGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCC‑CAAGGCATGTGCC‑TTgt                                                                                                                                    <  2:208244/144‑2 (MQ=255)             cAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGATTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGttt                                                                                                                     <  3:229960/149‑1 (MQ=255)             cAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTg‑gg                                                                                                                                  >  7:83022/1‑136 (MQ=255)             cAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTg‑gg                                                                                                                                  <  8:83022/136‑1 (MQ=255)                 ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACTGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAAt                                                                                                                 <  1:285612/149‑1 (MQ=255)                 ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAAt                                                                                                                 >  1:109170/1‑149 (MQ=255)                      ttCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCATGTGCC‑TTgtg                                                                                                                                   >  7:122528/1‑125 (MQ=255)                      ttCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCATGTGCC‑TTgtg                                                                                                                                   <  8:122528/127‑3 (MQ=255)                      ttCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGgc                                                                                                            <  4:59654/149‑1 (MQ=255)                           tCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTg‑g                                                                                                                                   <  1:77592/121‑1 (MQ=255)                           tCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTg‑g                                                                                                                                   >  2:77592/1‑121 (MQ=255)                           tCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGa                                                                                                                              <  2:30801/126‑1 (MQ=255)                           tCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGa                                                                                                                              >  1:30801/1‑126 (MQ=255)                                 aGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGaa                                                                                                 >  1:213900/1‑149 (MQ=255)                                  gAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGaa                                                                                                 <  2:213900/148‑1 (MQ=255)                                  gAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGaa                                                                                                 >  6:195492/1‑148 (MQ=255)                                    aTACCGGCCTGTAACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGa                                                                                              <  5:195492/149‑1 (MQ=255)                                       ccGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATAc                                                                                           <  2:109170/149‑1 (MQ=255)                                                 aCGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTa                                                                                  <  6:24032/148‑1 (MQ=255)                                                   gCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGtt                                                                               <  8:182467/149‑1 (MQ=255)                                                     cGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACACCGCACATGCC‑TTT‑GGCTGATCCCAGTTTCATTCTGGCGCAAGGCTGAACGATTCGCAGCGGTAGTTCa                                                                             <  6:286950/149‑1 (MQ=38)                                                           tCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGt                                                                       >  7:19586/1‑149 (MQ=255)                                                             gcgGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGt                                                                       <  4:168307/147‑1 (MQ=255)                                                             gcgGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGt                                                                       >  3:168307/1‑147 (MQ=255)                                                                                                       accggTGCCAGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTccc                           >  3:226812/4‑149 (MQ=255)                                                                                                       aTGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTccc                           <  8:19586/149‑1 (MQ=255)                                                                                                                    aCCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCa              >  7:140348/1‑149 (MQ=255)                                                                                                                         tCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTccc                           >  2:199700/1‑131 (MQ=255)                                                                                                                         tCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTccc                           <  1:199700/131‑1 (MQ=255)                                                                                                                            gAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCgg                                      >  7:96017/1‑117 (MQ=255)                                                                                                                            gAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCgg                                      <  8:96017/117‑1 (MQ=255)                                                                                                                                gcccaAGGCATGTGCCTTTT‑GGCTGATCCCTGTTTCAATCGGGCGCAAGCCTGAATGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTc                                  <  8:271492/115‑1 (MQ=11)                                                                                                                                ccACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTc    >  5:90348/1‑147 (MQ=255)                                                                                                                                ccACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTc    >  5:158686/1‑147 (MQ=255)                                                                                                                                ccACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCaa  <  8:140348/149‑1 (MQ=255)                                                                                                                                 cccaAGGCATGTGCC‑TTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAATGATACGCAGCGGTAGtt                                                                               <  8:106742/70‑1 (MQ=17)                                                                                                                                 cccaAGGCATGTGCC‑TTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTc    >  3:272000/3‑147 (MQ=255)                                                                                                                                 cccaAGGCATGTGCC‑TTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTc    >  2:199597/3‑147 (MQ=255)                                                                                                                                  acaAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTc         >  6:62428/1‑140 (MQ=255)                                                                                                                                  acaAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTc         <  5:62428/140‑1 (MQ=255)                                                                                                                                       gCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCaa  <  6:202256/142‑1 (MQ=255)                                                                                                                                       gCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCaa  >  5:202256/1‑142 (MQ=255)                                                                                                                                          |                                                                                                                                             GAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACAAGGCACATGCC‑TTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAA  >  NC_002947/2271529‑2271804

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