Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 4,889,556 G→A intergenic (+2/‑151) hyi → / → glxR hydroxypyruvate isomerase/tartronate semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,889,5560GA88.5% 67.4 / ‑3.1 26intergenic (+2/‑151)hyi/glxRhydroxypyruvate isomerase/tartronate semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (14/9);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGG‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGA‑TATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGATCGCTGCCGGTGCCGTGGCC  >  AE015451/4889269‑4889847
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                      
cGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGTCTGAAAACCCACAACCCAACCTGAAA‑A‑ACAACCCGGGAc                                                                                                                                                                                                                                                                                         >  2:171854/1‑301 (MQ=255)
         cAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGTGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGGTGAAGACCCCCTACTCACTCT‑AGGAA‑ACAACTCGGGACCCTGTGGGa                                                                                                                                                                                                                                                                                >  1:443353/1‑301 (MQ=255)
         cAGCGACAACCTGTTCCGGCAATACGAGATCTACAACATGCAGATCATGGAAGGTGACCTGGCTCGCACAATGGAAGCGAAACTGAAGATGATCAACCAGAACCCGCTGGCCGAGAACCCAGGCCGAAATGAACCAGGCACCGGCGAGTTCAAATACCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGa                                                                                                                                                                                                                                                                                <  1:32776/301‑1 (MQ=255)
                              aTACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACAcc                                                                                                                                                                                                                                                           <  1:171928/301‑1 (MQ=255)
                                                           aaggtaaactggacggtaCCATGGAAGCCAACCGGAATCTGATCAACCACCGCCAGCTGGACGACAACACAGGCCGCAATGAACCAGGCACCGGCGAGATCAACGACCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTGCCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGAGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTCAGA‑A‑ACAACTCAGGACCCTGTGGGAGCGGGCTCGCCCGCGAACAGCGGCAAAGCCGGTGCCATCCACCGAGTTGt                                                                                                                                                                                                                              <  2:52368/284‑1 (MQ=255)
                                                            aGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑A‑ACAACTCTGGACCCTGCGGGAGCGGGCTCGCCCGCGAACACCGGCAAAACCCGTGCCCTCCACCGAGTTGTc                                                                                                                                                                                                                             >  1:37306/1‑301 (MQ=255)
                                                                                                      cAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCCCCCACAGGGGGCCAAAgcttg                                                                                                                                                                                   >  1:335335/1‑301 (MQ=255)
                                                                                                        aCCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCCCCGAAGCGGGCCTGGGCTGGCTGAAAACCCACACCGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAACGCCGGTCCCATCCACCGCGTTGTCTGCTTCGCGCGCGAGCACGCCCCCACAGGTGCCCaaaccctgct                                                                                                                                                                                 >  1:23917/1‑294 (MQ=255)
                                                                                                            ctgacagcgagCCGACCACCCAGGCCGCAATGAACCAGGCGCGGGAGAGATCAACTACCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTgag                                                                                                                                                                             <  2:116573/291‑1 (MQ=255)
                                                                                                             aTCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGCCAAAGCCGGCGCCATCCACCGAGTTGCCTGCTTCGCGGGCGAGCCCGCTCCCATAGGGGGCCAAAGCTTGATTgagg                                                                                                                                                                            >  1:703161/1‑300 (MQ=255)
                                                                                                                       acGACAACCCAGGCCGCAATGAACCAGGCAACGGCGAGATCAACTACCGCTTCC‑GTGTAGAACACGTGGACCGCATTGGCTACCAGGGCTGGGTGGGGGCGGAGAACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGA‑TATCACAAtt                                                                                                                                                                 <  1:303274/300‑1 (MQ=255)
                                                                                                                                   gccgacaTGACCCAGGCACCGGCGAGATCAACTACCGCTTCCGGT‑TCGAACACCTGGACCGCATGGGCGACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCCCCGAAGCGGGCCTGGGCTGGCTGACAACCCACAACGCAATCTAAGG‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGA‑TATCACGATTACAAAGAGGTaa                                                                                                                                                     <  2:702861/295‑1 (MQ=255)
                                                                                                                                      gCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTG‑GACTATCACACTTACAAAGAGGTAAtgt                                                                                                                                                  >  2:21744/1‑299 (MQ=255)
                                                                                                                                               caggcccgcgcgggaTCAACTACCGCTGCCTGT‑TCGAAACCCTGGACCCCATTGGCTACCGGGGCTGTGTGGGCGCGGGTTACAGGCCGCTGACCGCCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑C‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGA‑TATCACAATTACAAAGAGGTAATTTCTCATGGCt                                                                                                                                         <  2:599027/288‑1 (MQ=255)
                                                                                                                                                                aaCTACCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGA‑TATCACAATTAAAAAGAGGAAATTTCACATGGCTAAAATCgggttactcgg                                                                                                                        >  1:342792/1‑292 (MQ=255)
                                                                                                                                                                  cTACCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑A‑ACAACTCGGGCCCCTGTGGGAGCGCGCTCGCCCGCGACAACCGGCAAAGCCGGTGCCATCCACCGAGTTCTCTGCTTCGCGGTCGAGCCAGCTCCCACAGGGGTCCAAAGCTTGCTTGCGA‑TATCACAATTACAAAGAGGTAATTTCTCATGGCTTACATCGGTTTCCTCGGCa                                                                                                                      >  2:542302/1‑301 (MQ=255)
                                                                                                                                                                    aCCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGA‑TATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCTGCAcc                                                                                                                    >  2:379106/1‑301 (MQ=255)
                                                                                                                                                                          tCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGA‑TATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatc                                                                                                              >  1:603612/1‑301 (MQ=255)
                                                                                                                                                                               t‑tCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGA‑TATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATggg                                                                                                         <  2:739499/301‑1 (MQ=255)
                                                                                                                                                                                                                                               accaccGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGA‑TATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACccaccaaga                                          >  1:543387/1‑298 (MQ=255)
                                                                                                                                                                                                                                                ccaccGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGA‑TATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCacgac                                         <  1:313175/301‑1 (MQ=255)
                                                                                                                                                                                                                                                         ggggCCGGGGCTGGCTGCAGACCCACAACGGAATC‑GAAG‑AGACAGCTCGGGACCCTGTGGGAGCGGGCTCGCCCGGGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGA‑TATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCa                                <  1:542532/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                           aaacacacgaagCAATCTGGGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCCGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCACGCTCCCACAGGGGGCCAAAGCTTGCTTGAGA‑TATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGATCGCTGcc              <  2:23902/290‑1 (MQ=255)
                                                                                                                                                                                                                                                                             aCCCACAACGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGA‑TATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTCCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCGGGTCACAGCCTGTTCATTTCCACCAACCACGACGCCGCGCCAGCCGACCAGCTCGCTGCCgg            >  2:607694/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                              cccACAACGCAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCTACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGTCAAAGCTTGCTTGAGA‑TATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCCCAGAACCTGCAAAAGGCAGGTCAAAGCATGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGATCTGACCGCTGCCggc           >  1:578732/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                       cAATCTGAGA‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGA‑TATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGATCGCTGCCGGTGCCGTGGcc  >  1:379873/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                      
CGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAACCACATCCAGCTGGCCGACAACCCAGGCCGCAATGAACCAGGCACCGGCGAGATCAACTACCGCTTCCTGT‑TCGAACACCTGGACCGCATTGGCTACCAGGGCTGGGTGGGCGCGGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGG‑A‑ACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGA‑TATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAGAACCTGCAAAAGGCAGGTCACAGCCTGTTCATTTCCACCCACCACGACGCCGCGCCAGCCGACCTGATCGCTGCCGGTGCCGTGGCC  >  AE015451/4889269‑4889847

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