Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 968,837 G→A D158N (GAT→AAT)  lpxK → tetraacyldisaccharide 4'‑kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913968,8370GA96.5% 91.0 / ‑3.3 30D158N (GAT→AAT) lpxKtetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (13/15);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01

GGTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGG  >  NC_000913/968591‑969075
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
ggTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAAtgt                                                                                                                                                                                                                                            <  1:298341/251‑1 (MQ=255)
 gTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAAtgtg                                                                                                                                                                                                                                           <  1:69266/251‑1 (MQ=255)
     tCGCGTGGGGGTCGTATCGCGTGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGaa                                                                                                                                                                                                                                        <  1:189835/250‑1 (MQ=255)
                         ggggATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGGTTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGGAATTGTTGTTTTTGATGGtgtg                                                                                                                                                                                                                   >  1:61796/1‑251 (MQ=255)
                         ggggATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTGTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCGGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTc                                                                                                                                                                                                                                  <  1:377462/236‑1 (MQ=255)
                         ggggATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGGTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGtgtg                                                                                                                                                                                                                   >  1:179179/1‑251 (MQ=255)
                         ggggATAAGGTGGTACGGCTGAATCTTATCCGCTGTTCTTGTCGGCAGATACCACAACAGCACAGGCGGGACATGCACCCGTTGTAATTTCTCAAAGCAATGACGCGGCTGCTGCGGTTTCTCCCGTTCGGTCTGGTGAGGAAAAAGCCATTTTGGCGCAACCCCCTGATGCGCAGATCATTGTAAACGCAGACGGTTTCCAGCATTACAGTCTGGGGCGGAATGTGGAAATTGTcagatagggagagcgg                                                                                                                                                                                                                   >  2:377462/1‑236 (MQ=255)
                                 ggtggtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTTTCAACCCACTGGTGCGCCCGTTGCGGTTTCTTCCGTTCGTTCTGGTGCGGTAAAAACCAATCTGGCGCAACACCCTGATGTGCAGAACAGCGTAACCGCCGGCGGTTTTCAACATTAACGGCTGGCGCGTAATGTGGAAATTGGCGTTATTGTTGGTGTGCGTCGCtt                                                                                                                                                                                                           >  1:226123/1‑251 (MQ=255)
                                                    aTCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGtt                                                                                                                                                                                         <  1:382625/250‑1 (MQ=255)
                                                    aTCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTACCAACGCACTGATGCGCCTGTTGCGGTTTCGCCCGTTCGTTCTGAGGCGGTAAAAGCCATGCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGGCTTGCGCGTAATTTGGAAATTGGCGTGATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGtt                                                                                                                                                                                         <  2:192383/250‑1 (MQ=255)
                                                          tGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTTTCAACGCACTGATGCGCCTGTTGCGGGTTCTTCCGTTCGTTCTGATGCGGGAAAAGCCCTTCTTGCGCAACCCCCTGATGTGCAGATCATCGTAACCGACGCCGGTTTACAGCATTACCGTCTTGCGCGTTATGTGGGAATTGTCGTTTTTGATGGTGTGCGTCGCTTTGGCCATGGCTGGTGGTTTCcggcg                                                                                                                                                                                  >  1:387706/1‑251 (MQ=255)
                                                          tGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAACCTGTGTTGATTTATCTACGCCATGTTGCGCCTGTTGAGGGTTTTCCCGATTGTTTTAATGCGGTAAAAACCATTCTTGCGCAAAACCCAGTTGTGCCGATCATCGTAACCGAATAAGGGTTACAGACTTACCGTCTGGCGCGTAAGGTGGGAATTGTTGTTCTTGATGGGGGGCGACGCTTTGGGAAAGGCTGGTGGtttgcggcg                                                                                                                                                                                  >  1:81943/1‑244 (MQ=255)
                                                                   cGGCAGATACCGCCACAGCACAGGCGGGTGATGATCCTGTTTTGATGTATCAACGCACTGATGCGCCTGTTGCTGTTCCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTTGCGCAACCCCCTGATGTGCAGATCATCGTAACCGATGACGGTTTACAGCATTACCGTCTGGCTCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAAt                                                                                                                                                                          <  1:27483/250‑1 (MQ=255)
                                                                       aGATACCACAACAGCACAGGCGGGAGGTGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGGTTTTTCCGTTCGTTCTGATGCGGTCAAAGCCATTCTTGCGCACCAACCCGGTGTGCAGATCATCGTAACCGACGGCGGTTTACAGCCTTACCGTCTGGCGCGTAATGGGGGACTTTTCGTTTTTGTTTGGGTGGGTCGCTTTGGCCATGGATGGGGGTTGCCCGCGGGGcccgaggggg                                                                                                                                                                     >  1:13667/1‑243 (MQ=255)
                                                                            ccACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTTATGCGCCTGTTGCGGGTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTgagcga                                                                                                                                                                <  2:350875/251‑1 (MQ=255)
                                                                                          gCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTaaa                                                                                                                                                  <  1:30911/251‑1 (MQ=255)
                                                                                            gggTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCCATGCGGGAGCGAGCGGGGCGCTTAAAGt                                                                                                                                                >  1:172824/1‑251 (MQ=255)
                                                                                                gatgaGCCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGTGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCgg                                                                                                                                             <  2:387706/250‑1 (MQ=255)
                                                                                                              gATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATc                                                                                                                               <  1:321370/250‑1 (MQ=255)
                                                                                                               aTTTAGCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACTCCCGGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAATGTCGGTTGATGCGGTAATCGt                                                                                                                             <  2:308832/251‑1 (MQ=255)
                                                                                                                        cGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGACAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGCGGTAATCTGCAAcggcgg                                                                                                                    >  2:413139/1‑251 (MQ=255)
                                                                                                                                       gTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGt                                                                                                                                                <  2:323080/208‑1 (MQ=255)
                                                                                                                                       gTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGt                                                                                                                                                >  1:323080/1‑208 (MQ=255)
                                                                                                                                                    ccGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTAGTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCTCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTTAAATCCCCATg                                                                                        <  2:21483/251‑1 (MQ=255)
                                                                                                                                                         cgttcTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCCGCGGGGGAATCCCCATGCATcc                                                                                   >  1:159004/1‑250 (MQ=255)
                                                                                                                                                          gttcTGCTGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATACTCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAATTGGGTGATGAGGGAATCGTCAACGGGGGAGTTCCGCGCTGCGGGGTAAACCCCATGcctatt                                                                                  >  1:269089/1‑246 (MQ=255)
                                                                                                                                                                     gTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGGCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCa                                                                       <  2:179179/251‑1 (MQ=255)
                                                                                                                                                                                     cgcAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGACATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCCATGCGTGGGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGTCCCTCGCAGCGATGAAATCACCAAGCACCTTCTGCCGGGGCcggcgggggattttcct                                                       >  1:342486/1‑240 (MQ=255)
                                                                                                                                                                                             ccTGATGTGCAGATCATCGGAACCGACGACGGTTTACAGCAGTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGGCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTTATGCGGTAATCGTCAGCGGCGGTGTCCCTCGCAGCGGTGAAATCCCCGTGCAGCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTAc                                               <  2:276306/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGGGAAATACCCATGCATCTTCTGCCGGGTCAGGCGGGGAATTTACGTACCGGTACGCGGTGTGACGGTTCTCAGCTTGAACCTTTAGTGGCGGTGGCggg  >  1:228191/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
GGTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGG  >  NC_000913/968591‑969075

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