Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 969,352 G→A *329* (TAG→TAA lpxK → lipid A 4'kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913969,3520GA100.0% 83.9 / NA 27*329* (TAG→TAAlpxKlipid A 4'kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/14);  total (13/14)

GACTGAAAAAGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGCGTCGTT  >  NC_000913/969217‑969480
                                                                                                                                       |                                                                                                                                
gACTGAAAAAGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTa                                                                                                                               <  2:24994/139‑1 (MQ=255)
          gATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCa                                                                                                                     <  2:868157/139‑1 (MQ=255)
                       gCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGAt                                                                                                        <  1:105281/139‑1 (MQ=255)
                        ccGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATg                                                                                                       <  2:315410/139‑1 (MQ=255)
                                                  ggtATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACct                                                                             >  1:487167/1‑139 (MQ=255)
                                                     aTTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTGATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCtctc                                                                          >  2:540157/1‑139 (MQ=255)
                                                              tAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGa                                                                 >  2:229950/1‑139 (MQ=255)
                                                                    cacaGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGcgcg                                                           <  1:229950/139‑1 (MQ=255)
                                                                      caGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTa                                                         <  2:121680/139‑1 (MQ=255)
                                                                                   atgaACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGc                                            >  1:491082/1‑139 (MQ=255)
                                                                                   atgaACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGc                                            >  1:286989/1‑139 (MQ=255)
                                                                                          aGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGc                                         <  1:873301/135‑1 (MQ=255)
                                                                                          aGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGc                                         >  2:873301/1‑135 (MQ=255)
                                                                                           gCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAg                                    >  1:415885/1‑139 (MQ=255)
                                                                                            cGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAgg                                   <  1:891105/139‑1 (MQ=255)
                                                                                             gAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGc                                  <  2:487167/139‑1 (MQ=255)
                                                                                                              tAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTaaacaaa                      >  2:384476/1‑134 (MQ=255)
                                                                                                              tAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTaaacaaa                      <  1:384476/134‑1 (MQ=255)
                                                                                                                   ttGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATg                                                                                                       >  1:5262/1‑48 (MQ=255)
                                                                                                                   ttGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATg                                                                                                       <  2:5262/48‑1 (MQ=255)
                                                                                                                     gCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGa          >  1:623021/1‑139 (MQ=255)
                                                                                                                     gCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCACGCCGTCGa          >  2:756755/1‑139 (MQ=255)
                                                                                                                       tGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGc        >  2:31736/1‑139 (MQ=255)
                                                                                                                           ttCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGcgtcg    <  2:213543/139‑1 (MQ=255)
                                                                                                                            tCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAAcc                    >  1:88544/1‑122 (MQ=255)
                                                                                                                            tCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAAcc                    <  2:88544/122‑1 (MQ=255)
                                                                                                                             cTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGCGTCGtt  <  2:491082/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                
GACTGAAAAAGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGCGTCGTT  >  NC_000913/969217‑969480

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