Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,962,311 C→A 14.6% R346S (CGC→AGC)  wecA → UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,962,3110CA14.6% 116.5 / 4.1 41R346S (CGC→AGC) wecAUDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
Reads supporting (aligned to +/- strand):  ref base C (18/17);  new base A (3/3);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

TGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCG  >  CP009273/3962072‑3962522
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                   
tGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTTACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTAGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTNTTTGCTAGCATTCTTCCTCTATGGAGATTGCATTAAGCGTGCCTGGAAAGTTGCTAGCTTTATTAAGCGCTTAAAACGCAGACTGCGTAg                                                                                                                                                                                   >  1:493718/1‑274 (MQ=255)
                                                       aaTGCGTGGGGGTGNTACTTCCCGGAAGGCGTTTGTGCTGATTTCCCTTGCCGCATCACTGCTAGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCTGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTAGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAAcc                                                                                                                           <  2:493710/274‑1 (MQ=255)
                                                                                        tgtgCTGATTACCCTTGCCGCAGCACTGCGCGCTTCCATTGGTGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGNAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGaaaa                                                                                          >  1:473197/1‑275 (MQ=255)
                                                                                                                    cTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGt                                                                                                                                                                                     >  2:547509/1‑156 (MQ=255)
                                                                                                                    cTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGaa                                                                                                                                                                                 >  1:597243/1‑160 (MQ=255)
                                                                                                                        cTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGaa                                                                                                                                                                                 <  1:547520/156‑1 (MQ=255)
                                                                                                                        cTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCg                                                                                                                                                                             <  2:597232/160‑1 (MQ=255)
                                                                                                                                                      aTTTTGGCCCGGTGTGGGTCATGCTGGTGCGCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCGGGAAAGTTGCGCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGAGATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCa                            <  1:93833/275‑1 (MQ=255)
                                                                                                                                                        gtggttacGGGGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTTCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCa                            <  2:473189/266‑1 (MQ=255)
                                                                                                                                                                      ggTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTAGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACAAATGCCTGGGAAACCGGCCGAAGACGCt                                                                                               >  1:640446/1‑192 (MQ=255)
                                                                                                                                                                          aTGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTAGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACAAATGCCTGGGAAACCGGCCGAAGACGCTGaaa                                                                                           <  2:640435/192‑1 (MQ=255)
                                                                                                                                                                                                          ctATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTaa                                                                                                                                                                                                           >  1:163533/1‑48 (MQ=255)
                                                                                                                                                                                                             tGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCATACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACAAAATAAATGAGGATGtgat                                                                                                                                    >  2:165293/1‑116 (MQ=255)
                                                                                                                                                                                                             tGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTggg                                                                                                                 >  1:396514/1‑135 (MQ=255)
                                                                                                                                                                                                              ggATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAgcgc                                                                                                                                                                                                       <  2:163530/48‑1 (MQ=255)
                                                                                                                                                                                                                 tatTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCATACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACAAAATAAATGAGGATGTGATGaca                                                                                                                                <  1:165296/116‑1 (MQ=255)
                                                                                                                                                                                                                 tatTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAAc                                                                                                             <  2:396509/135‑1 (MQ=255)
                                                                                                                                                                                                                  atTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGgttgt                                                                   >  1:103033/1‑176 (MQ=255)
                                                                                                                                                                                                                   ttGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCaaat                                                                                                                                                   >  1:111646/1‑95 (MQ=255)
                                                                                                                                                                                                                      cATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCg                                                               <  2:103032/176‑1 (MQ=255)
                                                                                                                                                                                                                       aTTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCaaataaat                                                                                                                                               <  2:111645/95‑1 (MQ=255)
                                                                                                                                                                                                                       aTTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTAcc                                                           >  2:508418/1‑179 (MQ=255)
                                                                                                                                                                                                                       aTTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAAt      >  2:35373/1‑232 (MQ=255)
                                                                                                                                                                                                                         tAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTaaaa                                                                                                                                                                                                 >  1:134797/1‑43 (MQ=255)
                                                                                                                                                                                                                           aGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCgtgg                                                                                                                                                                          >  2:589262/1‑64 (MQ=255)
                                                                                                                                                                                                                           aGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGAtg                                                                                                                                        >  2:316967/1‑98 (MQ=255)
                                                                                                                                                                                                                           aGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTtgt                                                       <  1:508428/179‑1 (MQ=255)
                                                                                                                                                                                                                           aGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATcgcg  <  1:35373/232‑1 (MQ=255)
                                                                                                                                                                                                                           aGCGTGCCTGGAAAGTTGCTAGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGc                                                                                                                      >  1:127397/1‑116 (MQ=255)
                                                                                                                                                                                                                             cGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTaaaa                                                                                                                                                                                                 >  2:121615/1‑39 (MQ=255)
                                                                                                                                                                                                                             cGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCa                                                                                                                                                                                             <  2:134796/43‑1 (MQ=255)
                                                                                                                                                                                                                               tGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGc                                                                                                                                                                      <  1:589273/64‑1 (MQ=255)
                                                                                                                                                                                                                               tGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGtgat                                                                                                                                    <  1:316972/98‑1 (MQ=255)
                                                                                                                                                                                                                               tGCCTGGAAAGTTGCTAGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTgg                                                                                                                  <  2:127396/116‑1 (MQ=255)
                                                                                                                                                                                                                                gCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACt                                                                                                                                                                                         >  2:323661/1‑44 (MQ=255)
                                                                                                                                                                                                                                 ccTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCa                                                                                                                                                                                             <  1:121616/39‑1 (MQ=255)
                                                                                                                                                                                                                                   tGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACt                                                                                                                                                                                         >  1:67468/1‑41 (MQ=255)
                                                                                                                                                                                                                                    ggAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGt                                                                                                                                                                                     <  1:323666/44‑1 (MQ=255)
                                                                                                                                                                                                                                       aaGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGt                                                                                                                                                                                     <  2:67468/41‑1 (MQ=255)
                                                                                                                                                                                                                                           tGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGtggtgg                                                                                                                                                                       >  2:117994/1‑51 (MQ=255)
                                                                                                                                                                                                                                            gCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATg                                                                                        >  2:177687/1‑129 (MQ=255)
                                                                                                                                                                                                                                               cGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGc                                                                                                                                                                   <  1:117995/51‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                   
TGCGTAAAGGCATGAGCCCATTCTCTCCTGACCGTCAGCATATTCACCATTTGATCATGCGTGCCGGGTTTACTTCCCGTCAGGCGTTTGTGCTGATTACCCTTGCCGCAGCACTGCTCGCTTCCATTGGCGTGCTGGCAGAATATTCTCATTTTGTCCCGGAGTGGGTCATGCTGGTGCTCTTTTTGCTAGCATTCTTCCTCTATGGATATTGCATTAAGCGTGCCTGGAAAGTTGCTCGCTTTATTAAGCGCGTAAAACGCAGACTGCGTAGAAATCGTGGTGGCAGCCCCAATTTAACCAAATAAATGAGGATGTGATGACACAACCAATGCCTGGGAAACCGGCCGAAGACGCTGAAAATGAACTGGATATTCGTGGGTTGTTTCGTACCTTGTGGGCTGGGAAGCTATGGATTATTGGCATGGGGCTGGCGTTTGCGTTAATCGCG  >  CP009273/3962072‑3962522

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑