Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,789,647 G→A G145G (GGC→GGT waaH ← UDP‑glucuronate:LPS(HepIII) glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,789,6470GA97.0% 93.9 / NA 33G145G (GGC→GGTwaaHUDP‑glucuronate:LPS(HepIII) glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (20/12);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAAC  >  NC_000913/3789415‑3789891
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
tCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCATTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACACCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCg                                                                                                                                                                                                                                     >  2:353798/1‑250 (MQ=255)
    tGACGTATTATGCCGATAATATTTATATAATGATTGCTCGGTATATCGCGCCCGCAGCGCATTAAACATGAATTCTGTTGTCCAGACAATATCCTGCTGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACCATACCATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGACAGCCCCATCCGCAGCCAGTCCGGGCCGGTTACTACACCGGTTGAGCGAAGGCGATCgc                                                                                                                                                                                                                                >  1:214358/1‑250 (MQ=255)
            tatGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCATTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGggatg                                                                                                                                                                                                                        >  1:129354/1‑251 (MQ=255)
               gCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCATTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGggatggat                                                                                                                                                                                                                     >  2:368637/1‑251 (MQ=255)
                    taataTTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCATTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACAGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGcc                                                                                                                                                                                                                 >  1:314761/1‑250 (MQ=255)
                             tataATGATTGCTCGGTATATCGAGCAAGCAGCGCCTTAAACCTGAATTCTGTTCTCCAGACAATATCCTGCTGATGTAATACGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCAAACGGCGCGAAGAAAGCCCCCACCACAGCCAGTCCGGGCAGGTTAATACACCGGATGAGCGCAACCGCTCGCTGTGGATGGCTTGCCAGGTTtcccc                                                                                                                                                                                                       >  2:382521/1‑248 (MQ=255)
                               taATGATTGCTCGGTATATCGCGCACGCAGCGCATTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCg                                                                                                                                                                                                     >  1:270491/1‑251 (MQ=255)
                                                                  tAAACATGAAGTCTGTTGTCCAGACAATAGCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAAGAACATCGCGGCGATAAACCCCCATCCAGCCAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACAGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTgcgc                                                                                                                                                                   <  1:403628/250‑1 (MQ=255)
                                                                        tGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCg                                                                                                                                                            >  1:118638/1‑251 (MQ=255)
                                                                                                ccTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGtc                                                                                                                                    <  1:290308/251‑1 (MQ=255)
                                                                                                      gatgTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGc                                                                                                                              >  2:312700/1‑251 (MQ=255)
                                                                                                                      aaaTTTAATGTTATTTTTCACAATAACAGCGCGGCGATAAACCCCCAGCCAGACAACGTGATGCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGTATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCTCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTGTCGTACAtggtgt                                                                                                              <  2:287523/251‑2 (MQ=255)
                                                                                                                            aaTGTTATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGa                                                                                                        <  1:353798/251‑1 (MQ=255)
                                                                                                                                ttATTTTTAACAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTc                                                                                                    >  2:321699/1‑251 (MQ=255)
                                                                                                                                          cAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGc                                                                                           >  2:330027/1‑250 (MQ=255)
                                                                                                                                          cAATAACATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGATAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGc                                                                                           >  1:143355/1‑250 (MQ=255)
                                                                                                                                              aaCATCGCGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCg                                                                                       >  1:373896/1‑250 (MQ=255)
                                                                                                                                                 atcgcGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACa                                                                                    >  1:224724/5‑250 (MQ=255)
                                                                                                                                                        cGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGAc                                                                            <  1:28225/251‑1 (MQ=255)
                                                                                                                                                         gATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGAc                                                                            >  2:365559/1‑250 (MQ=255)
                                                                                                                                                             aaCCCCCATCCAGACAACGGGAGTCCAACGGCTCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGTTCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTGCACGAAAACACCAGGCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTGACGCTATGGTCATCAGCGTTTCGTGCGGGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATAtt                                                                       <  1:227458/251‑1 (MQ=255)
                                                                                                                                                                  ccATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGccc                                                                  <  2:118638/251‑1 (MQ=255)
                                                                                                                                                                      ccAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAACAGCGACATATTTGCCCGTTg                                                              >  1:73411/1‑251 (MQ=255)
                                                                                                                                                                                   gTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGGCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATccc                                                 <  2:380229/251‑1 (MQ=255)
                                                                                                                                                                                                aagaaAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGcc                                                            >  2:137534/1‑227 (MQ=255)
                                                                                                                                                                                                aagaaAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGcc                                                            <  1:137534/227‑1 (MQ=255)
                                                                                                                                                                                                  gaaAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACat                                                                          >  2:203213/1‑211 (MQ=255)
                                                                                                                                                                                                  gaaAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACat                                                                          <  1:203213/211‑1 (MQ=255)
                                                                                                                                                                                                        cccATCCGCAGCCAGTCCGGGCCGGTTAATACAACGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAACGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCCTCGACAAATGCGCAATATTTGCCCGTTGCCACTTCAATCCCACGATTAAGCGCCACAGATgc                            >  2:248013/1‑251 (MQ=255)
                                                                                                                                                                                                                        ccGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTg             <  1:430624/250‑1 (MQ=255)
                                                                                                                                                                                                                         cGGGCCGGTTAATACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGAt           <  2:129354/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ttAATACACCGGTTGAGCGCAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGcaacaa    >  2:43393/1‑250 (MQ=255)
                                                                                                                                                                                                                                   aaTACACCGGTTGAGCGAAGGCGATCGCTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCTATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAaac  <  1:312700/250‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
TCACTGACGTATTATGCAGATAATATTTATATAATGATTGCTCGGTATATCGCGCACGCAGCGCGTTAAACATGAATTCTGTTGTCCAGACAATATCCTGATGATGTAATCCGGCAATAAATTTAATGTTATTTTTAACAATAACATCACGGCGATAAACCCCCATCCAGACAACGTGAGTCCAACGGCGCGAAGAAAGCCCCATCCGCAGCCAGTCCGGGCCGGTTAATACGCCGGTTGAGCGAAGGCGATCGGTGGGGATGGATTGCCAGGTTTCTCCCGTTTCACGAAAACACCAGTCAGCGTTGCACTGCGCCACGTCGAGGTCGTCCTCTAACGCCATGGTCATCAGCGTTTCGTACATGGTGGGATAGACTTCATCGTCAGCATCGACAAAAGCGACATATTTGCCCGTTGCCACTTCAATCCCACGATTACGCGCCACCGATGCGCCAGCATTCGCCTGATGCAACAAAC  >  NC_000913/3789415‑3789891

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