Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,122,363	C→G	8.2%	G20A (GGC→GCC)	PP_2739 ←	sensory box protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,122,363	0	C	G	8.2%	148.1 / 3.9	49	G20A (GGC→GCC)	PP_2739	sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (24/21);  new base G (2/2);  total (26/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 = 7.08e-01

CGGCTGGTCGTGATCCTCGCCCTGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCATGATGAACCAAGGAGG  >  NC_002947/3122221‑3122506                                                                                                                                                |                                                                                                                                                 cGGCTGGTCGTGATCCTCGCCCTGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGt                                                                                                                                            >  7:254954/1‑149 (MQ=255)       gTCGTGATCCTCGCCCTGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTACTGCTCGGCGACAAC‑GATGCCATCGTTGGAAt                                                                                                                                      >  6:109115/1‑149 (MQ=255)             aTCCTCGCCCTGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGa                                                                                                                                >  5:248759/1‑149 (MQ=255)             aTCCTCGCCCTGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGa                                                                                                                                <  4:86228/149‑1 (MQ=255)                cTCGCCCTGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCa                                                                                                                             >  1:106632/1‑149 (MQ=255)                      cTGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGtt                                                                                                                       <  3:264065/149‑1 (MQ=255)                       tGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCagatcgg                                                                                                                       >  7:31641/1‑142 (MQ=255)                       tGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTg                                                                                                                      >  7:57514/1‑149 (MQ=255)                       agaagaaATCGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGAc                                                                                                                               <  8:31641/139‑1 (MQ=255)                         aagaaAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGc                                                                                                                     <  8:61402/148‑1 (MQ=255)                         aagaaAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGc                                                                                                                     <  2:40799/148‑1 (MQ=255)                                   cAGTCCTGATAGAGAATATGGTCGGCGCAGTAGCCGTTCAGGCGCTCGAAGGCCGGTTTGACGTTGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAACTGATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTg                                                                                                           <  1:225289/149‑1 (MQ=255)                                     gTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑Gt                                                                                                       <  6:248759/149‑1 (MQ=255)                                        cttaTAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGa                                                                                                    <  8:41136/146‑1 (MQ=255)                                         tGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGAt                                                                                                   >  2:28270/1‑149 (MQ=255)                                         tGATAGAGAATATCGTCGGCGCAGTAGCCGGTAAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGAt                                                                                                   >  8:50034/1‑149 (MQ=255)                                               agaATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATggg                                                                                             >  5:271004/1‑149 (MQ=255)                                                   tatCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCAAAGCCATACGGCGGGGGGACAAC‑CGAGCCATGGTTGGAAAGTTGGGCCATAATTTGCAGGGGGTTTTC‑GTTGATGATgggggtg                                                                                         >  1:94981/1‑146 (MQ=255)                                                     tcgtcgGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCa                                                                                        <  2:106632/148‑1 (MQ=255)                                                        tcgGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTa                                                                                    <  7:50034/149‑1 (MQ=255)                                                                      ggTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTAt                                                                      <  2:147973/149‑1 (MQ=255)                                                                            gcgcTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGTATTTCGTATGCGAtg                                                                <  5:256464/149‑1 (MQ=255)                                                                            gcgcTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGTATTTCGTATGCGAtg                                                                >  6:256464/1‑149 (MQ=255)                                                                              gcTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGAtgtg                                                              >  2:242372/1‑149 (MQ=255)                                                                                      gCCGGGTTGACGTAGATAAGGATGCTCTCAATGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAg                                                      >  1:27159/1‑149 (MQ=255)                                                                                               aCGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCtt                                              <  8:254954/148‑1 (MQ=255)                                                                                                cGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGTCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTg                                             <  5:109115/148‑1 (MQ=255)                                                                                                 gTATATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGTGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGc                                            <  7:43944/148‑1 (MQ=255)                                                                                                       aaGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCt                                     >  2:236330/1‑149 (MQ=255)                                                                                                              ctctCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGgcg                              >  2:308614/1‑149 (MQ=255)                                                                                                                 tcATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGgcgc                             <  6:271004/147‑1 (MQ=255)                                                                                                                 tcATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCg                           >  5:134313/1‑149 (MQ=255)                                                                                                                 tcATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCg                           >  2:187862/1‑149 (MQ=255)                                                                                                                   aTTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAATGCTTCCTTGCAGAAGCTTTTGGCTCCGtt                         >  3:205932/1‑149 (MQ=255)                                                                                                                     tGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTttg                       >  3:60516/1‑149 (MQ=255)                                                                                                                       ccTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTttgtt                     <  1:187862/149‑1 (MQ=255)                                                                                                                       ccTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTttgtt                     <  1:28270/149‑1 (MQ=255)                                                                                                                       ccTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTAGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTttgtt                     <  8:57514/149‑1 (MQ=255)                                                                                                                        cTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTttgttg                    >  3:82831/1‑149 (MQ=255)                                                                                                                         ttCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTttgttg                    >  8:297661/1‑148 (MQ=255)                                                                                                                            cTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGCTGTTGATCATGGGGCTCACCGACCGGCATTTAGTATTCGAAGTGCCCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCat                 >  1:249509/1‑149 (MQ=255)                                                                                                                                  ggCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCATGATGAAc          <  4:60516/149‑1 (MQ=255)                                                                                                                                    cGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCATGATGAAcc         >  4:215382/1‑148 (MQ=255)                                                                                                                                       cAAC‑GATGGCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGt                          <  4:103048/128‑1 (MQ=255)                                                                                                                                       cAAC‑GATGGCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGt                          >  3:103048/1‑128 (MQ=255)                                                                                                                                       cAAC‑GATGGCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGt                          <  2:95435/128‑1 (MQ=255)                                                                                                                                       cAAC‑GATGGCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGt                          >  1:95435/1‑128 (MQ=255)                                                                                                                                        aaC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCATGATGAACCAagga    <  1:242372/149‑1 (MQ=255)                                                                                                                                          c‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCATGATGAACCAaggagg  >  1:66517/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                 CGGCTGGTCGTGATCCTCGCCCTGAAGAAAACGGCAGTCCTGATAGAGAATATCGTCGGCGCAGTAGCCGGTCAGGCGCTCGAAGGCCGGGTTGACGTAGATAAGGATGCTCTCATTGCCTTCCTGCTCGGCGACAAC‑GATGCCATCGTTGGAATGTTCGACCATCAGTTGCAGGAGTTTTGC‑GTTGATCATGGGGCTCACCTACAGGCATTTCGTATGCGATGTGCGCCAAAGGCTTCCTTGCAGAAGCTTTTGGCGCCGTTTGTTGCATGATGAACCAAGGAGG  >  NC_002947/3122221‑3122506

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