Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,173,448	C→T	intergenic (‑575/+834)	gatD ← / ← gatB	galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,173,448	0	C	T	97.0%	99.1 / ‑6.8	33	intergenic (‑575/+834)	gatD/gatB	galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (19/13);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATG  >  NC_000913/2173204‑2173688                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  aTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGc                                                                                                                                                                                                                                           >  1:564963/1‑251 (MQ=255)       gTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGACTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCATGACCGCCGAACTTCGCATGTAAACGACTACGAGCTTGCTTAGCGATggg                                                                                                                                                                                                                                     >  1:247780/1‑251 (MQ=255)            gAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATggg                                                                                                                                                                                                                                     >  1:25751/1‑246 (MQ=255)                    ttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGACGATACCACC‑‑GCCGTATGTCACCGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAACGATGGGCGTTAAAccaccc                                                                                                                                                                                                                        >  2:102906/1‑248 (MQ=255)                              cACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGCAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt                                                                                                                                                                                                              >  1:503102/1‑251 (MQ=255)                               cgggcccccGCCATCGCCACGAAGAAGCCGATGGGGGCAAGACCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                                                                                                                                                                                                             <  1:148879/244‑1 (MQ=255)                                                           cGATGGGGGCAAGACCGCCAAACGGCAGCCCCTGATTACCCGGCACAAAAACAGAAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACt                                                                                                                                                                                <  1:284379/251‑1 (MQ=255)                                                            gATGGTGGCAAGACCGCAAAACGGCAGCCCCTGATTCCCCGGCACACAAACAGCAATAAAAATGGTGAGTGGGATAAAAATCAGGCTGGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTg                                                                                                                                                                               <  1:79156/251‑1 (MQ=255)                                                               ggtggCAAGATCGCCAAAAGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGACGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACAATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTgccg                                                                                                                                                                            >  1:383737/1‑251 (MQ=255)                                                                                   gcagcaCCTGATTACCCGGTACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTgg                                                                                                                                                                                                >  1:347760/1‑211 (MQ=255)                                                                                         ccTGATTCCCCGGCACACAAAAAGCAATAAAAAGGGTGTGTGGGATAAAAATCGGGCTGCCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATaca                                                                                                                                                  <  1:102906/251‑1 (MQ=255)                                                                                          cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGCTGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATacac                                                                                                                                                 >  1:197814/1‑251 (MQ=255)                                                                                           tGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATCCAcc                                                                                                                                                >  1:372447/1‑251 (MQ=255)                                                                                                        cacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTa                                                                                                                                    >  1:4419/1‑250 (MQ=255)                                                                                                                gCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCCCTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTTACATCGAAACCcgcg                                                                                                                             >  2:442356/1‑249 (MQ=255)                                                                                                                      aaaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc                                                                                                                     >  1:376914/1‑251 (MQ=255)                                                                                                                                            tCAGGCATGCCGATACCACC‑‑GCCGAATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAACGCCGAACTTCGCCTGTAAACGAATACGAGCTTGCTTAGCGATGGGCGTTAAACCATCAATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCCGCTGCAATAAACCTTTGACATCGTAACCCGCGAGGATGACGATATTCAGCCCCATCCCaaa                                                                                                >  2:511340/1‑250 (MQ=255)                                                                                                                                                            ccccccccCGTCTGCCCCGGCAGCACGGCCGGAACAAGCCCAATCGGGACCTCCGTACCGCCGACCTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCCCTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAgg                                                                                 <  2:230914/249‑1 (MQ=255)                                                                                                                                                              caccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGcc                                                                               >  1:183599/3‑250 (MQ=255)                                                                                                                                                               cccGCCGAATGTCCCAGCAGCAACGCCTGATCAAGGCCAATCAGGAACCCCTGACCGCTGACCTTCGCTTGTAAACGACTGCGAGCTTGCTCAGCGATAGGCATTACACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTACCAGCCTCACTTCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCCCGGTGACAGGCTc                                                                              <  1:2041/251‑1 (MQ=255)                                                                                                                                                                       aTGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGCGACAGGCTCGCCAAATg                                                                      >  1:216769/1‑251 (MQ=255)                                                                                                                                                                                              cAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATAtc                                                <  1:290819/250‑1 (MQ=255)                                                                                                                                                                                              cAAGACCAATCAGGAACTACTGACCGCCGAACTTCGCCTATAAACGACTACGAGCTTGCTTAGCAATGGGCGTAAAACCATCCATGATGGGTTGAA‑TCACCCGTGGCATTATCAGCCTCACTGCCGCCGTTTTTACCGCCAGCTGCAATACAACTTTGAAATCGTACACCGCGAGGATGCAGATATTCAGCCCCATCACAAACCCCACGGTGACAGGAGCGCCAACTGGACCACATATGAGCTAAATAtc                                                >  2:302921/1‑250 (MQ=255)                                                                                                                                                                                               aaGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTCAACCATCCATGCTGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCCGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCACCATCACAAAACCCACGGGGACAGGCACGCCAATTGGACCAACTTTGCGCTGAAAAtcgt                                              >  2:228955/1‑251 (MQ=255)                                                                                                                                                                                               aaGGCCAATCAGAAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGTCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATAtcgt                                              <  1:228955/251‑1 (MQ=255)                                                                                                                                                                                                aGGCCAATCAGGACCTCTTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTAAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATAtcgtc                                             <  1:524790/251‑1 (MQ=255)                                                                                                                                                                                                 gtcaaaccaggaatcctgccccgcagAACTTCGCCTGTAAACGACTAAGAGCTTGCTTAGCGAGGGGCGTAAAACCATCCATGA‑GGGTTTAAATCACCCGTGGAATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATAtcgtcg                                            <  2:564963/231‑1 (MQ=255)                                                                                                                                                                                                                   gtcggccgAACGTCGCCTGGAAACTACTTCGAGCTTGCTCAGCGAGGGGAGTTAAACCATCCATGAGGGGTTGAA‑TCACCCGGGGCTTTAGCCTCATCACTGCCGACGTTTTTACGGCCAGCTGCAATACACCTATGACAGCGTAACGCGCGAGGATGCCGATAAGCAGCTCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAGATTGGCGCTTAGTATCGTTGGCGCTAAATTTAATTc                            <  2:545327/245‑1 (MQ=255)                                                                                                                                                                                                                         cgactttcgcttgtaaatgcttACGAGCTTGCTTATCGTGGGGCGTTAACCCATCCATGATGGGTTTAA‑TCCCCCGTGGCATTAGCAGCATCACGGCCGCCGTTTTTACCGCCAGCTGCAATCCACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATGCGGTTAACg                    <  2:372447/230‑1 (MQ=255)                                                                                                                                                                                                                                 ccctggaaacgcctcggAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTGTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAAGTCGGTTAACGCCTGGGAt            <  2:24541/235‑1 (MQ=255)                                                                                                                                                                                                                                        aaCGACTACGAGCCTGCTTATCGATGGGCGTTAAACCTGACAGGATGGGCTTAA‑ACACCCGGGGCTTTAGCAGCATCACCGCCTCCGTTTTTCCCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGTATGCCGATAAACAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGACTATCGTCGGCGCTAAATTTACTTCGGTTAACGCATGGGATTTTCTCg     <  2:445324/251‑1 (MQ=255)                                                                                                                                                                                                                                           gaataagAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCCCCCGTGGCATTAGCAGCATCACTGACGCCGTTTTTACAGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGCGGATGCCGATAATCAGCCCCATCACAAACCCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCACTGAATATCTTCGGCGCTCAATTTAATTCGGTTAACGCCTGGGATTTTCTAgatt  >  2:146112/7‑250 (MQ=255)                                                                                                                                                                                                                                           gACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCgatg  >  1:389318/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAA‑TCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATG  >  NC_000913/2173204‑2173688

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