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

GTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTG‑GCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGATTG  >  NC_000913/968592‑969080
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
gTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAAACTGTGTTGATTTATCAACGCCCTGAAGGGGCTGGTGCGGGTTCTCCCGTTCTTTCTTATTCGGGAAAAGCCCTTCTG‑GGGCAACACCCCGATGGGCAGATCATCGTAACCGACGACGGTGTACAGAATTACCCTCTGGCGCGTAAtgtg                                                                                                                                                                                                                                                >  2:288402/1‑251 (MQ=255)
        gTGGGTGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTG‑GCGCAACTCCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGt                                                                                                                                                                                                                                        <  1:189923/251‑1 (MQ=255)
        gTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCGCAACAGCACAGGCGGGGGATGAAACTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTG‑GCGCAACACCCCGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCAATACCGTCTGGCGCGGAATGTGGAAATTGt                                                                                                                                                                                                                                        >  1:353850/1‑251 (MQ=255)
            gggTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGTCGGTGGATGAACCTGTGTTGATTGATCAACGCACTGATGCGCCTTTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAATCCACTCTG‑GCGCAACTCCCTGAGGTTCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGtt                                                                                                                                                                                                                                    <  1:189687/251‑1 (MQ=255)
                        ggggATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGATGCTGAACCTGTGGTTATTTATTAACGAACTGATGCGCCTGTTGCGGGTTCTCCCGCTCGTTCTTATGCGGGAAAAACCCTT‑TGTGCGCCACACCCTCATGTGCAGATCATCGTTACCGACGGCGGTTTACAGCAATTCCGTCTGGGGCGTAATGTGGGCATTGTGGTTATTTATgggttg                                                                                                                                                                                                                        >  2:98682/1‑247 (MQ=255)
                          ggATATGGTGGTAAGGTTGAATCTTCTCCGCTTTTATTGTCGGCAGATACCACAACAGCGCAGGCGGGTGTTGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGGTCTGATGCGGTAAAATCCATTCTG‑GCGCAACCCCTTGATGTGCAGATCATCGTAACAGACGACGGTTTACATCATTACCGTCTGGCGCGGAATGTGGATTTTGTCGTTATTGATGGTGTGct                                                                                                                                                                                                                      <  2:283484/251‑2 (MQ=255)
                                         gCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGGACCTGTGTTGATTTATAAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGGTCTGATGCGGGAAAAGCCATTCTG‑GCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACCGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGGTATTGATGGTGTGGGTCGCTTTTGCAATgg                                                                                                                                                                                                       >  2:7818/1‑251 (MQ=255)
                                                   aTCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGGTTTATCAACGCACTGATGCGCCTGTTGCGGGTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTG‑GCGCAACACCCTGATGTGCCGATCATCGTAACCGAAGCCGGTTTACAGCATTACCGTCTGGCGCGTAATGGGGGAATTGGAGTTTTTGATGGGGTGCGTCGCTTTGGCAATGGcgggggggtg                                                                                                                                                                                             >  2:153260/1‑242 (MQ=255)
                                                    tCCGCTGTTATTGCCGGCAGATACCACAACAGCACAGGCGGGTGATGGAACTGTGTTGATTTATCAACACCCCGACCCGCCTGTTCCGGTTTCTCCCGTTCGTTCTGATGCGGCAAACACCATTCTT‑GCGCAAAACCCTCATGTGCAGATCCTCGTAACCCACGCAGGGTTACAGCATTCACGTCTGGCCCGTCATGTGGGAATTTTCGTTATTGATGGTGAGCGTCGCTTTCGCAAAGGCTGGGGGTTGc                                                                                                                                                                                            >  1:197857/1‑251 (MQ=255)
                                                       gCTGTTATTGTCGGCAGATTCCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTG‑GCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGcctg                                                                                                                                                                                         <  2:34720/251‑3 (MQ=255)
                                                                 tCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAAAGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTG‑GCGCCACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTTTCGTTTTTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCgggggcca                                                                                                                                                                                >  1:2064/1‑246 (MQ=255)
                                                                 ccGGCAGATACCACAACAGTGCAGGCGGGTGCGGAACCTGTGTTGATTTATCAACGCACTGAGGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGATAAAAGCCATTCTG‑GCTCAGCACCCTGATGTGCAGGTCATCGTAACCGACGACGGTTTACAGCACTACCGTCGGGCGCGTAATGGGGAAATTGTCGTTATTGATGGAGTGCGTCGCTTTGGCAATGGTTGGGGGTTGCCGGCGTGGCCAAt                                                                                                                                                                               <  2:148506/250‑1 (MQ=255)
                                                                       gattccccaacaGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCGG‑GCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAAGGTGGAAATTGTCGTTATTGATTGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTga                                                                                                                                                                         <  2:245785/244‑1 (MQ=255)
                                                                                      caGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGCTGCGGTAAAAGCCATTCTG‑GCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGGGAGCGGGCGGGGCGCtt                                                                                                                                                          >  1:34459/1‑251 (MQ=255)
                                                                                            ggTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTG‑GCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTc                                                                                                                                                    <  2:58866/251‑1 (MQ=255)
                                                                                                            ggattttgcAACGCACTGATGCGCCTGTTGTGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTG‑GCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTGTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTGGATGCGGTAATc                                                                                                                                    <  2:308862/243‑1 (MQ=255)
                                                                                                                                tctccggTTGCTGTTTCTCCCGTGCGTTCTGATGCGGTAAAAGCCATTCTG‑GCGCAACTCCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTGTATGCGGTAATCGTCAACGGCGGTGTCCCTc                                                                                                                 <  2:39065/244‑1 (MQ=255)
                                                                                                                                                  cccGTTCGTGCGGATGCTGTAAAAGCCATTCTG‑GCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCAt                                                                                              <  2:2064/251‑1 (MQ=255)
                                                                                                                                                                                                           cagcGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTGATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGtt                                       <  1:387304/247‑1 (MQ=255)
                                                                                                                                                                                                                                               ggCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTTCCTCGCCGCGGTGAAATTCCCATGCACCTGCTGGCGGGTCAACGGGTGGATTTAAGTACCGGGACGCGTTGGGACGTTTCTCAGCTTGAACATGGAGTgggggtggggggggttg  >  1:251371/1‑236 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
GTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTG‑GCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGATTG  >  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: ‑