Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	78,406	A→T	G262G (GGA→GGT)	setA →	sugar exporter SetA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	78,406	0	A	T	93.1%	81.5 / ‑3.6	30	G262G (GGA→GGT)	setA	sugar exporter SetA
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (4/23);  total (6/24)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTCTGGTATTGATTGCATTTATGCTTCCGTCTGTGGCGCGGGTAGAACTGCCGTCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGACTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTCTTTAATAGCCGTATGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATT  >  NC_000913/78169‑78651                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      cAGTCTGGTATTGATTGCATTTATGCTTCCGTCTGTGGCGCGGGTAGAACTGCCGTCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAg                                                                                                                                                                                                                                          <  2:55164/251‑1 (MQ=255)   gTCTGGTATTGATTGCATTTATGCTTCCGTCTGTGGCGCGGGTAGAACTGCCGTCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCa                                                                                                                                                                                                                                        <  1:364908/251‑1 (MQ=255)     cTGGTATTGATTGCATTTATGCTTCCGTCTGTGGCGCGGGTAGAACTGCCGTCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACAGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGGTTCCTTGTGGGGACGGGAGctggtctgggaaaaacagaca                                                                                                                                                                                                                                       >  2:169960/1‑238 (MQ=255)      gggTATTGATTGCATTTATGCTTCCGGCTGTGGCGCGGGTAGAACTGCCGTCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATg                                                                                                                                                                                                                                     <  2:379758/250‑1 (MQ=255)        gTATTGATTGCATTTATGCTTCCGTCTGTGGCGCGGGTAGAACTGCCGTCGGAAAATACTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGCCCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGGTTGCTGATGGGGACGGCAACTGGACAGGAAAAACCAGGAATgg                                                                                                                                                                                                                                    >  1:35414/1‑249 (MQ=255)                 gCATTTATGCTTCCGTCTGTGGCGCGGGTAGAACTGCCGTCGGAAAATGCTTGATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTAGTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTctggctgg                                                                                                                                                                                                                          <  1:311966/251‑1 (MQ=255)                      tATGCTTCCGTCTGTGGCGCGGGTAGAACTGCCGTCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACAGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGctac                                                                                                                                                                                                                      <  1:169960/250‑1 (MQ=255)                            tCCGTCTGGGGCGCGGGTAGAACTGCCGGCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTc                                                                                                                                                                                                                <  2:201904/250‑1 (MQ=255)                                           ggTAGAACTGGCGGCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGAGCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGTCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAgc                                                                                                                                                                                                <  2:366487/251‑1 (MQ=255)                                               gAACTGCCGTCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAgcggcg                                                                                                                                                                                            <  1:80681/251‑1 (MQ=255)                                                aaCTGCCGTCGGAAAATGCCTGATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGa                                                                                                                                                                                           <  1:10021/251‑1 (MQ=255)                                                             aaaTGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAg                                                                                                                                                                              <  1:367270/251‑1 (MQ=255)                                                             aaaTGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAg                                                                                                                                                                              <  1:94293/251‑1 (MQ=255)                                                                                      tggcAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGGTGGGGACGGCAGCCGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGCTGGGCATAGCAGTGGCGGCAGGAGGACTGttttt                                                                                                                                                     >  1:157678/1‑250 (MQ=255)                                                                                       ggcAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACTGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTa                                                                                                                                                     <  1:303989/250‑1 (MQ=255)                                                                                                       aCGGATGTTATTTGTCGCTGCGACGTTAATGTGGACCTGCAACTCCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTGTCCTGATGGGGACGGCAGCTGGTCGGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGAtttta                                                                                                                                     <  2:77247/250‑2 (MQ=255)                                                                                                       aCGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGAttttat                                                                                                                                    <  1:399207/251‑3 (MQ=255)                                                                                                               ttgttGTCGCCTCGCCGTTAATTTGTTCCTGCAACACCACGTACATTATGCATATTCCGTTGTGTATCAGTATCGATTTAGGATTGCCAGACACACTGGCTGGTGTCCTGATGGGGTCGGCTGCTGTTCTGGAAATTCCAGCACTGATTTTGTCTGGCTACTATGTCAAACGTTATGGTAAGCTGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTATTGTTTTACAcctgagttattttaattcatat                                                                                                                             <  2:346976/247‑21 (MQ=255)                                                                                                                      cGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAACTGGGCTGGGAAAACCAGCAATGATTCTGGCTGGCTACTATGTCCAACGTTATGGGAAGCGGCGCATGATGGTCATAGGCGGGGCGGCAGGGGTAATGTTTTACACCGGATTGATTTTATTTCATAGCCGTATgg                                                                                                                     >  1:242053/1‑251 (MQ=255)                                                                                                                              cGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTtgatg                                                                                                             <  1:287422/251‑1 (MQ=255)                                                                                                                                                                    atggcGTTGTGGATCAGTAGCGAGCTAGGAGTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGGTTGACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATTACGCTTCAATTTTTTAACGCTGGATTTATCGGCAttg                                                                        <  2:301340/246‑1 (MQ=255)                                                                                                                                                                      gCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCAttgtt                                                                      <  1:433860/250‑1 (MQ=255)                                                                                                                                                                               gATCAGTAGCGAGTTAGGAGTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGTCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTg                                                            <  2:29556/251‑1 (MQ=255)                                                                                                                                                                                   aGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGAt                                                        <  1:257666/251‑1 (MQ=255)                                                                                                                                                                                        cGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTAt                                                   <  1:435953/251‑1 (MQ=255)                                                                                                                                                                                                      cAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCCATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGGATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCCGTATTTATCGGGATTGTTGCGGGGATTGGGGTGCTATGGTTTCAGGAttt                                      >  1:221536/1‑250 (MQ=255)                                                                                                                                                                                                                  cGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAaga                         <  1:231498/251‑1 (MQ=255)                                                                                                                                                                                                                   gggTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTTGCGTTGGTGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTTTTTCGGGGTTTGGGATGCTATGGGTTTAGGATTTGATGCCTggtag                          >  2:221113/1‑246 (MQ=255)                                                                                                                                                                                                                     gTTTCCTGATGTGGACGTCAGCTTGACTGTAAATCCCAGCAATGATTCTGGCTGGCTACGAGGACAAACGTTATGGTAAGCGGCGAATGAGTGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTCTTTAATAGCCGTATGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGTCATTGTTGCGGGTATGGGGATGCTATGGTTTCAGGATTTAATTCCTGGAAGAGCg                      <  2:35414/251‑1 (MQ=255)                                                                                                                                                                                                                                   aCGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTaccac        <  1:414520/251‑1 (MQ=255)                                                                                                                                                                                                                                         gCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCGGCTACCACCTTAtt  >  1:77185/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      CAGTCTGGTATTGATTGCATTTATGCTTCCGTCTGTGGCGCGGGTAGAACTGCCGTCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGACTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTCTTTAATAGCCGTATGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATT  >  NC_000913/78169‑78651

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: ‑