Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,625,248	C→T	D118D (GAC→GAT)	serB →	phosphoserine phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,625,248	0	C	T	87.1%	66.8 / 1.3	32	D118D (GAC→GAT)	serB	phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (14/13);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

TCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGT  >  NC_000913/4625002‑4625493                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      ccGTAGCGGCTGGCTGCTGAATGGTCGTGGGCTGGATAAACATCGTCTGACCCAATCCCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCGGGTGCGTGAAAGATTAGCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGACACGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTT‑CTGGTGATGGACATGGATTCCa                                                                                                                                                                                                                                                   <  2:76472/250‑1 (MQ=255)   gTAGCGGCTGTCTGCTGTATGGTCGTGGGCGGGATAAACAGCGGTTGCCCCCATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCAcc                                                                                                                                                                                                                                                 <  2:314099/251‑1 (MQ=255)                             gggCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGTCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGAGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGCTGGAAATGGACTCCACCGCCATCCAGATTGAATGTATtggtgg                                                                                                                                                                                                                       >  2:163712/1‑247 (MQ=255)                             gggCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGAGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATtgatg                                                                                                                                                                                                                        <  1:163712/250‑1 (MQ=255)                                    aTAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGcc                                                                                                                                                                                                                <  1:284465/251‑1 (MQ=255)                                                                  aaCTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTATGGCAGGTTCACTCACCGCCCGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGG‑CAAATCCCGCACCTGCGCACGCCGGGCTTG‑CTGGGGAGGGATATAGACTCCACCGCCATCCAGATTGCATGTATTGCTGAAATTGCCAAACTGGCCGGGACGGGCgagttgtgagcg                                                                                                                                                                                  >  2:161433/1‑242 (MQ=255)                                                                      gtgTGCCGCGATGGTTATTGTTGTCGCCTGTTTCGTGTAAGATTTTCCTGTGATTCGTCTTGCAGGTGCACTATCCGCCCGGGCTTCGCGCCTGGCCCAGGAATAGCCTCTGGTTGTCGCCTCGATGTGG‑AAAATCCCGCACCTGCGCACGCCGGGTGTG‑CTGTTGATGGACATTTATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCATACTGGCCGGAACGGGCGAGAtgttgtcgta                                                                                                                                                                                <  2:324834/247‑9 (MQ=255)                                                                      gggTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTACCAGGTGATTCGTCTGGCAGGTGCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAg                                                                                                                                                                              <  1:56970/251‑1 (MQ=255)                                                                               gATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCAATGAAGCGCAACTGGATGTCGCCCCGCTGGGC‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGCAACCGAAc                                                                                                                                                                     >  1:364080/1‑251 (MQ=255)                                                                                  ggTGTTTGTTGCCGCCTGGGGCGTGGAAGATTATCAGGGTATTCGGCTGGCAGGTTCACTCCCCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCTCCTGCGCACGCCGGGTTTG‑CTGGTGAGGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACggg                                                                                                                                                                  <  2:378211/251‑1 (MQ=255)                                                                                     gATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCg                                                                                                                                                                >  1:118467/1‑250 (MQ=255)                                                                                                   ggTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTc                                                                                                                                                 <  1:273034/251‑1 (MQ=255)                                                                                                       cGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGAtt                                                                                                                                             >  1:291949/1‑251 (MQ=255)                                                                                                                    cAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCt                                                                                                                                >  1:46679/1‑251 (MQ=255)                                                                                                                                                     cGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTg                                                                                                >  1:108742/1‑250 (MQ=255)                                                                                                                                                                 cTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGAAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGGGGCGACGCTGAAAGGCGCTGACGCCACTATTCt                                                                                   >  1:276405/1‑251 (MQ=255)                                                                                                                                                                        agaaggagcactgggATGCCGCCCGGTTGGGG‑AAACCCCCCCCCCTGCGCACGC‑TGGTTTGTCTGGTGATGGACATGAATTCCCCCGCCATCCAGATTGAATGTATGGATAAAATTGCCAAACTGGCCGGAACGGGCGAGTGGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAg                                                                            <  2:69716/238‑1 (MQ=255)                                                                                                                                                                                  acTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGaa                                                                   >  1:82004/2‑250 (MQ=255)                                                                                                                                                                                     gtaTGTCGCCCCG‑TGGGGTAAAATCCCGCTCCGGCGCCCGCCGGGTTTG‑CTGGTGATGGGCATGGATTCCCCCGCCATCCAGAGTGATTGTATTGATGAAATTGCCAAACGGGCCGGAACGGGCGAGAGGGTGGCGGAAGTAACCGAACGGCCGATGCGCGGTGACCTCGATTTTCCCTCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGGGAAAATc                                                               <  1:242192/249‑1 (MQ=255)                                                                                                                                                                                       aTGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTg                                                             >  1:281338/1‑251 (MQ=255)                                                                                                                                                                                                    ggaggatttcccctactcTGCGCACGTCTGGTTTG‑CTGGTGAGTGACATGATATCCGTCGCCATCCAGATTGAATCTTTGGATAAAATTGCCAAATGGGCCGGACCGGGCGTGTTGGTGGCGGAAGTACCCGAACTGGCGTTGCGCTGCGAACTCGATTTTACCGCCATCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGGGCGTGAAAATCTGCCGCTGATGCCa                                                 <  2:10222/234‑1 (MQ=255)                                                                                                                                                                                                     ggg‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCAAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGGGGCGACGCTGAAATGCGCTGACGCCAATATGCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGGCAgg                                               >  2:287385/1‑251 (MQ=255)                                                                                                                                                                                                               ccGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGc                                                                                                                                                                 >  1:156430/1‑128 (MQ=255)                                                                                                                                                                                                               ccGCACCTGCGCACGCCGGGTGTG‑CTGGGGATGGACATGGATGCCACCGCCAGCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGATATGGTGGCGGAAGTAACCGAACGGGc                                                                                                                                                                 <  2:156430/128‑1 (MQ=255)                                                                                                                                                                                                                cGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCa                                      >  1:110475/1‑250 (MQ=255)                                                                                                                                                                                                                    ccTGCGCACGCCGGGTTTG‑CTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTgg                                 <  1:331062/251‑1 (MQ=255)                                                                                                                                                                                                                                           gTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgctggg           >  1:263125/1‑251 (MQ=255)                                                                                                                                                                                                                                            tGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAACATCTGCCGCTGCTGCCCGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgctgggc          >  2:198825/1‑251 (MQ=255)                                                                                                                                                                                                                                                  acatGGATTCCACGGCCAGCCAGATTGAATGTATGGATGAAATTGCCAAACTGGCCGGAACGGGAGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTGACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGaaa    <  1:281777/249‑1 (MQ=255)                                                                                                                                                                                                                                                   catGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACgg                                                                                                                                                                                                <  1:20148/61‑1 (MQ=255)                                                                                                                                                                                                                                                   catGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACgg                                                                                                                                                                                                >  2:20148/2‑62 (MQ=255)                                                                                                                                                                                                                                                    atGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGCAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCCGCCTGCGCAGCCGCGTGGCGTCGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGATAGt  >  1:18625/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      TCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTG‑CTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGT  >  NC_000913/4625002‑4625493

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