Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,392,071	G→C	10.7%	R186P (CGG→CCG)	PP_2097 →	sensory box protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,392,071	0	G	C	10.7%	135.4 / 5.7	47	R186P (CGG→CCG)	PP_2097	sensory box protein
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base C (3/2);  total (24/23)
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.95e-01

GTGACCAGGACGCGCAAGGCCAATGGCGCCCGGCTTCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGCAGAGCGCGGCCTGCTGATGGTGGCGCCGGTATTTTCCGATGCCGACCCCCGCGGTGCGGCAGTGGGCTATGTGA  >  NC_002947/2391924‑2392220                                                                                                                                                     |                                                                                                                                                       gTGACCAGGACGCGCAAGGCCAATGGCGCCCGGCTTCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCggg                                                                                                                                                      >  2:36958/1‑150 (MQ=255)   gACCAGGACGCGCAAGGCCAATGGCGCCCGGCTTCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCAAGGAGGATCTGGCGCGgccg                                                                                                                                                    >  2:28814/1‑150 (MQ=255)      cAGGACGCGCAAGGCCAATGGCGCCCGGCTTCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCATGCTGCCC‑GGCGCCagcgccgg                                                                                                                                                >  1:1012115/1‑150 (MQ=255)      cAGGACGCGCAAGGCCAATGGCGCCCGGCTTCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCCggcgcca                                                                                                                                                 <  2:939210/150‑3 (MQ=255)      cAGGACGCGCAAGGCCAATGGCGCCCGGCTTCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCCggcgcca                                                                                                                                                 >  1:939210/1‑148 (MQ=255)      cAGGAAGAGCCAGGCAAATGGCGACCGGCTTCGCAGCTCGACCATTACTGCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGggcgctg                                                                                                                                                 <  2:511169/150‑1 (MQ=255)                  ggCCAATGGCGCCCGGCTTCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGggcgctggcgccggcgctg                                                                                                                                     <  2:629764/150‑1 (MQ=255)                    ccAATGGCGCCCGGCTTCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCc‑ggcgccagcgccggcgccagcgcccgcg                                                                                                                                  <  1:265007/150‑6 (MQ=255)                     cAATGGCGCCCGGCTTCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGggcgctggcgccggcgctggcg                                                                                                                                  >  2:1075261/1‑150 (MQ=255)                     cAATGGCGCCCGGCTTCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGggcgctggcgccggcgctggcg                                                                                                                                  >  2:1013972/1‑150 (MQ=255)                                  cTTCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCATGCTGCCC‑GGCGCCAGCGCCGGCGCCAGCGCCCGCGCCAGGGCCGCCTgg                                                                                                                    <  2:1012115/150‑1 (MQ=21)                                   ttCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGc                                                                                                                   <  1:255282/150‑1 (MQ=255)                                    tCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCCATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCg                                                                                                                  <  2:1029409/150‑1 (MQ=255)                                    tCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGGGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCg                                                                                                                  <  2:95572/150‑1 (MQ=255)                                    tCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCT‑AACCTTTTTGGGTCAAAGCGAGCCCCAGGCGGGTATGGCGCTTGCTCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCg                                                                                                                  >  1:519365/1‑150 (MQ=255)                                         gcgcGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGc                                                                                                                   >  1:1066249/1‑144 (MQ=255)                                         gcgcGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGc                                                                                                                   <  2:1066249/144‑1 (MQ=255)                                         gcgcGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTc                                                                                                             <  1:107082/150‑1 (MQ=255)                                                        cccGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTgg                                                                                                  >  1:823944/1‑146 (MQ=255)                                                        cccGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTgg                                                                                                  <  2:823944/146‑1 (MQ=255)                                                         ccGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTg                                                                                             >  2:782303/1‑150 (MQ=255)                                                            gTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTc                                                                                          <  2:1117017/150‑1 (MQ=255)                                                                 ctATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATAc                                                                                     >  2:792708/1‑150 (MQ=255)                                                                  tataCCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGCGCCCTTGGCCCTGTTCGATAcc                                                                                    >  1:156201/1‑150 (MQ=255)                                                                       ccAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGGGTCCGAGCCATTGGCCCTGTTCGATACCCAGTc                                                                               >  2:108314/1‑150 (MQ=255)                                                                        cAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCg                                                                              >  2:1099665/1‑150 (MQ=255)                                                                         aaTCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTc                                                                                                             >  1:461055/1‑118 (MQ=255)                                                                         aaTCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTc                                                                                                             <  2:461055/118‑1 (MQ=255)                                                                          aTCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGc                                                                            <  2:206209/150‑1 (MQ=255)                                                                             gagcgagcGGCCCGGGTTGCCTTACGGCCTGGACCTCGAT‑AGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGCaga                                                                         <  1:851261/150‑1 (MQ=255)                                                                                  agcgGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGCAGAGcgcg                                                                    <  1:792708/150‑1 (MQ=255)                                                                                     ggCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGCAGAGCGCGGcc                                                                 >  1:358012/1‑150 (MQ=255)                                                                                      gCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCAt                                                                                                     <  1:118219/113‑1 (MQ=255)                                                                                      gCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCAt                                                                                                     >  2:118219/1‑113 (MQ=255)                                                                                                 cTTACGGCCTGGACCTCGCT‑GGTCAAAGCGACCCCCACGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCAGAGCCATTGGCCCTGTTCGATACCCAGTCGGCAAAGCGCGTCCTGCTGAtggttg                                                     >  2:255522/1‑148 (MQ=255)                                                                                                   tACTGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGCAGAGCGCGGCCTGCTGATGGTGgcg                                                   >  2:706649/1‑150 (MQ=255)                                                                                                          tGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGCAGAGCGCGGCCTGCTGATGGTGGCGCCGGTAt                                            >  1:675879/1‑150 (MQ=255)                                                                                                           ggACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGt                                                                                            <  2:1019011/101‑1 (MQ=255)                                                                                                           ggACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGt                                                                                            >  1:1019011/1‑101 (MQ=255)                                                                                                              ccTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGa                                                                                                          >  1:890910/1‑84 (MQ=255)                                                                                                              ccTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGa                                                                                                          <  2:890910/84‑1 (MQ=255)                                                                                                                   cT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGCAGAGCGCGGCCTGCTGATGGTGGCGCCGGTATTTTCCGATg                                   <  1:1075261/150‑1 (MQ=255)                                                                                                                        tCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATa                                                                                      <  2:1051081/95‑1 (MQ=255)                                                                                                                        tCAAAGCGAGCCCCAGGCGGCCCTGGCGCCGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATa                                                                                      >  1:1051081/1‑95 (MQ=255)                                                                                                                          aaaGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGCAGAGCGCGGCCTGCTGATGGTGGCGCCGGTATTTACCGATGCCGAcc                             >  1:260339/1‑150 (MQ=255)                                                                                                                           acgCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGCAGAGCGCGGCCTGCTGATGGTGGCGCCGGTATTTTCCGATGCCGAccc                            <  1:28814/148‑1 (MQ=255)                                                                                                                            aGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGCAGAGCGCGGCCTGCTGATGGTGGCGCCGGTATTTTCCGATGCCGAcccc                           >  2:677006/1‑150 (MQ=255)                                                                                                                                    ccAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGCAGAGCGCGGCCTGCTGATGGTGGCGCCGGTATTTTCCGATGCCGACCCCCGCGGTGc                   <  2:735111/150‑1 (MQ=255)                                                                                                                                                ggCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGCAGAGCGCGGCCTGCTGATGGTGGCGCCGGTATTTTCCGATGCCGACCCCCGCGGTGCGGCAGTGGGCTa       <  2:358012/150‑1 (MQ=255)                                                                                                                                                   gcgGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGCAGAGCGCGGCCTGCTGATGGTGGCGCCGGTATTTTCCGATGCCGACCCCCGCGGTGCGGCAGTGGGCTAtgt    <  1:1099665/150‑1 (MQ=255)                                                                                                                                                     gggCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGCAGAGCGCGGCCTGCTGATGGTGGCGCCGGTATTTTCCGATGCCGACCCCCGCGGTGCGGCAGTGGGCTATGTGa  <  1:36958/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       GTGACCAGGACGCGCAAGGCCAATGGCGCCCGGCTTCGCAGCGCGACCATTACTTCCCGGTGCTCTATACCCAATCGAGCGAGCGGCCCGGGTTGCCTTACGGCCTGGACCTCGCT‑GGTCAAAGCGAGCCCCAGGCGGCCCTGGCGCGGGCGCTGGCGCCGGCGCTGGCGCC‑GGGCAGCATGGCGGTGTCCGAGCCATTGGCCCTGTTCGATACCCAGTCGGCAGAGCGCGGCCTGCTGATGGTGGCGCCGGTATTTTCCGATGCCGACCCCCGCGGTGCGGCAGTGGGCTATGTGA  >  NC_002947/2391924‑2392220

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