Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,984,090	0	T	C	39.6%	76.9 / 16.5	53	F378S (TTC→TCC)	aslB	putative AslA‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base T (28/4);  new base C (1/20);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 9.75e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAA  >  NC_000913/3983951‑3984236                                                                                                                                             |                                                                                                                                                   tttAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTAc                                                                                                                                         >  2:367521/1‑151 (MQ=255)        aGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCat                                                                                                                                   >  1:567478/1‑150 (MQ=255)        aGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCat                                                                                                                                   <  1:675239/150‑1 (MQ=255)           tACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTNNGTGNNNNNNNTCAGCGTTATTTCCGCCATCTACCGCCATATCt                                                                                                                               >  2:119411/1‑151 (MQ=255)             ccGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGGGGCTGCCCGAAACCCCGCTTCATG‑CTCGATGCCAGCGGCAACCCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTa                                                                                                                             <  1:688685/151‑1 (MQ=255)             ccGGCGCAGTGTCCCAGTTGAAACGTTTTAAAAGCGTGGTGGGGGGGCTGCCCGAAACCCCGCTTCATG‑CGCGATGCCAGCGCCAACCCGGGACTAAATTATTTGTGTGCCGGGTATCAGCGTTATTCCCGCCTTCTACCGCCATATCtt                                                                                                                              <  2:39792/150‑1 (MQ=255)               ggCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTaa                                                                                                                            >  2:370987/1‑150 (MQ=255)                     gtgtCGCAGTTGTACCGTGTAAAAAGCGTGCTGGGGAGGCTGCCCGAAACCCCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATg                                                                                                                     <  1:870134/151‑1 (MQ=255)                     gtgtCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATg                                                                                                                     >  1:557059/1‑151 (MQ=255)                       gtCGCAGTTGTAACGTGTAAAAAGCGTGCTGGGGAGGCTGCCCAAAACCCCGCTTCATG‑CTCGATGCCAGCGCCAACCCGGGACTAAATTTTTTGTGTGCCGGGTATCAGCTTTATTCCCGCCTTCTACCGCCATATCTTAAAGCAATgg                                                                                                                    <  1:296126/150‑1 (MQ=255)                          gCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTg                                                                                                                 >  2:111299/1‑150 (MQ=255)                           cAGTTTTAACGTGTTAAAAGCGTGCTGGGGGGGCTGCCCGAACCCCCCTTTCATG‑CGCGACGCCTGCGGCAATCCGGGTCTAATTTATTTTTGTGCCGGGTATCAGCGTTATTCCCTCCATCTACCGCCATATCTTAAAGCAATGGc                                                                                                                   <  2:827733/147‑1 (MQ=255)                            agaagTAACGAGTTAAAAGCATGCTGCGGGAGATGTCCGAAATACTACTTAACG‑CCATATGCAAGCTGCAAACCGGGACTGAATTATTTGTGAGCCTTGTAACAGCCATATGCACGCCCTCTACCGACATATCCTAAAGCagagggtga                                                                                                                >  2:833025/5‑141 (MQ=255)                            aGTTGTACCGTGTTAAAGCCTTGCTGGGGAGGCTCCCCAAGCCACCG‑TTCCTGCCTCGTAGCCAGCGCCAAACCGGGACGAATTTTTTTTTGTGCCGGTTTTCATCGTTATTCCCGCCTTTTACCGCCATATCTTAAAGCAATGGc                                                                                                                   <  2:784748/146‑1 (MQ=255)                            aGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGAtt                                                                                                              >  2:710828/1‑151 (MQ=255)                              ttGTAACGTGTTAAAAGCGTGCTGGGGAGCCGCCCCGAAACCCCCCTTCATG‑CTCGATGCCAGCGCCAACCGGGAACTAAATTATTTGTGTGCGGGGTATCAGCGTTTTTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGAttt                                                                                                             <  2:974958/150‑1 (MQ=255)                                 taccttgtAAAAAGCGGGCGGGGGGGGCTGCCCGAACCCCCGCTTCAGG‑CTCGAGGCCAGCGCCAACCCGGGCCTGAATTATTTTTGTGCGGGTTATCATCTTTTTTTCCCCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCt                                                                                                          <  2:390732/145‑1 (MQ=255)                                  aaCGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTgg                                                                                                        >  2:1028135/1‑151 (MQ=255)                                   aCGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGgc                                                                                                       >  2:859244/1‑151 (MQ=255)                                    cGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGgcg                                                                                                      >  1:945655/1‑151 (MQ=255)                                     gtgtAAAAAGCGTGCGGGGGAGGCTCCCCGAACCCCCGTTTCATG‑CGCGAGCCCGGCGACAACCCGGGACTAAATTTTTTGTGGGCGGGGTTTCAGCTTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGgcg                                                                                                      <  2:136297/150‑1 (MQ=255)                                          aaaggctggcgggggATGCTGCCCAAAACCCCCCTTCATG‑CTCTGTGCCGGCTCCAACCAGGTACTAATTTTTTTGTGTCCCGGGTATCAGCTTTAGTCCCCCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGt                                                                                                <  2:564603/140‑1 (MQ=255)                                          aaaaGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGt                                                                                                >  2:575689/1‑151 (MQ=255)                                           aaaGCGTGCTGGGGAGGCTGCCCAAAACCCCGCTTCATG‑CTCGATGCCGGCGGCAACCCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTc                                                                                               <  1:121271/151‑1 (MQ=255)                                           aaaGCGTGCTGGGGAGGCGCCCCAAAACCCCGCTTCATG‑CTCGATGCAGGCGGCAACCCGGGACTAAATTTTTTGTGTGCCGGGTATCAGCGTTATTCCCGCCTTCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTc                                                                                               <  1:830713/151‑1 (MQ=255)                                               cGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGccc                                                                                           >  1:678214/1‑151 (MQ=255)                                               cGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGccc                                                                                           >  2:220978/1‑151 (MQ=255)                                                    tGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGcc                                                                                       >  2:160359/1‑150 (MQ=255)                                                          ggCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGAc                                                                                 >  1:706717/1‑150 (MQ=255)                                                               cccGAAACCCCGCTTCATG‑CTCGATGCCAGCGGCAACCCGGGACTGATTTATTTGTGTGCCGGGTATCAGCGTTTTTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACAtt                                                                              <  1:1065389/148‑1 (MQ=255)                                                               cccGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTatg                                                                           >  2:114375/1‑151 (MQ=255)                                                                   aaaCACCGCTCCATC‑CTCGCTCCCAGCGGCAAACCGGCACTGAATTATTGGTGTGCGGGGTTTCATCGTTTTTCCCCCCCTCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTatgc                                                                          <  2:260462/148‑1 (MQ=255)                                                                       acCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATgcgc                                                                    >  2:448394/1‑150 (MQ=255)                                                                           cTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTg                                                               >  2:955886/1‑151 (MQ=255)                                                                            ttCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGc                                                              >  1:44222/1‑151 (MQ=255)                                                                                 g‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCtggtg                                                         >  1:162029/1‑151 (MQ=255)                                                                                          ccggCGGAAAACCGGGAGTGAATTATTGGTGTGCGGGTTATCGGCGTTATTCCCGCCTTCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCtg                                                            <  2:1021960/137‑1 (MQ=255)                                                                                          ccAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGGGagta                                                  >  2:110828/1‑150 (MQ=255)                                                                                           cAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGagta                                                  >  1:595530/1‑149 (MQ=255)                                                                                                  aacccGGGACTGAATTTTTTGTGTGCCGGGTATCAGCTTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAAt                                         <  1:836633/148‑1 (MQ=255)                                                                                                  aaaCCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAAt                                         >  1:432805/1‑151 (MQ=255)                                                                                                      cGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATcggc                                     >  2:350506/1‑151 (MQ=255)                                                                                                       gggACTGAATTTTTTTTGTGCCGGTTATCAGCGTTATTCCCGCCATCTCCCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGGACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATcggcg                                    <  2:157705/151‑1 (MQ=255)                                                                                                       gggACTGAATTTTTTGTGTGCGGGGTAGCAGCGTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATcggc                                     <  1:452615/150‑1 (MQ=255)                                                                                                                   tttGTGTGCCGGGTATCAGCTTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCgg                         <  2:1039605/150‑1 (MQ=255)                                                                                                                     tgtgtgCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTg                      >  2:226527/1‑151 (MQ=255)                                                                                                                        gcgccgtttATCAGCTTGTTTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGAt                    <  2:366151/142‑1 (MQ=255)                                                                                                                         tgCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGAtt                   >  2:795024/1‑150 (MQ=255)                                                                                                                             gggTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCgg              >  1:179480/1‑151 (MQ=255)                                                                                                                                   cAGCTTTTTTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATgc          <  1:619976/149‑1 (MQ=255)                                                                                                                                     gtgTTATTCCCCCCATCTGCCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGAt            <  1:1028468/143‑1 (MQ=255)                                                                                                                                        ttATTCCCCCCATCAACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTa    <  2:508551/150‑1 (MQ=255)                                                                                                                                          aTTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTaaa  >  2:1042177/1‑150 (MQ=255)                                                                                                                                             |                                                                                                                                                   TTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATG‑CTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAA  >  NC_000913/3983951‑3984236

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