Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,424,538	Δ1 bp	8.8%	coding (89/915 nt)	PP_4764 →	histone deacetylase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,424,538	0	A	.	8.8%	167.0 / 11.0	57	coding (89/915 nt)	PP_4764	histone deacetylase family protein
Reads supporting (aligned to +/- strand):  ref base A (29/23);  new base . (2/3);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

CCATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAG  >  NC_002947/5424393‑5424680                                                                                                                                                  |                                                                                                                                                 ccATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGac                                                                                                                                             >  2:731632/1‑150 (MQ=255) ccATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGac                                                                                                                                             <  2:695645/150‑1 (MQ=255) ccATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGac                                                                                                                                             >  2:514901/1‑150 (MQ=255) ccATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGac                                                                                                                                             >  1:947891/1‑150 (MQ=255)  cATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGC‑‑CGacca                                                                                                                                           >  1:562593/1‑150 (MQ=255)    ttGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGaccac                                                                                                                                          <  1:1132600/150‑1 (MQ=255)    ttGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGaccac                                                                                                                                          <  2:1400941/150‑1 (MQ=255)    ttGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGC‑‑CGaccag                                                                                                                                          <  2:1070282/149‑2 (MQ=255)            aCCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCg‑g                                                                                                                                 >  2:56610/1‑150 (MQ=255)             cccGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATTTACAAGAAAGATACCACCCCGGATTTCCCCGCTGAACACCGCTTCCCATTGGCCAAGTT‑CGCCTTCTTCACCGACCACCTGGTCgtc                                                                                                                                 >  2:1399737/1‑148 (MQ=255)              ccGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGC‑‑CGACCAGGTGGTCgtgca                                                                                                                               <  2:562593/150‑5 (MQ=255)                gTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACAGTGATGACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGc                                                                                                                             <  1:246635/150‑1 (MQ=255)                   aTGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCCGCTGC‑‑CGACCAGGTGGTCGTGCAGCagg                                                                                                                          <  1:364019/150‑1 (MQ=255)                        aaGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTc                                                                                                                                    <  2:1062327/136‑1 (MQ=255)                        aaGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTc                                                                                                                                    >  1:1062327/1‑136 (MQ=255)                               gCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGAc                                                                                                              >  1:645166/1‑150 (MQ=255)                                  cTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAg                                                                                                           <  1:514901/150‑1 (MQ=255)                                                gACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGGCAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGccc                                                                                             >  1:178276/1‑150 (MQ=255)                                                gACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGACAGGATACACTCTCTGAACACCGCTGCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGccc                                                                                             >  1:739685/1‑150 (MQ=255)                                                  ccccATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCgg                                                                                           >  1:104240/1‑150 (MQ=255)                                                   cccATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGa                                                                                          >  2:171641/1‑150 (MQ=255)                                                    ccatccatGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGAc                                                                                         >  2:143955/1‑150 (MQ=255)                                                      atccatGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCttgc                                                                                                    <  1:149845/137‑1 (MQ=255)                                                      atccatGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCttgc                                                                                                    >  2:149845/1‑137 (MQ=255)                                                      atccatGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACAt                                                                                       >  2:789609/1‑150 (MQ=255)                                                           tGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGcc                                                                                  <  2:668175/150‑1 (MQ=255)                                                                           cTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGc                                                                  <  1:56610/150‑1 (MQ=255)                                                                                aCGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGGCAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTgg                                                             <  2:178276/150‑1 (MQ=255)                                                                                aCGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGaccac                                                                                                                                          <  1:583262/74‑1 (MQ=255)                                                                                aCGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGaccac                                                                                                                                          >  2:583262/1‑74 (MQ=255)                                                                                aCGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTgg                                                             >  1:720790/1‑150 (MQ=255)                                                                                  gAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGcc                                                           >  2:230968/1‑150 (MQ=255)                                                                                     gATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCAt                                                        >  2:1256250/1‑150 (MQ=255)                                                                                             cccGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAg                                                <  2:104240/150‑1 (MQ=255)                                                                                             cccGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAg                                                <  1:171641/150‑1 (MQ=255)                                                                                                      ccccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGa                                       <  2:947891/150‑1 (MQ=255)                                                                                                        ccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGaccac                                                                                                                                          >  2:1213551/1‑50 (MQ=255)                                                                                                        ccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGaccac                                                                                                                                          <  1:1213551/50‑1 (MQ=255)                                                                                                         cGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CTACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgcg                                    <  2:567236/150‑1 (MQ=255)                                                                                                              aaCACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACa                               >  2:1320689/1‑150 (MQ=255)                                                                                                                        cccGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGa                     <  1:666706/150‑1 (MQ=255)                                                                                                                          cGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGc                   >  2:1001214/1‑150 (MQ=255)                                                                                                                           gATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCt                  >  1:346293/1‑150 (MQ=255)                                                                                                                                 cAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgcg            >  1:362000/1‑150 (MQ=255)                                                                                                                                  aaGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgcgc           >  1:334397/1‑150 (MQ=255)                                                                                                                                    gTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGa         <  2:362000/150‑1 (MQ=255)                                                                                                                                     ttCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAg        <  2:645166/150‑1 (MQ=255)                                                                                                                                       ccGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACAt                                          <  2:252518/114‑1 (MQ=255)                                                                                                                                       ccGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACAt                                          >  1:252518/1‑114 (MQ=255)                                                                                                                                         gCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGAcc    >  2:1186513/1‑150 (MQ=255)                                                                                                                                         gCCCGCTGC‑‑CGACCAGGTGGTCGTGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATc                                                                       >  2:774062/1‑83 (MQ=25)                                                                                                                                          ccTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGc                                                  >  2:1087460/1‑103 (MQ=255)                                                                                                                                          ccTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGc                                                  <  1:1087460/103‑1 (MQ=255)                                                                                                                                          ccTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCa   <  1:1001214/150‑1 (MQ=255)                                                                                                                                           ctgctgCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAg  >  1:162131/1‑150 (MQ=255)                                                                                                                                              ctgcA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACAt                                          <  2:779208/107‑1 (MQ=255)                                                                                                                                              ctgcA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACAt                                          >  1:779208/1‑107 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CCATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCA‑CGACCACCTGGTCG‑GCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAG  >  NC_002947/5424393‑5424680

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