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,3920AC77.1% 72.9 / 13.0 35Y240D (TAC→GAC) lptDLPS assembly protein LptD
Reads supporting (aligned to +/- strand):  ref base A (5/3);  new base C (12/15);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

GCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGCTAT  >  NZ_CP009273/56306‑56475
                                                                                      |                                                                                   
gCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGtt                                                                                  <  2:724403/90‑1 (MQ=255)
   acgacgATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTtggt                                                                               <  2:87238/90‑1 (MQ=255)
        gaTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTtggtggtgg                                                                          <  1:238655/90‑1 (MQ=255)
        gaTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTtggtggtgg                                                                          <  2:351286/90‑1 (MQ=255)
        gaTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTtggtggtgg                                                                          <  1:29288/90‑1 (MQ=255)
               ataATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACtt                                                                   <  2:144373/90‑1 (MQ=255)
                  aTGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTCGTTGGTGGGGGGGTACTTGGc                                                                >  1:282368/1‑90 (MQ=255)
                       gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCg                                                           >  1:134341/1‑90 (MQ=255)
                       gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAAACGTTGGTGGGGGTGGACTTGGGGGTCg                                                           >  2:478856/1‑90 (MQ=255)
                       gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAATCGTTGGGGGGGGGGGACTTTGCGGTCg                                                           >  2:552646/1‑90 (MQ=255)
                        cgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACCCAAAGTCGTTGGTGGTGGGGTACTTGGCGTTCgg                                                          >  2:339013/1‑90 (MQ=255)
                          tGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCgaga                                                        <  2:653812/90‑1 (MQ=255)
                          tGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGa                                                        <  2:806669/90‑1 (MQ=255)
                            aTGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAGTCGTTGGGGGGGGGGGACTTTggggtcggggtc                                                      >  2:325192/1‑90 (MQ=255)
                                      ccATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGGGGGGGGGGACTTGGGGTTCGGGGTCAAGAAACCag                                            >  1:370318/1‑90 (MQ=255)
                                         tatTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGc                                         <  2:194568/90‑1 (MQ=255)
                                             tGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGGGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACg                                     >  1:762565/1‑90 (MQ=255)
                                               gcgcGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGc                                         >  1:805506/1‑84 (MQ=255)
                                               gcgcGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGc                                         <  2:805506/84‑1 (MQ=255)
                                               gcgcGATGTTCCAGTAATATGGCAGGTAGAACTCAAAAGAGTTGGGGGGGGGGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                                   >  2:152100/1‑90 (MQ=255)
                                                              aaTATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCa                    <  1:325192/90‑1 (MQ=255)
                                                              aaTATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCa                    <  1:478856/90‑1 (MQ=255)
                                                               atatGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCaa                   <  1:794569/90‑1 (MQ=255)
                                                                tatGGCAGGTAGAACTCAAAGTAGTTGGGGGGGGGGTACTTGGCGTTCGGGGGCAAGAAAACAGAGCGACGTTTGTCACCCACCGGCAAc                  >  2:532489/1‑90 (MQ=255)
                                                                  tGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTg                >  1:349407/1‑90 (MQ=255)
                                                                    gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCa              >  1:249655/1‑90 (MQ=255)
                                                                      aGGGAGAACTCAAAGTCGTTGGGGGTGGGGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCaaa            >  2:552238/1‑90 (MQ=255)
                                                                         tAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAg         <  2:385463/90‑1 (MQ=255)
                                                                         tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAg         <  1:532489/90‑1 (MQ=255)
                                                                         tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGACGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAg         <  2:319855/90‑1 (MQ=255)
                                                                         tAGAACTCAAAATCGTTGGTGGGGGGGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAg         >  1:269509/1‑90 (MQ=255)
                                                                             aCTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGc     <  1:249117/90‑1 (MQ=255)
                                                                                cAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGCtat  >  1:546300/1‑90 (MQ=255)
                                                                                  aaGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAg         >  2:754025/1‑81 (MQ=255)
                                                                                  aaGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAg         <  1:754025/81‑1 (MQ=255)
                                                                                      |                                                                                   
GCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAGGGGCTAT  >  NZ_CP009273/56306‑56475

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