Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927356,3790TG59.5% 8.2 / 35.1 37Y244S (TAC→TCC) lptDLPS assembly protein LptD
Reads supporting (aligned to +/- strand):  ref base T (4/11);  new base G (22/0);  total (26/11)
Fisher's exact test for biased strand distribution p-value = 1.60e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACATGATGTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAG  >  NZ_CP009273/56296‑56468
                                                                                   |                                                                                          
acaTGATGTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACt                                                                                      >  2:414418/1‑90 (MQ=255)
  atgatgTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGGACTCa                                                                                    >  2:46630/1‑90 (MQ=255)
   tgatgTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCaa                                                                                   <  2:723647/90‑1 (MQ=255)
     atgTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAag                                                                                 <  2:588664/90‑1 (MQ=255)
     atgTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAag                                                                                 >  1:689483/1‑90 (MQ=255)
       gTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAAATCAAagga                                                                               >  1:324330/1‑88 (MQ=255)
       gTTGCCACGACGATGCACATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCCAagta                                                                               >  1:78992/1‑90 (MQ=255)
          gCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAATAGtt                                                                           >  2:618153/1‑90 (MQ=255)
               gacgaTGCATATAATGCGGCGTGATGGGGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGTGGAATCAAAATATGTtggtg                                                                       >  2:72946/1‑90 (MQ=255)
                   aTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTggtggtggt                                                                  <  2:114609/90‑1 (MQ=255)
                             tGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAATAGGTGGGGGTGGGGGACtttggg                                                        >  1:337602/1‑86 (MQ=255)
                             tGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAGTAGTTGGTGGGGGTGGACTTgggg                                                        >  2:310609/1‑88 (MQ=255)
                             tGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAATAAGTGGTGGTGGGGGACTTgggg                                                        >  2:178387/1‑88 (MQ=255)
                             tGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAAAAGTTGGTGGGGGGGTACtttggg                                                        >  2:365137/1‑86 (MQ=255)
                               cggcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGtt                                                      <  1:282453/90‑1 (MQ=255)
                               cggcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGGAATCAAAATAGTTGGGGGGGGTGGACTTTGCggt                                                      >  2:648985/1‑88 (MQ=255)
                               cggcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAGGAGTTGGGGGGGGGGTACTTgggggt                                                      >  2:438865/1‑86 (MQ=255)
                                ggcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTc                                                     <  1:438865/90‑1 (MQ=255)
                                ggcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAAAAGGTGGTGGGGGGGTACTTTGCGTTc                                                     >  2:379874/1‑90 (MQ=255)
                                 gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGGAATCAAAATAGTTGGGGGTGGTGTTCTTGGGGTTCg                                                    >  2:262529/1‑90 (MQ=255)
                                  cgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAAATCAAAAAAGTTGGTGGTGGTGTACTTGGGGGTCgg                                                   >  1:341209/1‑90 (MQ=255)
                                    tGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGa                                                 <  2:501254/90‑1 (MQ=255)
                                           ggCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGGGAACTCAAAAAAGTTGGGGGGGGTGTACTTGGGGCTCGGGGTCAaaaa                                          >  2:488216/1‑90 (MQ=255)
                                                   tatTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGc                                  <  2:41663/90‑1 (MQ=255)
                                                      ttGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGAc                               <  1:414418/90‑1 (MQ=255)
                                                         gcgcGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGGGTACTTGGCGGTTGGGAGCAAGAAAACAAAGCGACggt                            >  2:619587/1‑88 (MQ=255)
                                                         gcgcGATGTTACAGTAATATGGCAGGTGGAAATCAAAAAAGTTGGTGGTGGGGGACTTGGGGGTCGGGGGCAAGAAACCAGAGCGACGtg                            >  2:640175/1‑89 (MQ=255)
                                                             gATGTTCCAGTAATATGGCAGGGAGAACTCAAAATAGTTGGGGGTGGGGTACTTGGCGTTCGGGGTCAAGAAACCAGAGCGACGTTTGTc                        >  2:295383/1‑90 (MQ=255)
                                                             gATGTTCCAGTAATATGGCAGGGAGAACTCAAAAGAGTTGGGGGGGGGGGACTTGGGGGTCGGGGTTAAGAAACCAGAGCGACGGTTTTc                        >  2:101372/1‑90 (MQ=255)
                                                                        aaTATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCa             <  2:768699/90‑1 (MQ=255)
                                                                          tatGGCAGGGAGAAATCAAAGTAGTTGGGGGGGGGGTACTTGGGGTTCGGGATCAAAAAACCAAAGCGACGTTTGTCACCCACCGGCAAc           >  2:135818/1‑90 (MQ=255)
                                                                           atGGCAGGGAGAACTCAAAGAAGGTGGGGGTGGGGTACTTGGCGGTCGGGGGCAAAAAAACAAAGCGGCGTTTGTCACCCACCGGCAAcc          >  2:406156/1‑89 (MQ=255)
                                                                             ggCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCCCCCACCGGCAACTGc        <  1:331717/90‑1 (MQ=255)
                                                                                aGGGAGAAGTCAAAGCAGTTGGGGGTGGGGTACTTGGGGTTCGGGGTCAAAAAAAAAGAGGGGCGTTTTTCACCCACCGGCAACGGCgaa     >  2:116509/1‑90 (MQ=255)
                                                                                aGGGAGAACTCAAAGTAGTTGGTGGTGGGGTACTTGGCGTTCGGGGTCAAGAAAACAGAGCGACGTTTGTCACCCACCGGCAACTGGaaa     >  2:678960/1‑90 (MQ=255)
                                                                                aGGGAGAACTCAAAGGAGTTGGGGGTGGGGGACTTGGCGTTCGGGGTCAAGAAACCAGAGCGGCGTTTGTCACCCACCGGCAACTGCaaa     >  2:470890/1‑90 (MQ=255)
                                                                                  gTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATa   <  1:648985/90‑1 (MQ=255)
                                                                                   tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCCCCGGCAACTGCAAATAg  <  2:341209/90‑1 (MQ=255)
                                                                                   tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAg  <  1:116994/90‑1 (MQ=255)
                                                                                   |                                                                                          
ACATGATGTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAG  >  NZ_CP009273/56296‑56468

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