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,3920AC94.7% 45.8 / NA 19Y240D (TAC→GAC) lptDLPS assembly protein LptD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (9/9);  minor base G (1/0);  total (10/9)
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

GACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGCTATAA  >  NZ_CP009273/56311‑56477
                                                                                 |                                                                                     
gacgaTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTtggt                                                                                 <  2:395909/88‑1 (MQ=255)
gacgaTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAAATCAAAGTCGTtggt                                                                                 >  1:395909/1‑88 (MQ=255)
  cgaTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTtggtggtg                                                                             <  1:412757/90‑1 (MQ=255)
    aTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGGGAAATCAAAGGGGTtggtgggggg                                                                           >  2:7843/1‑89 (MQ=255)
             aTGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAAGCGTTGGTGGGGGGGGACTTGGc                                                                  >  2:527787/1‑90 (MQ=255)
              tGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGGGAACTCAAAATCGTTGGGGGGGGTGTAGTTGGCg                                                                 >  2:436488/1‑90 (MQ=255)
                  gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCg                                                             >  1:259354/1‑90 (MQ=255)
                     tGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGa                                                          <  1:578394/90‑1 (MQ=255)
                              cATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAAATCAAAATCGTTGGGGGGGGGGGACTTGGGGTTCGGGGTCAAGAAAc                                                 >  2:578394/1‑90 (MQ=255)
                                    tatTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGc                                           <  1:336785/90‑1 (MQ=255)
                                        tGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAGTCGTTGGGGGGGGGGTACTTGGCGTTCGGGATCAAAAAAACAGAGCGACg                                       >  2:180130/1‑90 (MQ=255)
                                          gcgcGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGGGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                                     >  1:29860/1‑90 (MQ=255)
                                          gcgcGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGGGGGGTACTTGGCGTTCGGGAGCAAAAAACCAGAGCGACggt                                     >  1:590907/1‑88 (MQ=255)
                                              gATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTc                                 <  1:527787/90‑1 (MQ=255)
                                                          atatGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCaa                     <  1:102313/90‑1 (MQ=255)
                                                           tatGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAAc                    <  1:403631/90‑1 (MQ=255)
                                                                    tAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAg           <  2:259354/90‑1 (MQ=255)
                                                                    tAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAg           <  2:279701/90‑1 (MQ=255)
                                                                             aaGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGCTATaa  >  2:630361/1‑90 (MQ=255)
                                                                                 |                                                                                     
GACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGCTATAA  >  NZ_CP009273/56311‑56477

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