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

CCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGT  >  NC_000913/3813122‑3813422                                                                                                                                                        |                                                                                                                                                         ccGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTc                                                                                                                                                          <  1:813992/151‑1 (MQ=255)      gCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGccc                                                                                                                                                      <  2:581231/150‑1 (MQ=255)             agcaTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTTGCCCGTCaa                                                                                                                                                 >  1:772588/1‑148 (MQ=255)                aTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGGG‑TTTTTTCCCCCTCAaccacc                                                                                                                                            >  2:451613/1‑150 (MQ=255)                aTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCC‑TTTTTTTCCCCTCAACGAcc                                                                                                                                            >  1:50780/1‑150 (MQ=255)                   aTTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCC‑TTTTTTGCCCGTCAACGAcc                                                                                                                                            >  1:583454/1‑147 (MQ=255)                      tAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGGGCC‑ATTTTTGCCCGTCCACGAcccgcaa                                                                                                                                       >  2:344365/1‑144 (MQ=255)                       agagGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTTGCCCGTCAACGACCTGaaa                                                                                                                                       >  1:249183/1‑148 (MQ=255)                              cttgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCA‑TTTTTTGCCCGGCAACGACCTGAAACATCTGaa                                                                                                                               >  1:92001/1‑149 (MQ=255)                               ttgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCC‑TTTTTTGCCCGTCAACGACCTGAAACATCTGAACa                                                                                                                             >  1:557940/1‑150 (MQ=255)                                tgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCC‑TTTTTTTCCCGTTAACGACCTTAAACATTTGAAAA‑tt                                                                                                                          >  1:92104/1‑151 (MQ=255)                                   ttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGCGAGCGCA‑ACTTTTGCCCGGCAACGACCTGGAACTTCTGGACA‑tt                                                                                                                          >  2:674607/1‑148 (MQ=255)                                     cccGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGa                                                                                                                      <  2:602733/150‑1 (MQ=255)                                     cccGTGGTTTGATCATCTGGGGGCCAGACAGTGGCAGCCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTCGCGCA‑CTTTTTGCCCGCCCACCACCTCAAAccactttaaattatg                                                                                                                        >  2:530323/1‑135 (MQ=255)                                       cGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTa                                                                                                                   <  2:191781/151‑1 (MQ=255)                                                  cATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTTGCCCGTCAACGACCTGAAAAATCTGAACA‑TTATGATTAACGCCGGccc                                                                                                         >  1:844624/1‑150 (MQ=255)                                                       gggggCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCC‑TTTTTTGCCCGTCAAAGACCTGGAACATCTGAACA‑TTATGATTACCGCCGGccccgcgc                                                                                                    >  2:347194/1‑145 (MQ=255)                                                       gggggCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGAcgc                                                                                                    <  1:485691/150‑1 (MQ=255)                                                         gggCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCATTTTTTTGCCCG‑CAAAGAACTGAAACATCTGAACATTTTTGTTAACCCCCGGCCCgcgc                                                                                                     >  2:485277/1‑145 (MQ=255)                                                           gCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGa‑a                                                                                              <  2:853728/151‑1 (MQ=255)                                                                aCAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCC‑TTTTTTGCCCGGCAAACACCTGGAACAACTTGAAA‑ATATGATTAACCCCCGCCCGCCGCGGGG‑Acc                                                                                            >  2:404451/1‑148 (MQ=255)                                                                 cAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTc                                                                                        <  1:664423/151‑1 (MQ=255)                                                                         gTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCgg                                                                                 <  1:59896/150‑1 (MQ=255)                                                                         gTCAGGCTTGTGGTGATAACGGTCCTGGGCGAATGCTCGACCCTTTAACCATTGTGGATATGGCGGTAGCGCT‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGt                                                                                <  2:765408/151‑1 (MQ=255)                                                                             ggCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCg                                                                             <  2:277050/150‑1 (MQ=255)                                                                              gCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGtt                                                                           <  1:595642/151‑1 (MQ=255)                                                                              gCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGtt                                                                           <  1:574035/151‑1 (MQ=255)                                                                              gCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGtt                                                                           <  1:392757/151‑1 (MQ=255)                                                                                        atatCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATCTCTa                                                                  <  1:98413/150‑1 (MQ=255)                                                                                         tatCGTTCCTGGGCGAATTCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TAAGGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATctct                                                                   <  2:687980/148‑1 (MQ=255)                                                                                         tatCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAACGCA‑TTTTTTGCCCCTCAACGAACTGAAACAACCTAACA‑TTATTATTACCGCCCGCCCCGCGCGGGA‑ACCGCTCCATCCCGTGCGGTTTATctctta                                                                 >  1:162192/1‑148 (MQ=255)                                                                                          atCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAGCA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATCTCTAAt                                                                <  2:638411/150‑1 (MQ=255)                                                                                                        aaTGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAAAATCTGAACA‑TTATGATTACCCCCGGCCCGACGCGTGG‑ACCGCCCGATCCGGGGCGGTTTATCTCTAATCACAGCTCCgg                                                     >  1:69280/1‑147 (MQ=255)                                                                                                            cTCGTTCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATg                                              <  2:425936/150‑1 (MQ=255)                                                                                                            cTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCC‑TTTTTTTCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCCCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATgg                                             >  1:284933/1‑151 (MQ=255)                                                                                                              cGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATgg                                             <  2:49722/149‑1 (MQ=255)                                                                                                              cGATCCGTTAACCATTGTGGATATGGCGGGCGCGCC‑TTTTTTTGCCGTCTACGCCCCTCAACACCTGAACC‑ATTTTATTTCCCCCGGCCCGACTCGTGG‑CCCCCTCCAACAGGGGCGTTTTATCTCTAATCACCACGCCGGCAAGATg                                              >  2:171902/1‑148 (MQ=255)                                                                                                                aTCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTTCCCCGTCCAACAACTGGAACATCTGAACA‑TTATGATTTACGCCGGCCCCAACCGTGA‑AACGCCCCATCCGGGGGGGTTTATTTCTAAACACACCTCCCGGAAGATggg                                            >  2:2838/1‑148 (MQ=255)                                                                                                                   cGTTAACCATTGTGGATATGGCGGTAGCGCC‑TTTTTCGCCCCTCCAAGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGGCGCGTGA‑ACCGCTCCATCCGGGGCGGTATATCTCTAATCACAGCTCCGGCAAGATGGGtttt                                        >  2:345978/1‑149 (MQ=255)                                                                                                                   cGTTAACCATTGTGGATATGGCGGCTGCGCA‑CTTTTTGCCCCTCAACACCCTGAAACATCTGTACA‑TTTTGTTTAACGGCGGCCCCCAGCGGGAGACCG‑ACCATCCGGTGCGCTATTTCCCTTATCACAGCCCCGGCAAGATGTGGTTTg                                       >  2:608286/1‑150 (MQ=255)                                                                                                                         ccATTGTGGATATGGCGGTAGCGCA‑TTTTTTGCCCGTCAACGACCTGAAACATCTGAACA‑TTATTATTACCGCCCGCCCCACGCGTGA‑ACCGCCCGATCCCGGGCGGTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATcgc                                >  1:9547/1‑151 (MQ=255)                                                                                                                          cATTGTGGATATGGCGGTAGCGCA‑TTTTTTGCCCCGCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑AACGCCCGATCCGGGGCGTTATATCTCTAATCACAGCCCCGGCAAGATGGGGTTTTCTATtgcc                               >  1:441069/1‑151 (MQ=255)                                                                                                                          cATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATcgcc                               <  1:67229/151‑1 (MQ=255)                                                                                                                             tgtgGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATcgccgc                             <  2:583315/150‑1 (MQ=255)                                                                                                                                ggATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATcgccgccgc                          <  1:222454/150‑1 (MQ=255)                                                                                                                                 gATATGGCGGTAGCGCA‑TTTTGCCCCCGCCAACGACCTGAAACATCGGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATcgccgccgcc                         <  2:72906/150‑1 (MQ=255)                                                                                                                                    atGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATcgccgccgccgc                       >  1:14868/1‑149 (MQ=255)                                                                                                                                        cGGTAGCTCA‑TTTTTCGCCCGTCACCGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGccc                  <  2:92083/150‑1 (MQ=255)                                                                                                                                            aGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCcgtcg              >  2:700819/1‑150 (MQ=255)                                                                                                                                                   tttttCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGcgcg       >  2:521409/1‑151 (MQ=255)                                                                                                                                                       tCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAAc    <  1:702841/150‑1 (MQ=255)                                                                                                                                                        cGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCAGGCCCGACGCGTGA‑ACCGCTCGATCCGGCGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCAGCCGCCGCTGCCCGTCGTGGCGCGAACGt  >  2:490218/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                         CCGCTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCA‑TTTTTCGCCCGTCAACGACCTGAAACATCTGAACA‑TTATGATTACCGCCGGCCCGACGCGTGA‑ACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGT  >  NC_000913/3813122‑3813422

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