Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,144,368 G→C 9.8% Q239H (CAG→CAC)  lpxK → tetraacyldisaccharide 4'‑kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,144,3680GC 9.8% 115.0 / 4.1 41Q239H (CAG→CAC) lpxKtetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base G (21/16);  new base C (2/2);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

CTGACGCGGTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGCCACTGCCGCTGG  >  NC_002947/2144238‑2144516
                                                                                                                                  |                                                                                                                                                    
cTGACGCGGTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGcc                                                                                                                                    <  5:325626/149‑1 (MQ=255)
 tGACGCGGTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCg                                                                                                                                   >  1:266433/1‑149 (MQ=255)
     gcgGTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCACCGCATGCAGGCGCTGGCCggcgg                                                                                                                               <  1:135411/149‑4 (MQ=255)
       ggTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCg                                                                                                                             <  5:75846/149‑1 (MQ=255)
           cTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCaa                                                                                                                         >  7:336792/1‑149 (MQ=255)
                 aaCGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGc                                                                                                                    >  5:152360/1‑148 (MQ=255)
                 aaCGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGc                                                                                                                    >  3:304528/1‑148 (MQ=255)
                      cgcCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAgcgc                                                                                                                                                   >  7:187151/1‑112 (MQ=255)
                      cgcCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAgcgc                                                                                                                                                   <  8:187151/112‑1 (MQ=255)
                            cGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCACCGCATGCAGGCGCTGGCCggcgggg                                                                                                                             <  3:224212/128‑6 (MQ=39)
                            cGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCACCGCATGCAGGCGCTGGCCggcgggg                                                                                                                             >  4:224212/1‑123 (MQ=39)
                             gAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCaa                                                                                                       >  4:82847/1‑149 (MQ=255)
                                        cGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTgg                                                                                            <  8:336792/149‑1 (MQ=255)
                                         gATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTggg                                                                                           <  3:243744/149‑1 (MQ=255)
                                                   gtttcTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTgg                                                                                 <  3:82847/149‑1 (MQ=255)
                                                              cAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGcc                                                                      >  5:109468/1‑149 (MQ=255)
                                                              cAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGc                                                                       <  4:304528/148‑1 (MQ=255)
                                                                     cGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATccc                                                               <  6:109468/149‑1 (MQ=255)
                                                                        cccTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACAACCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCATc                                                            >  7:86675/1‑149 (MQ=255)
                                                                          cTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGc                                                          >  6:261544/1‑149 (MQ=255)
                                                                          cTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGc                                                          >  8:132712/1‑149 (MQ=255)
                                                                              tCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTggg                                                                                           <  8:79285/112‑1 (MQ=255)
                                                                              tCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTggg                                                                                           >  7:79285/1‑112 (MQ=255)
                                                                                         aCTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGtt                                           >  4:15628/1‑149 (MQ=255)
                                                                                           tGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGtt                                           >  2:157313/1‑147 (MQ=255)
                                                                                           tGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGtt                                           >  5:370501/1‑147 (MQ=255)
                                                                                                         cgcTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTg                           <  4:78859/149‑1 (MQ=255)
                                                                                                         cgcTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTg                           <  7:297929/149‑1 (MQ=255)
                                                                                                         ccatctaCCCTTTCCCCGCCGGCCAGCGCCTGCAGGCGGGGGCCGGTATCCGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGCCCATGCGCAGTTCAGTGCCCGGAGCCTg                           <  4:230004/143‑1 (MQ=255)
                                                                                                          gcTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACg                                                                                                                >  1:301466/1‑63 (MQ=255)
                                                                                                          gcTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACg                                                                                                                <  2:301466/63‑1 (MQ=255)
                                                                                                            ttgACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGcc                        <  2:159427/147‑1 (MQ=255)
                                                                                                            tCGACCATTTCC‑‑‑CCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCtt                      >  7:167088/1‑148 (MQ=255)
                                                                                                              gACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCt                       >  4:108517/1‑148 (MQ=255)
                                                                                                                ccATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTg                           >  3:230004/1‑142 (MQ=255)
                                                                                                                 cATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAg                   <  1:142755/149‑1 (MQ=255)
                                                                                                                      ccccGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGc              >  8:181084/1‑149 (MQ=255)
                                                                                                                       gataccgGCCACCGCATGCAGGCGCTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGcc             >  2:295767/5‑149 (MQ=255)
                                                                                                                       cccGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGc              <  8:86675/148‑1 (MQ=255)
                                                                                                                           ccgGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGCCACTg         <  6:370501/149‑1 (MQ=255)
                                                                                                                               ccAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGCCACTGCCGc     >  1:295503/1‑149 (MQ=255)
                                                                                                                                cAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGCCACTGCCGCt    >  3:147070/1‑149 (MQ=255)
                                                                                                                                 aGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGCCACTGCCGCTg   >  1:312988/1‑149 (MQ=255)
                                                                                                                                  gcgcCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGCCACTGCCGCTgg  <  7:132712/149‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                                    
CTGACGCGGTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGCCACTGCCGCTGG  >  NC_002947/2144238‑2144516

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑