Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,920,270	A→G	9.4%	L499P (CTG→CCG)	acs ←	acyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,920,270	0	A	G	9.4%	139.1 / 5.3	53	L499P (CTG→CCG)	acs	acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base G (2/3);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

AATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATG  >  NC_002947/3920135‑3920410                                                                                                                                        |                                                                                                                                              aaTTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCGGCTGCACCCgggg                                                                                                                                  <  1:4780/149‑2 (MQ=255)         cACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCTGGGAGCGGCTGCACCCGGGGCatgct                                                                                                                            <  7:55253/147‑5 (MQ=255)               gCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCa                                                                                                                     <  6:234956/148‑1 (MQ=255)                    tcaccgtcaTCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TAc                                                                                                              <  7:200764/149‑1 (MQ=255)                           caTCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTtgctg                                                                                                        >  1:91245/1‑148 (MQ=255)                            aTCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCa                                                                                                                                     <  5:96753/119‑1 (MQ=255)                            aTCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCa                                                                                                                                     >  6:96753/1‑119 (MQ=255)                                       gtcgACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCaa                                                                                           >  7:64614/1‑149 (MQ=255)                                          gACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCa                                                                                        >  6:228046/1‑149 (MQ=255)                                           aCAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCAc                                                                                       <  8:53506/149‑1 (MQ=255)                                             aaaaCGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCt                                                                                     >  2:60791/1‑149 (MQ=255)                                               aaCGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCA‑TACTTGCTGCAGGGCACGCCAATCACCTGc                                                                                   <  4:72916/149‑1 (MQ=255)                                                aCGGAAATACCGGGGCACCTTGAAGTGCGCTATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTg                                                                                    >  1:30893/1‑147 (MQ=255)                                                aCGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTg                                                                                    <  2:30893/147‑1 (MQ=255)                                                aCGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGca                                                                                  >  7:112784/1‑149 (MQ=255)                                                        aCCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTg                                                                                    <  7:284448/139‑1 (MQ=255)                                                        aCCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTg                                                                                    >  8:284448/1‑139 (MQ=255)                                                          cGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGcc                                                                                             <  2:97341/128‑1 (MQ=255)                                                          cGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGcc                                                                                             >  1:97341/1‑128 (MQ=255)                                                          cGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCAc                                                                        >  6:154759/1‑149 (MQ=255)                                                                  cTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGc                                                                                   >  8:180054/1‑130 (MQ=255)                                                                  cTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCgggt                                                                >  6:236283/1‑149 (MQ=255)                                                                      aaGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCgggtggg                                                             <  8:64614/148‑1 (MQ=255)                                                                         tGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCa                                                                                        >  8:313036/1‑118 (MQ=255)                                                                         tGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCa                                                                                        <  7:313036/118‑1 (MQ=255)                                                                         tGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGcc                                                                          <  7:328629/132‑1 (MQ=255)                                                                         tGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGcc                                                                          >  8:328629/1‑132 (MQ=255)                                                                             cAATCCTTGCCCTCGCCCACTCATGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTgaagaa                                                      >  7:315044/1‑148 (MQ=255)                                                                              aaTCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTgaagaag                                                     >  1:34440/1‑148 (MQ=255)                                                                               aTCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTgaagaagaa                                                   <  5:228046/149‑1 (MQ=255)                                                                                 ccTTGTCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTgaagaag                                                     >  2:38535/1‑145 (MQ=255)                                                                                 ccTTGTCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTgaagaag                                                     <  1:38535/145‑1 (MQ=255)                                                                                  cTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAActc                                                >  4:27360/1‑149 (MQ=255)                                                                                               aCTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTcgcgc                                   >  4:327646/1‑149 (MQ=255)                                                                                                   aCGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAActcctc                                             <  4:91703/135‑1 (MQ=255)                                                                                                   aCGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAActcctc                                             >  3:91703/1‑135 (MQ=255)                                                                                                       aaCTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTTAATTGATGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTaaa                            >  1:286106/1‑149 (MQ=255)                                                                                                              ccTCACTGACAGCATGGCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTtctc                     <  7:156746/148‑1 (MQ=255)                                                                                                              ccTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCgggtgggt                                                            <  2:186040/109‑1 (MQ=255)                                                                                                              ccTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCgggtgggt                                                            >  1:186040/1‑109 (MQ=255)                                                                                                              ccTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCgg                                 <  8:41898/136‑1 (MQ=255)                                                                                                              ccTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCgg                                 >  7:41898/1‑136 (MQ=255)                                                                                                               cTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTcgc                   <  8:112784/149‑1 (MQ=255)                                                                                                                      aCAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTcgccgccgc             <  5:154759/148‑1 (MQ=255)                                                                                                                      aCAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTcgccgccgc             <  2:34440/148‑1 (MQ=255)                                                                                                                       cAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGa           <  3:327646/149‑1 (MQ=255)                                                                                                                              ccccGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATc     <  1:60791/148‑1 (MQ=255)                                                                                                                                ctgggTGCGGCTGCACCCGGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATg  >  1:226947/3‑149 (MQ=255)                                                                                                                                ctgggTGCGGCTGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTaa                             >  4:88817/3‑122 (MQ=255)                                                                                                                                ctgggTGCGGCTGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTaa                             <  3:88817/120‑1 (MQ=255)                                                                                                                                ccGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATg  >  8:288870/1‑149 (MQ=255)                                                                                                                                    gTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATg  >  4:191996/1‑145 (MQ=255)                                                                                                                                    gTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATg  <  3:191996/145‑1 (MQ=255)                                                                                                                                        |                                                                                                                                              AATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCG‑TACTTGCTGCAGGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATG  >  NC_002947/3920135‑3920410

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