Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,144,374 G→T 21.4% G227V (GGC→GTC)  lpxK → tetraacyldisaccharide 4'‑kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,144,3740GT21.4% 52.8 / 7.8 28G227V (GGC→GTC) lpxKtetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base G (9/13);  new base T (3/3);  total (12/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 = 8.39e-01

GGGCCGTTGCGTGAGCCGGCCGACCGCTTGCGTGAGGCTGACGCGGTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCA  >  NC_002947/2144244‑2144520
                                                                                                                                  |                                                                                                                                                   
gggCCGTTGCGTGAGCCGGCCGACCGCTTGCGTGAGGCTGACGCGGTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGAc                                                                                                                                 <  2:16659/150‑1 (MQ=255)
                             gCGTGAGGCTGACGCGGTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGc                                                                                                    <  2:150943/150‑1 (MQ=255)
                               gTGAGGCTGACGCGGTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCgg                                                                                                  <  1:130553/150‑1 (MQ=255)
                                tGAGGCTGACGCGGTGCCGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTCCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGATCGCCGTACGCTCGACCATTTCCCCGCCGGCCAGTCCCTGCATGCgtg                                                                                                 >  2:130553/1‑148 (MQ=255)
                                          gcgGTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGTGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATc                                                                                       >  2:41539/1‑150 (MQ=255)
                                           gtgTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCg                                                                                      <  1:381893/148‑1 (MQ=255)
                                                  gTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTcccc                                                                                                                        >  1:141400/1‑109 (MQ=255)
                                                  gTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTcccc                                                                                                                        <  2:141400/109‑1 (MQ=255)
                                                  gTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAAccc                                                                               >  2:398033/1‑150 (MQ=255)
                                                   ttCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCACCGGCCACCGCATGCAGGCGCTGGCCggcgg                                                                                        >  1:239310/1‑137 (MQ=255)
                                                   ttCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCACCGGCCACCGCATGCAGGCGCTGGCCggcgg                                                                                        <  2:239310/140‑4 (MQ=255)
                                                       aCGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAgcgcacg                                                                                                                                           <  2:89279/86‑4 (MQ=255)
                                                       aCGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAgcgcacg                                                                                                                                           >  1:89279/1‑83 (MQ=255)
                                                                                        ttttcTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGc                                         <  2:334098/149‑1 (MQ=255)
                                                                                               cctgcAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGt                                  >  2:108851/1‑150 (MQ=255)
                                                                                                        gTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATcc                         <  1:397789/150‑1 (MQ=255)
                                                                                                             cccTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCg                    >  1:94472/1‑150 (MQ=255)
                                                                                                              ccTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGAGCCGCATCCCTTCGc                   >  2:117512/1‑150 (MQ=255)
                                                                                                                tggccggcggggaaatgGTCGAGCGCACGGCGCTCGCCCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTgggg                                                   >  2:265089‑M1/18‑116 (MQ=255)
                                                                                                                tggccggcggggaaatgGTCGAGCGCACGGCGCTCGCCCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTgggg                                                   <  1:265089‑M1/99‑1 (MQ=255)
                                                                                                                     ggcggggaaatgGTCGAGCGCACGGCGCTCGCCCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTgggg                                                   >  2:360595‑M1/13‑111 (MQ=255)
                                                                                                                     ggcggggaaatgGTCGAGCGCACGGCGCTCGCCCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTgggg                                                   <  1:360595‑M1/99‑1 (MQ=255)
                                                                                                                       cGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATgc          <  1:117512/150‑1 (MQ=255)
                                                                                                                       cGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCAGGTGCCGCATCCCTTCGCCGACCATgc          <  1:68577/150‑1 (MQ=255)
                                                                                                                          gcgcACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCa       <  1:398033/150‑1 (MQ=255)
                                                                                                                              atgGTCGAGCGCACGGCGCTCGCCCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTgg                                                     >  1:479897‑M1/4‑100 (MQ=255)
                                                                                                                              atgGTCGAGCGCACGGCGCTCGCCCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTgg                                                     <  2:479897‑M1/97‑1 (MQ=255)
                                                                                                                               cTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCa  <  2:73725/150‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                                   
GGGCCGTTGCGTGAGCCGGCCGACCGCTTGCGTGAGGCTGACGCGGTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCA  >  NC_002947/2144244‑2144520

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