Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	783,388	G→T	11.4%	T847T (ACC→ACA)	PP_0672 ←	sensory box protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	783,388	0	G	T	11.4%	114.7 / 3.0	35	T847T (ACC→ACA)	PP_0672	sensory box protein
Reads supporting (aligned to +/- strand):  ref base G (21/10);  new base T (2/2);  total (23/12)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

TGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCTTCGTCGTCCAGTACCGCGGTGATGCCGACCCAG  >  NC_002947/783270‑783520                                                                                                                       |                                                                                                                                     tGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCgtagta                                                                                                                     >  2:60279/1‑136 (MQ=255)      ggTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCa                                                                                                                <  1:130918/136‑1 (MQ=255)                    gaggtcgaggAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATgcgc                                                                                                  >  2:9441/1‑136 (MQ=255)                     aggtcgaggAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCt                                                                                                 >  2:419556/1‑136 (MQ=255)                              aaCATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTgg                                                                                        >  1:181059/1‑136 (MQ=255)                                  tCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCtt                                                                                    <  2:25963/136‑1 (MQ=255)                                   cAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTg                                                                                   >  2:59881/1‑136 (MQ=255)                                           cccAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCt                                                                           <  2:41821/136‑1 (MQ=255)                                              aGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAggg                                                                                                                               >  2:594778/1‑81 (MQ=255)                                              aGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAggg                                                                                                                               <  1:594778/81‑1 (MQ=255)                                                  cTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACCGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCg                                                                    <  1:575875/136‑1 (MQ=255)                                                        ccGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGTCTTGCGTTCGCTGATGTCGGTgaag                                                              >  1:558352/1‑136 (MQ=255)                                                        ccGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTAGGTgaag                                                              >  1:48802/1‑136 (MQ=255)                                                                 ctgctgCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACg                                                     <  1:60279/136‑1 (MQ=255)                                                                             gTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAgg                                         >  1:301106/1‑136 (MQ=255)                                                                                   cAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCt                                   <  1:419556/136‑1 (MQ=255)                                                                                   cAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCt                                   <  1:55353/136‑1 (MQ=255)                                                                                         gTCCTGGAACAGTGTGCGGTTGGGCAGGTGTGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCTtcgtcg                             >  1:161315/1‑136 (MQ=255)                                                                                         gTCCTGGAACAGTGTGCGGTTGGGCAGGTGTGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCTtcgtcg                             >  1:164578/1‑136 (MQ=255)                                                                                         gTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCTtcgtcg                             >  2:61865/1‑136 (MQ=255)                                                                                              ggAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCTTCGTCGTCCAg                        >  2:62775/1‑136 (MQ=255)                                                                                                       tgCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCTTCGTCGTCCAGTACCGCGGt               >  2:268104/1‑136 (MQ=255)                                                                                                       tgCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCTTCGTCGTCCAGTACCGCGGt               >  1:123415/1‑136 (MQ=255)                                                                                                       tgCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCTTCGTCGTCCAGTACCGCGGt               >  1:642236/1‑136 (MQ=255)                                                                                                       tgCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCTTCGTCGTCCAGTACCGCGGt               >  1:447529/1‑136 (MQ=255)                                                                                                         cGGTTGGGCAGGTGTGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCTTCGTCGTCCAGTACCGCGGTGa             <  2:164578/136‑1 (MQ=255)                                                                                                         cGGTTGGGCAGGTGTGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCTTCGTCGTCCAGTACCGCGGTGa             <  2:161315/136‑1 (MQ=255)                                                                                                          ggTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCt                                   <  1:672628/113‑1 (MQ=255)                                                                                                          ggTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCt                                   >  2:672628/1‑113 (MQ=255)                                                                                                                  aGGTGGGTAAGGGCGTCGTAGTATGCCAGGCGTTGGATGCGATGTACGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGCAGCCCACGTAGCTGGCCAGGTCGCCTTCGTCGTCGAGTACCGCTGTGATGCAGAacc    >  2:156685/1‑136 (MQ=255)                                                                                                                  aGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTTGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTg                                               >  2:331281/1‑93 (MQ=255)                                                                                                                  aGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTTGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTg                                               <  1:331281/93‑1 (MQ=255)                                                                                                                  aGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCTTCGTCGTCCAGTACCGCGGTGATGCCGAccc    >  2:230963/1‑136 (MQ=255)                                                                                                                  aGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCTTCGTCGTCCAGTACCGCGGTGATGCCGAccc    >  1:587279/1‑136 (MQ=255)                                                                                                                    gTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCTTCGTCGTCCAGTACCGCGGTGATGCCGACCCAg  >  2:648273/1‑136 (MQ=255)                                                                                                                       |                                                                                                                                     TGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGTAAGGGCGTCGTAGTAGGCCAGGCGGTGGATGCGCTGTTCGCTGGCCTTGCGTTCGCTGATGTCGGTGAAGAAGCACACGTAGCTGGCCAGGTCGCCTTCGTCGTCCAGTACCGCGGTGATGCCGACCCAG  >  NC_002947/783270‑783520

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