Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 672,110 G→A 8.2% T26I (ACC→ATC)  lptE ← lipopolysaccharide assembly protein LptE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913672,1100GA 8.2% 115.1 / 2.7 49T26I (ACC→ATC) lptElipopolysaccharide assembly protein LptE
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base A (2/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

TGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACT  >  NC_000913/671968‑672246
                                                                                                                                              |                                                                                                                                         
tGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGGTCATGACCTTCATAGTGGAAGGAACgtgggtgg                                                                                                                                           >  1:444116/1‑143 (MQ=255)
 ggcgcGTCCTTTCGCGTGTTTTCTTTCTCAGGCCACGCGCCTCCATTCGGACGTAGCGGGTTACGCACCGCACGGCTTATTGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCTTTGTGGAAGGAACCTGCGTGGt                                                                                                                                          <  2:700858/139‑1 (MQ=255)
 gaaaCCTCCTTCCGCGTCGTTTTTTTATCAACCAACCAAACACCATTCAGTCGTATCTGTTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGt                                                                                                                                          <  1:166922/141‑1 (MQ=255)
         cTTACGCGTGGTGTCTTTATCCAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTGAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCa                                                                                                                                  <  2:624697/143‑1 (MQ=255)
          ttACGCGTGTTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAg                                                                                                                                 <  2:31708/143‑1 (MQ=255)
          ttACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAg                                                                                                                                 <  1:653612/143‑1 (MQ=255)
          ttACGCGTGGTTTCTGTCTCAGGCACGTGCACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAg                                                                                                                                 <  2:655411/143‑1 (MQ=255)
            aCGCGTGGTTTCTTTATCAAGCAACTCGCCACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGAt                                                                                                                               <  1:649184/143‑1 (MQ=255)
                 tGGTTTCTTTATCCAGCGACTCGACACCATTCAGACGTAACTGGAAACCCACCGCACCACTTAATGGCCCGTTCGGATCGCCTGAGTCCAGCATCATGACCTTCATATTGGAAGGAACCTGCGTGGTATCACGCAGATgccag                                                                                                                          <  1:126827/143‑1 (MQ=255)
                 tGGTTTCTTTAGCAAGCACCTCGACACCAGTCAGACGGAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACGTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgccag                                                                                                                          <  1:394966/143‑1 (MQ=255)
                      tCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGGGGTATCACGCAGATgccagccgcg                                                                                                                     >  2:429841/1‑140 (MQ=255)
                        tttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCCGCGACAcc                                                                                                                   >  1:760661/1‑143 (MQ=255)
                              gaGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGTCCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcgg                                                                                                             <  1:471006/142‑1 (MQ=255)
                                 caccTCGACACCATTCCGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAATGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                          <  2:425773/140‑1 (MQ=255)
                                  aaCTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGGTCGGCTGAGTGCCGCGTCAAGACCTTCATAGCGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCGACACAcggcggcacc                                                                                                         >  2:322069/1‑139 (MQ=255)
                                       gACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGCTGCCAGCCACACCCGGCGGTGCTTAa‑cc                                                                                                   >  1:852517/1‑142 (MQ=255)
                                            cATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAGCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGcc                                                                                              >  1:642056/1‑143 (MQ=255)
                                             aTTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGCGTCCAGGATCATGCCCTTCATAGTGGAAGGCACCTGCGTGGTATCACGCCGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCa                                                                                             >  2:429293/1‑143 (MQ=255)
                                               ccAGACGTAACTGGTTGCGCCCGCCACGGCTTGACGGCCCGTGCGGATCGCCTGAGGCCAGGAGCAGGACCTTCAGAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCaga                                                                                           <  2:797182/142‑1 (MQ=255)
                                                cAGCCGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCTGGCGGTG‑TTAATCACCGccgggt                                                                                          <  2:592613/143‑5 (MQ=255)
                                                 agccGTCACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTAGCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAgaga                                                                                         <  2:115598/140‑1 (MQ=255)
                                                 aGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCCGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCAGGCAGATGCCAGCCAGCCCCGGCGGTGGTCAC‑GACCGCCAgaga                                                                                         >  2:348475/1‑143 (MQ=255)
                                                  gACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGAGAGGTATCACGCAGATGCCAGCCACACCCGGCGGCGAAGCA‑CACAGCCAGAGAt                                                                                        >  2:688479/1‑143 (MQ=255)
                                                      tAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGGCCAGGATCATGAGCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGCGATaaca                                                                                    >  2:96083/1‑143 (MQ=255)
                                                              tACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGtt                                                                            <  2:763897/143‑1 (MQ=255)
                                                               aCGCTCCGCAGGGCTGAAGGGCCCGTTGGGGTCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCACCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTg                                                                           <  2:409748/143‑1 (MQ=255)
                                                                                 tGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtc                                                         <  1:726755/143‑1 (MQ=255)
                                                                                 tGGCCCGGTCGGATCGCCTGAGGCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtc                                                         <  1:500966/143‑1 (MQ=255)
                                                                                 cggCCCGTTCGGATCGCCTGAGCCCAGGATCATGACCTTCATAGTGGAGGGAACCTGCGTGGTATCACGCAGATGCCACCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtc                                                         <  1:635688/142‑1 (MQ=255)
                                                                                       gtgctgATCGCCTGAGTCCAGGATCATGACCTTCATCGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCcg                                                   <  2:204070/138‑1 (MQ=255)
                                                                                             aTCGCCTGAGTCCAGGATCATGACCTTCATCGTGGAAGGAACCTGCCTGGTATCACGCAGGTGCCCGCCCCACCCGGCGGCGCTTAA‑CACCGCCAGAGATAACAACAAGGTTGCCAGATATCGCACGCTCCCCCCCGCGCTTa                                             >  1:341191/1‑143 (MQ=255)
                                                                                                  ctggtcccAGGATGAGGACCTTCATCGTGTAATGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATCACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCaa                                        <  2:896288/137‑1 (MQ=255)
                                                                                                  cTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCaa                                        >  1:198148/1‑143 (MQ=255)
                                                                                                    gAGCCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTAGCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACg                                      <  2:581836/143‑1 (MQ=255)
                                                                                                         cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCCGATGCCAGCCACACCCGGCGGTGCTGAA‑CACCGCCAGAGATAACAACAATGTTGCGAGATAGCGCGCGATTCCGCAGGCGCGTAGGCAAGGACCAg                                 >  2:639488/1‑143 (MQ=255)
                                                                                                         cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCAACACCCGCCAGTGGCTTT‑TAGCTTCAGAGATAACAACAATGTTGCCCGACATCGCACGCCTCCTCCCGCGCGTCGCCACCGCCCAg                                 >  1:672608/1‑143 (MQ=255)
                                                                                                           ggATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGAGAGGAA‑CACCGCCAGAGATAACAACAATGTTGCCAGCGACCGCACGCTTCCTCCCGCACCTCACCAACgcctgggg                               >  2:541122/1‑136 (MQ=255)
                                                                                                               cATGACCTTCATAGTGGCAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATGAA‑CACCGCCAGAGATAACAACAATGGCGCCAGATAGCGCACGCCTCCCCCCCCGCATAGCCAGCGACCCGaccgaa                           >  2:383290/1‑138 (MQ=255)
                                                                                                                 tGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACCACCAGACCgagga                         >  2:395680/1‑143 (MQ=255)
                                                                                                                         aTAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCt                 <  2:456390/143‑1 (MQ=255)
                                                                                                                           aGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTgg               >  1:508187/1‑143 (MQ=255)
                                                                                                                            gTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATACCAACAATGTTGGCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGACTGAAGAGGTTACCTGGt              >  2:548597/1‑143 (MQ=255)
                                                                                                                             tGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGCTa             >  1:546406/1‑143 (MQ=255)
                                                                                                                                   gaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCCCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTAACCCCCGGATTCGCCACCGAACAGATTGTGGATGTTACCTGGAACGTcaa       >  1:125358/1‑143 (MQ=255)
                                                                                                                                    aaCCTGCGTGGTATCACGCATATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGGTGCCAGATATCGCACGCTTCCGCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAt      <  1:115321/143‑1 (MQ=255)
                                                                                                                                      tcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCa    >  2:560882/2‑143 (MQ=255)
                                                                                                                                      ccTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCa    <  2:793423/143‑1 (MQ=255)
                                                                                                                                       cTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCGTGGCCGAAGACCAGATTGAGGAGTTTACCTGGTACGTAAATCAc   >  1:126236/1‑143 (MQ=255)
                                                                                                                                       cTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCAAGAAATGACAACACTGTTGCCAGATGTCGGACGCTTTCTCCCGCGAGTAGCCAGCGGCCAGGTTGAGCAGTGCACGTGATACGTAAATcgc   >  2:512997/1‑141 (MQ=255)
                                                                                                                                        tGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCCGGTACGTCGATCACt  >  2:191916/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                         
TGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAA‑CACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACT  >  NC_000913/671968‑672246

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