Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	401,822	C→A	7.3%	intergenic (+349/+80)	PP_0334 → / ← PP_0335	transposase/glutathione S‑transferase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	401,822	0	C	A	7.3%	166.2 / 3.8	55	intergenic (+349/+80)	PP_0334/PP_0335	transposase/glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base C (25/26);  new base A (2/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.97e-01

TTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATAC  >  NC_002947/401680‑401956                                                                                                                                                |                                                                                                                                       ttCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCAGCGAACACACATAATGGGGTGACAG‑GCACCGCCTTAGCCCGTGTTCGCTGGCACCCCCg                                                                                                                                   >  2:180984/1‑148 (MQ=37)    cTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGCGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTc                                                                                                                                <  8:21197/148‑1 (MQ=255)     tATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCTGAGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTcc                                                                                                                               >  1:89997/1‑149 (MQ=255)        gTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCaa                                                                                                                            >  5:7530/1‑148 (MQ=255)         ttCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTc                                                                                                                                <  8:116569/143‑1 (MQ=255)         ttCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTc                                                                                                                                >  7:116569/1‑143 (MQ=255)         ttCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCaaaa                                                                                                                          >  4:126129/1‑149 (MQ=255)              aaaaGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCcca                                                                                                                     <  5:156406/149‑1 (MQ=255)                      cTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGAt                                                                                                             <  1:75422/149‑1 (MQ=255)                                tgGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACt                                                                                                   >  6:79778/1‑149 (MQ=255)                                 gggCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAacacag                                                                                                  >  1:66002/1‑147 (MQ=255)                                   gCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACtgtg                                                                                                >  3:157243/1‑149 (MQ=255)                                         ccATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGTCACGCCCCCTCCCAAAAGCACACCACAGATAGAGAACACTGTGGaaaa                                                                                           >  8:23991/1‑148 (MQ=255)                                              cataGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCTGTTTTCGCGGGCACGCCCGCTTCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTg                                                                                      <  5:189530/148‑1 (MQ=255)                                               ataGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTg                                                                                      >  8:8208/1‑147 (MQ=255)                                               ataGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTa                                                                                    <  3:126129/149‑1 (MQ=255)                                                  gACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGTACCCTGTGGAAAACCTTGGAg                                                                                   >  4:215572/1‑147 (MQ=255)                                                  gACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAg                                                                                   >  3:41619/1‑147 (MQ=255)                                                      gtgtGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCgg                                                                             <  8:74185/149‑1 (MQ=255)                                                      gtgtGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCgg                                                                             <  5:79778/149‑1 (MQ=255)                                                       tgtgGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCgg                                                                             <  8:74189/148‑1 (MQ=255)                                                          gggCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCg                                                                                                         >  8:41344/1‑117 (MQ=255)                                                          gggCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCg                                                                                                         <  7:41344/117‑1 (MQ=255)                                                            gCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAT‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGACCACCACAGATAGCGAACAATGTGGAAAACCTTGTAGGAGCGGGTTTag                                                                       <  8:127883/149‑2 (MQ=255)                                                             cGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTAcc                                                                      <  1:180984/149‑1 (MQ=255)                                                               aGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCcg                                                                    <  4:7825/149‑1 (MQ=255)                                                                gTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCcgc                                                                   >  5:109711/1‑149 (MQ=255)                                                                      cttttgggagcgggagtgcccgcgaacaccggcgaagccggtgcctgacaccccattatgtGTGTTCGCGGGCAAGCCCGCTCCTAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAgg                                                              <  7:78268‑M1/88‑1 (MQ=255)                                                                        ataGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGc                                                           <  3:243498/149‑1 (MQ=255)                                                                         taGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGcc                                                          >  5:177851/1‑149 (MQ=255)                                                                         taGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGcc                                                          <  6:177851/149‑1 (MQ=255)                                                                           gggagcgggcgtgcccgcgaacaccggcgaagccggtgcctgacaccccattatgtGTGTTCGCGGGCAAGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCggg                                                                            <  8:109769‑M1/74‑1 (MQ=255)                                                                           gggagcgggcgtgcccgcgaacaccggcgaagccggtgcctgacaccccattatgtGTGTTCGCGGGCAAGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCggg                                                                            >  7:109769‑M1/57‑130 (MQ=255)                                                                               cGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAgg                                                    <  1:100504/149‑1 (MQ=255)                                                                                gggCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCTGGTTTACCCGCGAAGGGGCCCTCAGGc                                                   <  2:89997/149‑1 (MQ=255)                                                                                gggCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGc                                                   <  6:109711/149‑1 (MQ=255)                                                                                 ggCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGc                                                   >  8:30573/1‑148 (MQ=255)                                                                                  gCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGACACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCa                                                 >  8:118/1‑149 (MQ=255)                                                                                            gAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGt                                       <  2:270644/149‑1 (MQ=255)                                                                                                  acaTAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGAc                                 <  6:7530/149‑1 (MQ=255)                                                                                                   cggcgaagccggtgcctgacaccccattatgtGTGTTCGCGGGCAAGCCCGCTCCTAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGAc                                 >  8:78268‑M1/33‑149 (MQ=255)                                                                                                    aTAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCa                               >  5:231643/1‑149 (MQ=255)                                                                                                    aTAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCa                               <  4:41619/149‑1 (MQ=255)                                                                                                        tGGGGTGTCAG‑GCAGCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGa                           <  7:23991/149‑1 (MQ=255)                                                                                                         ggggTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGAt                          >  3:28970/1‑149 (MQ=255)                                                                                                            gtgtCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGc                       >  5:198958/1‑149 (MQ=255)                                                                                                             tgtCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCggg                                                                            >  7:252463/1‑95 (MQ=255)                                                                                                             tgtCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCggg                                                                            <  8:252463/95‑1 (MQ=255)                                                                                                               cgcgaacACCGGCGAAGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGGAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGCCCATCGATGGCCg                     >  5:150318/7‑149 (MQ=255)                                                                                                                 aG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCAATGGCGCTGTAGCGACCATCGATGGCCGCtt                  >  1:71296/1‑149 (MQ=255)                                                                                                                          gCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCt          >  8:217567/1‑149 (MQ=255)                                                                                                                              ctccGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCgg       <  7:17154/146‑1 (MQ=255)                                                                                                                              cGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCg        <  7:30573/147‑1 (MQ=255)                                                                                                                               gCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCa     <  4:157243/149‑1 (MQ=255)                                                                                                                                   gtgtTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATAc  >  2:256916/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                       TTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAG‑GCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATAC  >  NC_002947/401680‑401956

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