Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,144,372	C→A	11.3%	L241M (CTG→ATG)	lpxK →	tetraacyldisaccharide 4'‑kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,144,372	0	C	A	11.3%	132.3 / 5.1	44	L241M (CTG→ATG)	lpxK	tetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base C (21/18);  new base A (2/3);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CTGACGCGGTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGCCACTGCCGCTGGTC  >  NC_002947/2144238‑2144518                                                                                                                                       |                                                                                                                                                   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)                             cGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCACCGCATGCAGGCGCTGGCCggcgggg                                                                                                                               >  4:224212/1‑123 (MQ=39)                             cGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCACCGCATGCAGGCGCTGGCCggcgggg                                                                                                                               <  3:224212/128‑6 (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                                                                                             >  7:79285/1‑112 (MQ=255)                                                                               tCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTggg                                                                                             <  8:79285/112‑1 (MQ=255)                                                                                          aCTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGtt                                             >  4:15628/1‑149 (MQ=255)                                                                                            tGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGtt                                             >  2:157313/1‑147 (MQ=255)                                                                                            tGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGtt                                             >  5:370501/1‑147 (MQ=255)                                                                                                          cgcTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTg                             <  7:297929/149‑1 (MQ=255)                                                                                                          cgcTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTg                             <  4:78859/149‑1 (MQ=255)                                                                                                          ccatctaCCCTTTCCCCGCCGGCCAGCGCCTGCAGGCGGGGGCCGGTATCCGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGCCCATGCGCAGTTCAGTGCCCGGAGCCTg                             <  4:230004/143‑1 (MQ=255)                                                                                                           gcTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACg                                                                                                                  <  2:301466/63‑1 (MQ=255)                                                                                                           gcTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACg                                                                                                                  >  1:301466/1‑63 (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)                                                                                                                                  accgcATGCAGGCGCTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGCCACTGCCGCTg     <  1:295767/147‑1 (MQ=255)                                                                                                                                  aGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGCCACTGCCGCTg     >  1:312988/1‑149 (MQ=255)                                                                                                                                   gcgcCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGCCACTGCCGCTgg    <  7:132712/149‑1 (MQ=255)                                                                                                                                     gcCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGCCACTGCCGCTGGTc  <  2:295503/149‑1 (MQ=255)                                                                                                                                      ccTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTCCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGCCACTGCCg        <  7:240592/142‑1 (MQ=255)                                                                                                                                      ccTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTCCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGCCACTGCCg        >  8:240592/1‑142 (MQ=255)                                                                                                                                       |                                                                                                                                                   CTGACGCGGTGCTGTTCAACGGCGCCAGCGAAGACCGCACCGATGGTTTCGGTTTCTGCCTGCAGCCGTCGGCCCTGGTCAACGTGCGCACTGGCGAGCGCCGTGCGCTCGACCATTTCCCCGCCGGCCAGCGCCTGCATGCGGTGGCCGGTATCGGCAACCCGCAACGTTTCTTCAACACCCTGCTGGGGCTAAACTGGCAGCCGGTGCCGCATCCCTTCGCCGACCATGCGCAGTTCAGTGCCCGGAGCCTGGCCTTCAGCCCGCCACTGCCGCTGGTC  >  NC_002947/2144238‑2144518

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: ‑