Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,580,143	C→T	T216T (ACG→ACA)	ydeM ←	putative anaerobic sulfatase maturation enzyme YdeM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,580,143	0	C	T	94.6%	113.8 / ‑3.7	37	T216T (ACG→ACA)	ydeM	putative anaerobic sulfatase maturation enzyme YdeM
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (18/17);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTAC  >  NC_000913/1579905‑1580387                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     gTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGcaca                                                                                                                                                                                                                                          <  2:498280/251‑1 (MQ=255)     aGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGaa                                                                                                                                                                                                                                       >  1:559949/1‑250 (MQ=255)       attaattaAAGTTTCCAATTTTGTACTGTGGATAGACAAAATGGGCGCATTCGTAAATGTATAAATTACTTTCAACAACCAGATTATCATTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTACCGGATGAAAGTTTCACCCACATCGATTTTAACCCATTGCATAAAAATGGATGACATAAACTTCCCATACGCTGTGGGGGGAACAGaaaa                                                                                                                                                                                                                                     >  2:652878/5‑250 (MQ=255)       aTTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGaaaa                                                                                                                                                                                                                                     >  1:508217/1‑250 (MQ=255)             taATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAAt                                                                                                                                                                                                                               <  1:454911/250‑1 (MQ=255)                                cTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCACAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCAc                                                                                                                                                                                                           >  1:143693/1‑251 (MQ=255)                                   tgGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTAt                                                                                                                                                                                                        >  1:412589/1‑251 (MQ=255)                                             aaaTGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGa                                                                                                                                                                                               >  2:375927/1‑250 (MQ=255)                                                                     ccATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTc                                                                                                                                                                      >  1:620413/1‑251 (MQ=255)                                                                     ccATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTc                                                                                                                                                                      <  1:233908/251‑1 (MQ=255)                                                                      cATTACTTTCAACAACCAGATTATCCTTGCAGGAATCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGCATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCt                                                                                                                                                                     >  1:148872/1‑251 (MQ=255)                                                                       aTTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCt                                                                                                                                                                     >  2:222727/1‑250 (MQ=255)                                                                                aacaacCAGATTATCCTTGCAGGACTCCTGGAAAATAAAACTGGTATGCCCATTCCCCCAAAAACGGCTGACAAAGCTTTCCAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATACACTTGACATAAGCTGTGGGAGACACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCg                                                                                                                                                            >  1:442160/1‑250 (MQ=255)                                                                                   aacCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACGATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATaa                                                                                                                                                        <  1:103576/251‑1 (MQ=255)                                                                                                     ggACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGc                                                                                                                                      <  2:636961/251‑1 (MQ=255)                                                                                                            gggAAAATACAACTGGAGTGCCCGTTCCCCAAACACCAGCAGACAAAGCGATCAGACTGAAGGAGGAACATATCACCCACATCGTTTTTAAACAATTGCACAAAAGAGAGTGACATAAACGTGCCATAAGCCGAGGGAGGCACAGAAAAATAAATGATACGGACTGTGTTCTCACTATGACCACTGAAATCAATATTGGGAGTCCCGGTTTTTAGCAATTCGATAAAATGCATATGTTTACTGCCGAtggt                                                                                                                               >  2:330286/2‑248 (MQ=255)                                                                                                                                   aTTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGAt                                                                                                                                  >  2:473816/1‑225 (MQ=255)                                                                                                                                   aTTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGAt                                                                                                                                  <  1:473816/225‑1 (MQ=255)                                                                                                                                                            cTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGttatt                                                                               >  1:8836/1‑251 (MQ=255)                                                                                                                                                                 aaaCTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATg                                                                           >  2:150810/1‑250 (MQ=255)                                                                                                                                                                      gACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACgg                                                                       <  1:335607/249‑1 (MQ=255)                                                                                                                                                                      gACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGtt                                                                     >  2:147852/1‑251 (MQ=255)                                                                                                                                                                                aatttcacgtacatacttcttAACTCATTGCAAAAAAATCGATTACACAAACTTGCCATAAGCCGTGGGAGGCTAAGAAAAACTAAAGACACGGAAAGTGATCTCACTTTGACCACTGAAATAAATATTCGGCGACACTGTTCCTAGTAATTCGATAAATTGAATATGTTTACTGCCGACAGATTTTAAAACATGCTAAACCTAAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTaa                                                           <  1:330286/232‑1 (MQ=255)                                                                                                                                                                                      ccccACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTa                                                     <  2:364504/250‑1 (MQ=255)                                                                                                                                                                                              cGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGatat                                             <  2:203860/251‑1 (MQ=255)                                                                                                                                                                                                  tttAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGAt                                          <  1:375927/250‑1 (MQ=255)                                                                                                                                                                                                      aCCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCACGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTc                                       >  2:462184/1‑249 (MQ=255)                                                                                                                                                                                                         cATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGTGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCTATAAATTGCATATGTTTACTGCCGAGAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTAGTAAACTCTACTTGATATGATTTCAGACg                                  <  1:615811/251‑1 (MQ=255)                                                                                                                                                                                                          aTTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCg                                                                                                                                                            <  2:448921/128‑1 (MQ=255)                                                                                                                                                                                                          aTTGCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCg                                                                                                                                                            >  1:448921/1‑128 (MQ=255)                                                                                                                                                                                                             gCATAAACATGGGTGAAATAGACTTGCCATAAGCTGTGGGAGGAACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTAGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTAAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTAAGACGCTCg                              >  2:409525/1‑251 (MQ=255)                                                                                                                                                                                                             gCATAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTc                               <  2:10371/250‑1 (MQ=255)                                                                                                                                                                                                               aTAAAAATGGTTGACATAAACTTGCCATAAGCTGTTGGAGGCACAGAAAAATCAATGAGACGGAATGTGTTCTCACTATGACCACTGAAATCTATATTCGGCGTCCCGGTTTCTAGCAATGCGATAAAGTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAAGGTGATTAATGACGGTTAACGTATTAAACTCTACTTGATATTATTTCAGACGCTCGAt                            <  1:632354/251‑1 (MQ=255)                                                                                                                                                                                                               aTAAAAATGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGAt                            >  1:510667/1‑251 (MQ=255)                                                                                                                                                                                                                 aaaaaTGGTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATg                           <  2:513362/250‑1 (MQ=255)                                                                                                                                                                                                                        gTTGACATAAACTTGCCATAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTAt                   <  2:17180/251‑1 (MQ=255)                                                                                                                                                                                                                                         aTAAGCTGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTAc  <  2:559949/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     GTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTTACTGCCGATAGATTTTAAAAAATGATAAACCTCAAGAGGGTAATGGACATTAACGTTATTAATGACGGTTAACGTATTAAACTCTACTTGATATGATTTCAGACGCTCGATGGCTGCTATCACTTTTGCAAAAGTAC  >  NC_000913/1579905‑1580387

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