Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,144,392 A→C 22.3% D233A (GAC→GCC)  lpxK → tetraacyldisaccharide 4'‑kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,144,3920AC22.3% 34.9 / 12.1 27D233A (GAC→GCC) lpxKtetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base A (8/13);  new base C (3/3);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

GGGCCGTTGCGTGAGCCGGCCGACCGCTTGCGTGAGGCTGACGCGGTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCC  >  NC_002947/2144244‑2144525
                                                                                                                                                     |                                                                                                                                     
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                                                                                                                             <  2:141400/109‑1 (MQ=255)
                                                  gTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTcccc                                                                                                                             >  1:141400/1‑109 (MQ=255)
                                                  gTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAAccc                                                                                    >  2:398033/1‑150 (MQ=255)
                                                   ttCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCACCGGCCACCGCATGCAGGCGCTGGCCggcgg                                                                                             <  2:239310/140‑4 (MQ=255)
                                                   ttCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCACCGGCCACCGCATGCAGGCGCTGGCCggcgg                                                                                             >  1:239310/1‑137 (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)
                                                                                                                                    gAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGcc  <  1:108851/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                     
GGGCCGTTGCGTGAGCCGGCCGACCGCTTGCGTGAGGCTGACGCGGTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCC  >  NC_002947/2144244‑2144525

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: