Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 552,699 T→G 6.4% R155R (CGA→CGC lpxH ← UDP‑2,3‑diacylglucosamine pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb552,6990TG 6.4% 148.7 / 3.2 61R155R (CGA→CGClpxHUDP‑2,3‑diacylglucosamine pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base T (23/34);  new base G (2/2);  total (25/36)
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

CGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGG  >  W3110S.gb/552632‑552768
                                                                   |                                                                     
cgTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTcgc                                                                    <  1:6765880/71‑1 (MQ=255)
cgTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTcgc                                                                    <  1:5652929/71‑1 (MQ=255)
 gTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTcgcg                                                                   <  1:3012005/71‑1 (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:5055276/1‑71 (MQ=255)
        ttAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAAt                                                            >  1:3580176/1‑71 (MQ=255)
        ttAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTCCGCGCGAATGcgcg                                                                   <  1:7907390/64‑1 (MQ=255)
        ttAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTCCGCGCGAATGcgcg                                                                   <  1:7898436/64‑1 (MQ=255)
          aaCGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTCCGCGCGAATGcgcgcg                                                                 >  1:2605646/1‑64 (MQ=255)
          aaCGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTCCGCGCGAATGcgcgcg                                                                 >  1:3645133/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:7254085/71‑1 (MQ=255)
              tCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGttt                                                      <  1:393276/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:2902101/1‑50 (MQ=255)
                    aTCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTcg                                                                     >  1:4126424/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:616993/70‑1 (MQ=255)
                                   ctgctgTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAggg                                  <  1:1247436/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:3718899/69‑1 (MQ=255)
                                                   tGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCt                   <  1:7722934/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)
                                                            cgcgcATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGt        <  1:6341078/71‑1 (MQ=255)
                                                             gcgcATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGtt       >  1:5324710/1‑71 (MQ=255)
                                                             gcgcATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGtt       <  1:1486936/71‑1 (MQ=255)
                                                              cgcATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGttt      <  1:7487762/71‑1 (MQ=255)
                                                                cATTCGCGCGGAAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtgt    <  1:7527626/71‑1 (MQ=255)
                                                                 aTTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtgtg   <  1:4480393/71‑1 (MQ=255)
                                                                 aTTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtgtg   <  1:3993459/71‑1 (MQ=255)
                                                                 aTTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtgtg   <  1:4875675/71‑1 (MQ=255)
                                                                 aTTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtgtg   >  1:5733672/1‑71 (MQ=255)
                                                                 aTTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtgtg   <  1:7784596/71‑1 (MQ=255)
                                                                 aTTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtgtg   <  1:1239118/71‑1 (MQ=255)
                                                                 aTTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTtgtg   <  1:1128668/71‑1 (MQ=255)
                                                                  ttCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTgg  <  1:1826075/71‑1 (MQ=255)
                                                                   |                                                                     
CGTTTTGGTTAACGTCCATGATCGCCAGCGATTTACTGCTGTTGGCTTCTTTGCTGTTCGCGCGCATTCGCGCGGCAATGCGTTTGCGCACAAACAACGGCAGGGCGAGGAATAGCGTCTGCAGCCAGGGTTTGTGG  >  W3110S.gb/552632‑552768

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: