Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,671,562	C→T	7.4%	P123S (CCG→TCG)	PP_3234 →	HSP20 family heat shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,671,562	0	C	T	7.4%	155.3 / 2.5	54	P123S (CCG→TCG)	PP_3234	HSP20 family heat shock protein
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base T (2/2);  total (28/26)
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.74e-01

TGGCGGCCCTTCCAGCAATTGCGCCGCCAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCGTCAATGGCAACCTGTTGATCCGGGGCGAGAAAAAGGAGGAAGTGGAC  >  NC_002947/3671421‑3671699                                                                                                                                              |                                                                                                                                           tgGCGGCCCTTCCAGCAATTGCGCCGCCAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGTCGGCCAt                                                                                                                                     >  8:315256/1‑149 (MQ=255) tgGCGGCCCTTCCAGCAATTGCGCCGCCAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGTCGGCCAt                                                                                                                                     >  8:310330/1‑149 (MQ=255)   gcggcCCTTCCAGCAATTGCGCCGCCAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATg                                                                                                                                    <  2:282007/148‑1 (MQ=255)      gcCCTTCCAGCAATTGCGCCGCCAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACa                                                                                                                                <  3:291843/149‑1 (MQ=255)        ccTTCCAGCAATTGCGCCGCCAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACAt                                                                                                                               >  5:48758/1‑148 (MQ=255)          ttCCAGCAATTGCGCCGCCAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAg                                                                                                                            >  4:29486/1‑149 (MQ=255)          ttCCAGCAATTGCGCCGCCAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAg                                                                                                                            >  3:230722/1‑149 (MQ=255)           tCCAGCAATTGCGCCGCCAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCCGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGc                                                                                                                           >  2:57196/1‑149 (MQ=255)             cAGCAATTGCGCCGCCAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCCGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCg                                                                                                                          <  1:57196/148‑1 (MQ=255)              aGCAATTGCGCCGCCAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAg                                                                                                                        >  3:83927/1‑149 (MQ=255)                       gccgccAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGc‑c                                                                                                                 >  1:167052/1‑146 (MQ=255)                       gccgccAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGa                                                                                                               <  2:167052/148‑1 (MQ=255)                        ccgccAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAgg                                                                                                             <  6:78/149‑1 (MQ=255)                          gccAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGaa                                                                                                           <  4:230722/149‑1 (MQ=255)                              gATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATAcc                                                                                                       <  6:60863/149‑1 (MQ=255)                                       ccTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGc                                                                                                                   <  4:25781/129‑1 (MQ=255)                                       ccTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGc                                                                                                                   >  3:25781/1‑129 (MQ=255)                                           ttCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGcc                                                                                            >  7:3054/1‑147 (MQ=255)                                             cGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAg                                                                                         <  3:147568/148‑1 (MQ=255)                                              gAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGtt                                                                                       >  8:7720/1‑149 (MQ=255)                                                 gATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGcc                                                                                    >  7:19873/1‑149 (MQ=255)                                                   ttttGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTg                                                                                  >  2:95138/1‑149 (MQ=255)                                                   ttttGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTg                                                                                  <  4:83927/149‑1 (MQ=255)                                                    tttGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTgg                                                                                 >  5:17135/1‑149 (MQ=255)                                                        gACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCAATGTCGACCCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTg                                                                             <  2:72841/149‑1 (MQ=255)                                                         aCGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTg                                                                             >  5:175269/1‑148 (MQ=255)                                                           ggcggcCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGAt                                                                          >  3:46213/1‑149 (MQ=255)                                                               gcCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCa                                                                      <  6:48758/149‑1 (MQ=255)                                                                 cGATGCGCATGCCGTTCAGCCACACACGGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAg                                                                    <  8:19873/149‑1 (MQ=255)                                                                 cGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTATGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAg                                                                    >  1:74356/1‑149 (MQ=255)                                                                  gATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAg                                                                    >  3:185374/1‑148 (MQ=255)                                                                          tGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCgag                                                           <  4:185374/149‑1 (MQ=255)                                                                            ccGTTCATCCACACACCGTTCCATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCAAGGCATGCCGGCCATGGACATCATCGAGCTGGCTCGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCgaga                                                          <  1:95138/149‑1 (MQ=255)                                                                            ccGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGTCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCgaga                                                          <  7:310330/148‑1 (MQ=255)                                                                            ccGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGTCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCACCGAGTTGCCTGGCGTGGATGCCAAGGATATCgaga                                                          <  7:315256/148‑1 (MQ=255)                                                                                  aGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGc                                                    <  7:7720/148‑1 (MQ=255)                                                                                    ccACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGcc                                                                                    >  1:160276/1‑114 (MQ=255)                                                                                    ccACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGcc                                                                                    <  2:160276/114‑1 (MQ=255)                                                                                    ccACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCg                                                 >  1:289708/1‑149 (MQ=255)                                                                                            cGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATc                                                        >  6:333950/1‑134 (MQ=255)                                                                                            cGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATc                                                        <  5:333950/134‑1 (MQ=255)                                                                                                cGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCGTCAATGGCAAc                                      <  3:29486/148‑1 (MQ=255)                                                                                                cGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCGTCAATGGCAAc                                      <  2:289708/148‑1 (MQ=255)                                                                                                 gATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCGTCAATGGCAACCt                                    >  2:244776/1‑149 (MQ=255)                                                                                                      cGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCGTCAATGGCAACCTGTTGa                               >  8:232872/1‑149 (MQ=255)                                                                                                            gTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCGTCAATGGCAACCTGTTGATCCggg                         <  8:3054/149‑1 (MQ=255)                                                                                                             ttCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCGTCAATGGCAACCTGTTGATCCgggg                        >  7:16484/1‑149 (MQ=255)                                                                                                             ttCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCGTCAATGGCAACCTGTTGATCCgggg                        <  4:46213/149‑1 (MQ=255)                                                                                                                   cggcgTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCGTCAATGGCAACCTGTTGATCCggg                         <  1:314/142‑1 (MQ=255)                                                                                                                   cggcgTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCGTCAATGGCAACCTGTTGATCCggg                         >  2:314/1‑142 (MQ=255)                                                                                                                              gTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCGTCAATGGCAACCTGTTGATCCGGGGCGAGAAAAAGGAGGAAg       >  4:320501/1‑149 (MQ=255)                                                                                                                                tCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCGTCAATGGCAACCTGt                                  >  2:18531/1‑120 (MQ=255)                                                                                                                                tCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCGTCAATGGCAACCTGt                                  <  1:18531/120‑1 (MQ=255)                                                                                                                                    ccACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTAGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCGTCAATGGCAACCTGTTGATCCGGGGCGAGAAAAAGGAGGAAGTGGAc  >  3:186406/1‑148 (MQ=255)                                                                                                                                              |                                                                                                                                           TGGCGGCCCTTCCAGCAATTGCGCCGCCAGATCGACAGCCTGTTCGAGGATTTTGGACGGCGGCCGATGCGCATGCCGTTCAGCCACACACCGTTCGATGTCGAACCGTTCTGGCGGCGTGACCTGTTCACCCACGGCATGCCGGCCATGGACATCAGCGAGCTGGC‑CGAGGAATACCGCATCAGCGCCGAGTTGCCTGGCGTGGATGCCAAGGATATCGAGATCAAGCTCGTCAATGGCAACCTGTTGATCCGGGGCGAGAAAAAGGAGGAAGTGGAC  >  NC_002947/3671421‑3671699

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