Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 552,689 T→C 5.8% N159D (AAC→GAC)  lpxH ← UDP‑2,3‑diacylglucosamine pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048552,6890TC 5.8% 191.2 / 2.4 69N159D (AAC→GAC) lpxHUDP‑2,3‑diacylglucosamine pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base C (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAG  >  AP009048/552620‑552755
                                                                     |                                                                  
cACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTc                                                                   >  1:7459948/1‑71 (MQ=255)
 aCTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTcg                                                                  <  1:7935941/71‑1 (MQ=255)
  cTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTcgc                                                                 <  1:1227738/71‑1 (MQ=255)
  cTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTcgc                                                                 <  1:1784499/71‑1 (MQ=255)
    gACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTcgcgc                                                               >  1:6625022/1‑71 (MQ=255)
    gACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTcgcgc                                                               >  1:1940296/1‑71 (MQ=255)
    gACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTcgcg                                                                >  1:7135079/1‑70 (MQ=255)
    gACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTc                                                                   >  1:7882603/1‑67 (MQ=255)
    gACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTc                                                                   >  1:5087390/1‑67 (MQ=255)
     accaccGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGtt                                                                    >  1:6805748/1‑65 (MQ=255)
     accaccGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTcgcgcg                                                              >  1:162621/1‑71 (MQ=255)
     accaccGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTcgcgcg                                                              >  1:7081419/1‑71 (MQ=255)
     accaccGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTcgcgcg                                                              >  1:2784148/1‑71 (MQ=255)
     accaccGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTcgcgcg                                                              >  1:4318660/1‑71 (MQ=255)
            cgTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTcgc                                                       <  1:6765880/71‑1 (MQ=255)
            cgTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTcgc                                                       <  1:5652929/71‑1 (MQ=255)
             gTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTcgcg                                                      <  1:3012005/71‑1 (MQ=255)
              ttttGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTcgc                                                                 <  1:851484/59‑1 (MQ=255)
              ttttGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTcgc                                                                 <  1:1137482/59‑1 (MQ=255)
               tttGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTcgc                                                                 <  1:2433178/58‑1 (MQ=255)
                ttGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTcgcgc                                                               >  1:3931654/1‑59 (MQ=255)
                ttGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCgg                                                   >  1:7458232/1‑71 (MQ=255)
                ttGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGAGCGCATTcgcgcg                                                    <  1:6592843/70‑1 (MQ=255)
                   gTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGAGGCa                                                 <  1:471126/70‑1 (MQ=255)
                    ttAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAAt                                               >  1:3580176/1‑71 (MQ=255)
                    ttAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAAt                                               >  1:5055276/1‑71 (MQ=255)
                    ttAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTCCGCGCGAATGcgcg                                                      <  1:7898436/64‑1 (MQ=255)
                    ttAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTCCGCGCGAATGcgcg                                                      <  1:7907390/64‑1 (MQ=255)
                     tAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTcgcgc                                                               <  1:3924219/54‑1 (MQ=255)
                      aaCGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTCCGCGCGAATGcgcgcg                                                    >  1:3645133/1‑64 (MQ=255)
                      aaCGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTCCGCGCGAATGcgcgcg                                                    >  1:2605646/1‑64 (MQ=255)
                        cGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTcg                                                        <  1:2866831/58‑1 (MQ=255)
                         gTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCaa                                                >  1:4252341/1‑65 (MQ=255)
                         gTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGc                                             <  1:1323444/68‑1 (MQ=255)
                          tCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGc                                             <  1:3062142/67‑1 (MQ=255)
                          tCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGttt                                         <  1:393276/71‑1 (MQ=255)
                          tCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGttt                                         <  1:7254085/71‑1 (MQ=255)
                            cATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTgc                                       <  1:7340435/71‑1 (MQ=255)
                             aTGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTgcg                                      <  1:3994471/71‑1 (MQ=255)
                              tGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGc                                                  <  1:2033591/58‑1 (MQ=255)
                              tGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTgcg                                      >  1:3908374/1‑70 (MQ=255)
                               gATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCTCGCGCATTCGCGCGGCAATGCGTTTgcg                                      <  1:4608735/69‑1 (MQ=255)
                                aTCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCAtt                                                          >  1:3242682/1‑48 (MQ=255)
                                aTCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCAtt                                                          >  1:6474855/1‑48 (MQ=255)
                                aTCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTcg                                                        >  1:4126424/1‑50 (MQ=255)
                                aTCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTcg                                                        >  1:2902101/1‑50 (MQ=255)
                                aTCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGca                                    >  1:3566630/1‑70 (MQ=255)
                                     cAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCAtt                                                          <  1:4520274/43‑1 (MQ=255)
                                      aGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCacaaacaa                             >  1:7600056/1‑71 (MQ=255)
                                      aGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCacaaacaa                             >  1:7747204/1‑71 (MQ=255)
                                      aGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCacaaacaa                             >  1:7848802/1‑71 (MQ=255)
                                       gCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTcgc                                                       >  1:7721051/1‑44 (MQ=255)
                                        cGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCaa                                                >  1:6284690/1‑50 (MQ=255)
                                        cGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCaa                                                >  1:5678433/1‑50 (MQ=255)
                                        cGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACg                           >  1:7587482/1‑71 (MQ=255)
                                           tttACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACgg                          >  1:5596694/1‑69 (MQ=255)
                                             tACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTgcg                                      >  1:3708174/1‑55 (MQ=255)
                                              aCTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAgg                      <  1:5505991/70‑1 (MQ=255)
                                               ctgctgTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAggg                     <  1:1247436/70‑1 (MQ=255)
                                               ctgctgTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAggg                     <  1:616993/70‑1 (MQ=255)
                                                tgctgTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCacaaacaa                             >  1:4990967/1‑61 (MQ=255)
                                                    gTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAgg                >  1:7221497/1‑70 (MQ=255)
                                                     ttGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAgg                <  1:7378882/69‑1 (MQ=255)
                                                     ttGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGaa              <  1:3831776/71‑1 (MQ=255)
                                                       ggCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAAt             <  1:4502916/70‑1 (MQ=255)
                                                             tttGCTGTTCGCGCGCATTCGCGCGGCAATGCGttt                                         <  1:4162465/36‑1 (MQ=255)
                                                               tGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCt      <  1:7722934/69‑1 (MQ=255)
                                                               tGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCt      <  1:3718899/69‑1 (MQ=255)
                                                                gCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGc    <  1:1824578/70‑1 (MQ=255)
                                                                 cTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAg  >  1:6805371/1‑71 (MQ=255)
                                                                     tCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAg  >  1:1450841/1‑67 (MQ=255)
                                                                     |                                                                  
CACTGACCACCGCGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAG  >  AP009048/552620‑552755

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