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,8370GA95.0% 54.3 / ‑5.9 20D158N (GAT→AAT) lpxKtetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (8/11);  minor base C (1/0);  total (9/11)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGATTG  >  NC_000913/968592‑969080
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
gTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAAACTGTGTTGATTTATCAACGCCCTGAAGGGGCTGGTGCGGGTTCTCCCGTTCTTTCTTATTCGGGAAAAGCCCTTCTGGGGCAACACCCCGATGGGCAGATCATCGTAACCGACGACGGTGTACAGAATTACCCTCTGGCGCGTAAtgtg                                                                                                                                                                                                                                                >  2:288402/1‑251 (MQ=255)
        gTGGGTGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACTCCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGt                                                                                                                                                                                                                                        <  1:189923/251‑1 (MQ=255)
        gTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCGCAACAGCACAGGCGGGGGATGAAACTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCCGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCAATACCGTCTGGCGCGGAATGTGGAAATTGt                                                                                                                                                                                                                                        >  1:353850/1‑251 (MQ=255)
            gggTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGTCGGTGGATGAACCTGTGTTGATTGATCAACGCACTGATGCGCCTTTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAATCCACTCTGGCGCAACTCCCTGAGGTTCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGtt                                                                                                                                                                                                                                    <  1:189687/251‑1 (MQ=255)
                        ggggATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGATGCTGAACCTGTGGTTATTTATTAACGAACTGATGCGCCTGTTGCGGGTTCTCCCGCTCGTTCTTATGCGGGAAAAACCCTTTGTGCGCCACACCCTCATGTGCAGATCATCGTTACCGACGGCGGTTTACAGCAATTCCGTCTGGGGCGTAATGTGGGCATTGTGGTTATTTATgggttg                                                                                                                                                                                                                        >  2:98682/1‑247 (MQ=255)
                          ggATATGGTGGTAAGGTTGAATCTTCTCCGCTTTTATTGTCGGCAGATACCACAACAGCGCAGGCGGGTGTTGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGGTCTGATGCGGTAAAATCCATTCTGGCGCAACCCCTTGATGTGCAGATCATCGTAACAGACGACGGTTTACATCATTACCGTCTGGCGCGGAATGTGGATTTTGTCGTTATTGATGGTGTGct                                                                                                                                                                                                                      <  2:283484/251‑2 (MQ=255)
                                         gCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGGACCTGTGTTGATTTATAAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGGTCTGATGCGGGAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACCGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGGTATTGATGGTGTGGGTCGCTTTTGCAATgg                                                                                                                                                                                                       >  2:7818/1‑251 (MQ=255)
                                                   aTCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGGTTTATCAACGCACTGATGCGCCTGTTGCGGGTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCCGATCATCGTAACCGAAGCCGGTTTACAGCATTACCGTCTGGCGCGTAATGGGGGAATTGGAGTTTTTGATGGGGTGCGTCGCTTTGGCAATGGcgggggggtg                                                                                                                                                                                             >  2:153260/1‑242 (MQ=255)
                                                    tCCGCTGTTATTGCCGGCAGATACCACAACAGCACAGGCGGGTGATGGAACTGTGTTGATTTATCAACACCCCGACCCGCCTGTTCCGGTTTCTCCCGTTCGTTCTGATGCGGCAAACACCATTCTTGCGCAAAACCCTCATGTGCAGATCCTCGTAACCCACGCAGGGTTACAGCATTCACGTCTGGCCCGTCATGTGGGAATTTTCGTTATTGATGGTGAGCGTCGCTTTCGCAAAGGCTGGGGGTTGc                                                                                                                                                                                            >  1:197857/1‑251 (MQ=255)
                                                       gCTGTTATTGTCGGCAGATTCCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGcctg                                                                                                                                                                                         <  2:34720/251‑3 (MQ=255)
                                                                 tCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAAAGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCCACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTTTCGTTTTTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCgggggcca                                                                                                                                                                                >  1:2064/1‑246 (MQ=255)
                                                                 ccGGCAGATACCACAACAGTGCAGGCGGGTGCGGAACCTGTGTTGATTTATCAACGCACTGAGGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGATAAAAGCCATTCTGGCTCAGCACCCTGATGTGCAGGTCATCGTAACCGACGACGGTTTACAGCACTACCGTCGGGCGCGTAATGGGGAAATTGTCGTTATTGATGGAGTGCGTCGCTTTGGCAATGGTTGGGGGTTGCCGGCGTGGCCAAt                                                                                                                                                                               <  2:148506/250‑1 (MQ=255)
                                                                       gattccccaacaGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCGGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAAGGTGGAAATTGTCGTTATTGATTGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTga                                                                                                                                                                         <  2:245785/244‑1 (MQ=255)
                                                                                      caGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGCTGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGGGAGCGGGCGGGGCGCtt                                                                                                                                                          >  1:34459/1‑251 (MQ=255)
                                                                                            ggTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTc                                                                                                                                                    <  2:58866/251‑1 (MQ=255)
                                                                                                            ggattttgcAACGCACTGATGCGCCTGTTGTGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTGTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTGGATGCGGTAATc                                                                                                                                    <  2:308862/243‑1 (MQ=255)
                                                                                                                                tctccggTTGCTGTTTCTCCCGTGCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACTCCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTGTATGCGGTAATCGTCAACGGCGGTGTCCCTc                                                                                                                 <  2:39065/244‑1 (MQ=255)
                                                                                                                                                  cccGTTCGTGCGGATGCTGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCAt                                                                                              <  2:2064/251‑1 (MQ=255)
                                                                                                                                                                                                           cagcGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTGATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGtt                                       <  1:387304/247‑1 (MQ=255)
                                                                                                                                                                                                                                               ggCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTTCCTCGCCGCGGTGAAATTCCCATGCACCTGCTGGCGGGTCAACGGGTGGATTTAAGTACCGGGACGCGTTGGGACGTTTCTCAGCTTGAACATGGAGTgggggtggggggggttg  >  1:251371/1‑236 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
GTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGATTG  >  NC_000913/968592‑969080

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