Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 56,392 A→C Y240D (TAC→GAC)  lptD ← LPS assembly protein LptD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927356,3920AC93.1% 75.5 / ‑0.4 29Y240D (TAC→GAC) lptDLPS assembly protein LptD
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (14/13);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-02

GCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGCTATAAA  >  NZ_CP009273/56306‑56478
                                                                                      |                                                                                      
gCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGtt                                                                                     <  1:722117/90‑1 (MQ=255)
   acgacgATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTtggt                                                                                  <  2:708650/90‑1 (MQ=255)
   acgacgATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTtggt                                                                                  <  2:284888/90‑1 (MQ=255)
      acgaTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGCtggtggt                                                                               <  1:333696/90‑1 (MQ=255)
         aTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTCGTtggtggtggg                                                                            >  2:206224/1‑89 (MQ=255)
         aTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAGTCGTtggggggggg                                                                            >  1:72925/1‑89 (MQ=255)
           gCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGgtgt                                                                          <  1:382018/90‑1 (MQ=255)
              tataATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAATCGTTTGGGGTGGGGTGCt                                                                       >  1:472772/1‑90 (MQ=255)
                 aaTGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTgg                                                                    <  1:462592/90‑1 (MQ=255)
                      ggcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAATCGTTGGTGGGGGTGTACTTGGCGTTc                                                               >  1:261567/1‑90 (MQ=255)
                       gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGGAATCAAAATCGGTGGGGGGGGGGTACTTGGGGGTCg                                                              >  1:10273/1‑90 (MQ=255)
                       gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTCGTTGGTGGGGGTGTACTTGGCGTTCg                                                              >  2:215915/1‑90 (MQ=255)
                       gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAATCGTTGGGGGGGGTGTACTTGGCGTTCg                                                              >  1:198442/1‑90 (MQ=255)
                                    aTCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGGGGGGGGGTACTTTGCGTTCGGGGTCAAGAAAcc                                                 >  2:449267/1‑90 (MQ=255)
                                            ttGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGAc                                         <  2:72925/90‑1 (MQ=255)
                                               gcgcGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGGGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                                      >  2:109764/1‑90 (MQ=255)
                                               gcgcGATGTTCCAGTAATATGGCAGGTAGAACACAAAATCGTTGGTGGGGGTGTACTTTGGGTTTGGGGTCAAGAAAACAGAGAGAAGtt                                      >  1:593533/1‑90 (MQ=255)
                                                              aaTATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCa                       <  1:667761/90‑1 (MQ=255)
                                                              aaTATGGCAGGTAGAACTCAAAATCGGTGGTGGTGGGGTACTTGGGGTTTGGGGTCAAGAAAACAGGGCGACGTTTTTCACCCCCCCGcg                       >  1:207208/1‑89 (MQ=255)
                                                               atatGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCaa                      <  1:984/90‑1 (MQ=255)
                                                                tatGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAAc                     <  2:697205/90‑1 (MQ=255)
                                                                  tGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTg                   <  1:191499/90‑1 (MQ=255)
                                                                    gCAGGTAGAACTCAAAGTCGTTGGTGGTGGGGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCa                 >  2:304102/1‑90 (MQ=255)
                                                                      aGGTAGAACTCAAAGGCGTTGGGGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCaaa               >  1:105637/1‑90 (MQ=255)
                                                                           gAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAggg          <  2:636197/90‑1 (MQ=255)
                                                                             aCTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGc        <  1:318089/90‑1 (MQ=255)
                                                                             aCTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGc        <  1:304102/90‑1 (MQ=255)
                                                                                cAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGAc                                         >  1:330763/1‑54 (MQ=255)
                                                                                cAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGAc                                         <  2:330763/54‑1 (MQ=255)
                                                                                   agtagtTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGCTATaaa  <  1:611495/90‑1 (MQ=255)
                                                                                      |                                                                                      
GCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGCTATAAA  >  NZ_CP009273/56306‑56478

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: