Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,813,272	0	C	T	43.6%	47.0 / 26.1	55	S181L (TCG→TTG)	dfp	fused 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (7/24);  new base T (24/0);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 1.11e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGC  >  NC_000913/3813134‑3813416                                                                                                                                           |                                                                                                                                                   agcaTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCaa                                                                                                                                           >  2:349113/1‑148 (MQ=255)      aaTTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACG‑ACCTg                                                                                                                                   >  2:518754/1‑150 (MQ=255)      aaTTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTg                                                                                                                                   <  1:464096/150‑1 (MQ=255)      aaTTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACATg                                                                                                                                   <  2:491627/150‑1 (MQ=255)           agagGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCCTCAACG‑ACCCG‑aaa                                                                                                                               >  1:527581/1‑148 (MQ=255)           agagGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACa                                                                                                                             <  1:55546/150‑1 (MQ=255)                  cttgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTTGTGCATTTTTTTCCCGCCACCC‑ACCCC‑CCACATCTGaa                                                                                                                       >  1:90535/1‑149 (MQ=255)                   ttgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTCGCCGGCAACC‑ACCTC‑AAACATCTgta                                                                                                                       >  2:157151/1‑146 (MQ=255)                     gcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTa                                                                                                                  <  1:188615/151‑1 (MQ=255)                        tCCCGTGGTTTGCTCATCTGTGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAAcatcat                                                                                                                 <  1:143192/149‑1 (MQ=255)                          ccGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGAtt                                                                                                             <  2:662464/151‑1 (MQ=255)                           cGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTa                                                                                                            <  1:676000/151‑1 (MQ=255)                            gTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTTCCCCTCAACG‑ACCCG‑AAAAATCTGAACATTATGAtt                                                                                                             >  1:255899/1‑149 (MQ=255)                              ggTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTTGCGGCTTTTTTGCCCCTCCACG‑ACCCT‑AAAAAACTGGAAATTTttattatc                                                                                                          >  1:169420/1‑148 (MQ=255)                                   gCTCATCTGGGGGCCAGACAGTGGCCGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAAGGAATTTTTTGCCCGGCCACGCTCCTG‑AAA‑ATTTTAACAttcttaattgaccc                                                                                                       >  1:470253/1‑136 (MQ=255)                                   gCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCgg                                                                                                     <  1:184729/150‑1 (MQ=255)                                       aTCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCCTCAACG‑ACCTG‑AAACATCTGAACAATATGATTAACCCCGGccc                                                                                                  >  2:353933/1‑149 (MQ=255)                                        tCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGCCAACG‑ACCTG‑AAAAATCTGAACATTATGATTACCGCCGGccc                                                                                                  >  2:267092/1‑148 (MQ=255)                                           gggggCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGAc                                                                                               <  1:354478/148‑1 (MQ=255)                                                ccAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAAc                                                                                       <  2:169442/151‑1 (MQ=255)                                                 cAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCATTTTTTGCCCCCCAAAC‑ACCTG‑AAAAAACTGAACAATATGATTACCGGCGGCCCGAAGCGGGaa                                                                                        >  1:641558/1‑149 (MQ=255)                                                 cAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTTCCCGTCAACG‑ACCCG‑AAAAATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAAcc                                                                                      >  2:316710/1‑151 (MQ=255)                                                    aCAGTGGCAGTCAGGTTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCTTTACCCATTGTGGATTTGGCGGTAGCGCATTTTTCCCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGc                                                                                    <  2:85223/150‑1 (MQ=255)                                                        tGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACC‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGAcgcggggaacgccccct                                                                               >  1:473055/1‑138 (MQ=255)                                                         ggccgtggggCTTGTGCTGATTGCGGTTCTGGGCGAATGCTCGATCCGTTAACCAGTGTGGATATCGCTGTCGGTAACTTGTCGCCCGTCAGCG‑CCCTG‑AAACTGCTCACAATTGTGTTTACCGCCGGCCAGACGCGTTCACCGCTTGa                                                                                <  1:193011/141‑1 (MQ=255)                                                          gCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATcc                                                                             <  1:349090/151‑1 (MQ=255)                                                            aGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCATTTTTTGCCCCTCCACG‑CCCTG‑TAACATCTGAACATTATTATTACCGCCGGCCCCACCCGTGCACCGCTCGATcc                                                                             >  1:707840/1‑149 (MQ=255)                                                            aGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCCTCAACG‑ACCCG‑AAACATCTTAACATTATGATTACCCCCCGCCCCACGCGTGAACCGCCCGATCCgg                                                                           >  1:425641/1‑151 (MQ=255)                                                               cAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTg                                                                         <  2:40810/150‑1 (MQ=255)                                                               cAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATT‑CCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCg                                                                       <  2:538309/151‑1 (MQ=255)                                                                      gtggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACG‑ACCCG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACCCGTGGACCGCCCGATCCCGTGCGGTatat                                                                 >  2:638745/1‑151 (MQ=255)                                                                        ggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCATTTTTTGCCCGTCAAAG‑ACCTG‑AAACATCTGAACATTATGATTTACGCCCGCCCCACGCGTGAACCCCTCGATCCGGTGCGTTATTtct                                                               >  2:35994/1‑151 (MQ=255)                                                                             tatCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCCCGATCCGGGGCGGTTTATCTCTaa                                                           >  2:278552/1‑150 (MQ=255)                                                                             tatCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCATTTTTCGCCCGTCAACC‑ACCTG‑AAAAATCTGAACAATATGATTACCGCCGGCCCCACACGTGAAACGCTCCATCCGGGGCGTTATATCCCTaa                                                           >  2:58112/1‑150 (MQ=255)                                                                              atCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAAt                                                          <  2:288400/150‑1 (MQ=255)                                                                                cGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGACGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATca                                                        <  2:743419/150‑1 (MQ=255)                                                                                 ggTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTTCCCGTCAACG‑ACCTG‑AAAAATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCCCGATCCCGTGCGGTTTATCTCTAAt                                                          >  1:399567/1‑147 (MQ=255)                                                                                 ggTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGCAGCGCATTTTTTTCCCCTCCACG‑ACCCG‑AAACATCTGAACATTATTATTACCGCCCGCCCGACGCGTGAACCGCCCGCTCCGGGGCGGTTTATCTCTaa                                                           >  1:480446/1‑146 (MQ=255)                                                                                      tGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGGGAACCGCTCGATCCGGTGCGTTTTATCTCTAATCACAGCTcc                                                 >  2:248869/1‑151 (MQ=255)                                                                                       gggCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTcc                                                 <  1:367871/150‑1 (MQ=255)                                                                                       gggCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTcc                                                 <  2:692710/150‑1 (MQ=255)                                                                                          cGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGACTTTTTCGCCCGTCCAAG‑AACCG‑AAAAATCTTAAAACTATGAATACCGCCCCCCCGCGACGCGGAACGCTCGATCACGAGCGTTATATTTCTAATCACAGGTCCgg                                               >  1:85199/1‑149 (MQ=255)                                                                                           gAATGCTCGATCCGTTAACCATTGTGGATATGGTGGTAGCGCATTTTTCGCCCGTCAACG‑CCCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCgg                                               <  2:52919/148‑1 (MQ=255)                                                                                           gAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAAACGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCgg                                               >  2:554181/1‑148 (MQ=255)                                                                                               gCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCCTTTTTTGCCCGTCAACC‑ACCTT‑ACACATCTGAAAATTATGATTAACCCCGGCCCGACGCGTGGAACCCTCTATCCGGGGCGTTATATCTCTTATTACAACTCCGGCCAGAt                                         >  1:770261/1‑150 (MQ=255)                                                                                                cTCGATCCGTTAACCATTGTGGATATGGCGGTATCGCCTTTTTCCCCCGTCAACA‑ACCTGTTAACAT‑TGTAAAATATGTTTAACCTCCGCTCGCCGCCGTAACCGCTCGCCACGGTTGGTTTTTTTTCTAATCACAggctcccctagtt                                         >  2:346651/1‑137 (MQ=255)                                                                                                cTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTTCCCGCCAACA‑ACCTG‑AAACATCTGAACATTATTATTACCGCCGGCCCGCCCCGGGGACCGCTCGATCCGGTGCGTTATATTTCTAATCACAACCCCGGCAAGATgg                                       >  1:470426/1‑151 (MQ=255)                                                                                                     tCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAAAC‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGGACCCCCCCATCCCGGGCGGTATATTTCTAATCAAAGCCCCGGCAAGATGGGgttt                                  >  1:367978/1‑151 (MQ=255)                                                                                                       cGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGc                                <  1:432303/151‑1 (MQ=255)                                                                                                          tAACCATTGTGGATATGGCGGTAGCGCCTTTTTCGCCCGCCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTAt                             >  2:502546/1‑151 (MQ=255)                                                                                                                 tgtgGATATGGCGGTAGCGCATTTTTTGCCCCTCAACG‑ACCTG‑AAAAATCTGAACATTATGATTACCCCCGGCCCCACGCGGGAACCGCCCGATCCGGGGCGTTTTATTTCTAATTACAGCCCCGGCCAGATGGGGTTTGCTATcgccccc                      >  1:410952/1‑151 (MQ=255)                                                                                                                     gATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑CCCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATcgccgccgccg                  <  1:329390/151‑1 (MQ=255)                                                                                                                              gTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTCTTGCTATCGCCGCCGCCGCTGCCcgtc         <  2:433160/151‑1 (MQ=255)                                                                                                                                aGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGACCATTATGGTTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCcgtcg        <  1:163815/150‑1 (MQ=255)                                                                                                                                     aTTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGgcgc  >  1:779166/1‑151 (MQ=255)                                                                                                                                           |                                                                                                                                                   AGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACG‑ACCTG‑AAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGC  >  NC_000913/3813134‑3813416

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