Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,030,279	G→T	8.9%	intergenic (+34/‑59)	PP_3553 → / → PP_3554	AMP‑binding domain‑containing protein/acyl‑CoA dehydrogenase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,030,279	0	G	T	8.9%	156.7 / 4.7	52	intergenic (+34/‑59)	PP_3553/PP_3554	AMP‑binding domain‑containing protein/acyl‑CoA dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base G (27/20);  new base T (3/2);  total (30/22)
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.90e-01

CTGGCCGGGTTCAAGGTGCCGAAGACCGTGGTGTTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGCATGACGTGTTCAACGCTTCGGGCCTGTGGGCCCGACTGCCCGCC  >  NC_002947/4030138‑4030427                                                                                                                                              |                                                                                                                                                      cTGGCCGGGTTCAAGGTGCCGAAGACCGTGGTGTTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTgcgc                                                                                                                                                >  3:274170/1‑149 (MQ=255)     ccGGGTTCAAGGTGCCGAAGACCGTGGTGTTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTg                                                                                                                                             >  2:119657/1‑148 (MQ=255)      cGGGTTCAAGGTGCCGAAGACCGTGGTGTTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCCGTGTGCTGCGCCTg                                                                                                                                             >  1:250624/1‑147 (MQ=255)        ggTTCAAGGTGCCGAAGACCGTGGTGTTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGa                                                                                                                                         <  7:113966/149‑1 (MQ=255)        ggTTCAAGGTGCCGAAGACCGTGGTGTTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGa                                                                                                                                         <  7:21618/149‑1 (MQ=255)             aaGGTGCCGAAGACCGTGGTGTTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCCGTGTGCTGCGCCTGAAggggagc                                                                                                                                    <  2:250624/149‑7 (MQ=38)                  gCCGAACACCACGGTCTTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTACAGGCGCAGCACACg                                                                                                                                                       <  4:50800/125‑1 (MQ=255)                  gCCGAACACCACGGTCTTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTACAGGCGCAGCACACg                                                                                                                                                       >  3:50800/1‑125 (MQ=255)                   ccGAAGACCGTGGTGTTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCAc                                                                                                                               <  6:140157/148‑1 (MQ=255)                       aGACCGTGGTGTTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGccc                                                                                                                           >  1:160303/1‑148 (MQ=255)                        gACCGTGGTGTTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCaa                                                                                                                         <  5:329299/149‑1 (MQ=255)                           cGTGGTGTTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAAccc                                                                                                                      >  7:245444/1‑149 (MQ=255)                             tggtgTTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACTGGGCCAAAGCCCTCTGATGCCATAGCCCAACCCCTACAGCCGCAGCCCACGCATCCGCGTGAACATAAAGCGCAC‑CCCAAccccac                                                                                                                   >  2:281451/1‑146 (MQ=255)                               gtgtTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAAccctccc                                                                                                                  <  8:245444/149‑1 (MQ=255)                                  ttCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCtt                                                                                                                <  4:274170/148‑1 (MQ=255)                                              cccAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAagag                                                                                                   >  6:19585/1‑149 (MQ=255)                                              cccAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAagag                                                                                                   <  7:303432/149‑1 (MQ=255)                                                cAAGACCTCGACCGGCAATATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCa                                                                                                                                     <  8:226158/113‑1 (MQ=255)                                                cAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCt                                                                                                 >  4:153909/1‑149 (MQ=255)                                                  aGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGc                                                                                               >  6:130732/1‑149 (MQ=255)                                                   ttccgcGACCGGCAACATCCATCAGTACGTGCTGCGGGCATTGGCCAAAGCCCGTTGATGCCATACCCGCTCCCCTTCAGGCGCAGCACACTGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGcc                                                                                              <  1:281451/144‑1 (MQ=255)                                                    aCCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGcc                                                                                              >  1:45390/1‑148 (MQ=255)                                                    aCCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCg                                                                                             <  5:130732/149‑1 (MQ=255)                                                      cTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAag                                                                                                     >  4:141219/1‑139 (MQ=255)                                                      cTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAag                                                                                                     <  3:141219/139‑1 (MQ=255)                                                              gCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACt                                                                                    >  8:29858/1‑148 (MQ=255)                                                               cAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTa                                                                                  >  7:198410/1‑149 (MQ=255)                                                                  gATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTc                                                                               >  3:167612/1‑149 (MQ=255)                                                                    tCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCg                                                                             >  6:47956/1‑149 (MQ=255)                                                                      cAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑c                                                                           <  2:45390/148‑1 (MQ=255)                                                                         aaGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGc                                                                                               >  3:228681/1‑126 (MQ=255)                                                                         aaGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGc                                                                                               <  4:228681/126‑1 (MQ=255)                                                                           gTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCt                                                                     >  2:97023/1‑149 (MQ=255)                                                                                      gcgACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCCGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCCACCCTCCCTTACTTTTGCATGAGCTGCCGATGTCCACTTATTCCG‑CCCCCCTGCGCGACATgc                                                          >  6:300474/1‑149 (MQ=255)                                                                                        gACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTCTATGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATgcgc                                                        <  1:119657/149‑1 (MQ=255)                                                                                          ccGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCtt                                                      >  5:173249/1‑149 (MQ=255)                                                                                           cGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCGCCCCCCCTGCGCGACATGCGCtt                                                      >  7:132381/1‑149 (MQ=255)                                                                                             ggCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGt                                                   >  2:133475/1‑149 (MQ=255)                                                                                                cAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCt                                                >  6:263435/1‑149 (MQ=255)                                                                                                  aaGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGc                                              <  7:29858/149‑1 (MQ=255)                                                                                                    gCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCAATTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGCAt                                            >  7:338792/1‑149 (MQ=255)                                                                                                     cccTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCAc                                                                                                                               >  5:106238/1‑66 (MQ=255)                                                                                                     cccTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCAc                                                                                                                               <  6:106238/66‑1 (MQ=255)                                                                                                          tGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGCATGACg                                        <  5:47956/147‑1 (MQ=255)                                                                                                           gATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGCATGACgtgt                                     >  4:78379/1‑149 (MQ=255)                                                                                                            aTGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCTCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGCATGACGTGtt                                    <  8:198410/149‑1 (MQ=255)                                                                                                            aTGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGCATGACGTGtt                                    <  6:173249/149‑1 (MQ=255)                                                                                                              gCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGCATGACGTGTTc                                   >  3:299450/1‑148 (MQ=255)                                                                                                                       gCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGCATGACGTGTTCAACGCTTCgg                         >  3:110335/1‑149 (MQ=255)                                                                                                                           cccTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGTCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGCATGACGTGTTCAACGCTTCGGGCCt                     >  8:298384/1‑149 (MQ=255)                                                                                                                             cTTCAGGCGCAGCCGTGTGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGCATGACGTGTTCAACGCTTCGGGCCtgt                   >  4:310017/1‑149 (MQ=255)                                                                                                                             cTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGCATGACGTGTTCAACGCTTCGGGCCtgt                   <  4:167612/149‑1 (MQ=255)                                                                                                                              ttCAGGCGCAGCCGTGTGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGCATGACGTGTTCAACGCTTCGGGCCtgtg                  >  1:345125/1‑149 (MQ=255)                                                                                                                                cAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGCATGACGTGTTCAACGCTTCGGGCCTGTgg                 <  4:299450/148‑1 (MQ=255)                                                                                                                                           ataggCTGCGCTTGAAGATTCTTCGATCTCACAAACTTCCCTTAATATTGCAAGAGATGTCGATGTCCACTTATTCCG‑CCCCCCTGAGCGACATGCGCTTCTTCCTGCATGACGTGTTAAAAGCTTCGGTCCTGTGGGCCCGACTgccc     <  8:338792/146‑1 (MQ=255)                                                                                                                                              ggCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGCATGACGTGTTCAACGCTTCGGGCCTGTGGGCCCGACTgcccgcc  <  4:110335/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                      CTGGCCGGGTTCAAGGTGCCGAAGACCGTGGTGTTCGGCGAGCTGCCCAAGACCTCGACCGGCAAGATCCAGAAGTACGTGCTGCGCGACCGGGCCAAAGCCCTCTGATGCCATAGCCGCTCCCCTTCAGGCGCAGCACACGGCTGCGCCTGAAGATTCAGCGCACGCCCAACCCTCCCTTACTATTGCAAGAGCTGCCGATGACCACTTATTCCG‑CCCCCCTGCGCGACATGCGCTTCGTCCTGCATGACGTGTTCAACGCTTCGGGCCTGTGGGCCCGACTGCCCGCC  >  NC_002947/4030138‑4030427

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