Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927356,3790TG57.6% 25.4 / 42.0 33Y244S (TAC→TCC) lptDLPS assembly protein LptD
Reads supporting (aligned to +/- strand):  ref base T (1/13);  new base G (19/0);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 2.44e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-07
Rejected as consensus: Frequency below/above cutoff threshold.

CATGATGTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAG  >  NZ_CP009273/56297‑56468
                                                                                  |                                                                                         
caTGATGTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTc                                                                                    <  2:801786/90‑1 (MQ=255)
 atgatgTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCa                                                                                   <  2:359458/90‑1 (MQ=255)
     tgTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAagt                                                                               <  1:857206/90‑1 (MQ=255)
     tgTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGGACGCaaaat                                                                               >  2:274854/1‑88 (MQ=255)
      gTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAagta                                                                              >  1:610388/1‑90 (MQ=255)
        tGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAagtagt                                                                            >  1:68940/1‑90 (MQ=255)
         gCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGtt                                                                           <  2:185841/90‑1 (MQ=255)
              gacgaTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTAGTtggtgg                                                                      >  1:477099/1‑90 (MQ=255)
               acgaTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTtggtggt                                                                     <  1:908003/90‑1 (MQ=255)
                  aTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTggtggtggt                                                                  <  1:256564/90‑1 (MQ=255)
                  aTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTAGTTggtgggggt                                                                  >  1:581072/1‑90 (MQ=255)
                  aTGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAAACAAAGTAGTtggggggggg                                                                  >  2:374972/1‑89 (MQ=255)
                  aTGCATATAATGCGGCGTGATGGTGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAAGAGTtgggggggggg                                                                 >  1:542377/1‑88 (MQ=255)
                          aaTGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAAGAGTTGGGGGGGGTGTACTTgg                                                          >  2:202274/1‑90 (MQ=255)
                               ggcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGTGGTGGTGTACTTGGCGTTCg                                                    <  1:508441/90‑1 (MQ=255)
                                gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTAGTTTGGGGTGGTGGACTTTGCGTTCg                                                    >  1:29986/1‑90 (MQ=255)
                                gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTAGTTGGGGGTGGTGTACTTGGCGTTCg                                                    >  1:326195/1‑90 (MQ=255)
                                gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAATAGTTGGGGGTGGGGTACTTGGCGGTCg                                                    >  1:831301/1‑90 (MQ=255)
                                gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAATAGTTGGGGGGGGTGTACTTGGCGTTCg                                                    >  2:319314/1‑90 (MQ=255)
                                gcgTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAATAGTTGGGGGGGGGGTACTTGGGGTTCg                                                    >  2:242873/1‑90 (MQ=255)
                                   tGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGa                                                 <  1:202274/90‑1 (MQ=255)
                                            cATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAAc                                        <  2:2669/90‑1 (MQ=255)
                                              tCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTCGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCa                                      <  2:899555/90‑1 (MQ=255)
                                               ccATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAGTAGTTGGTGGGGGGGTACTTGGCGTTCGGGGTCAAGAAACCag                                     >  1:514457/1‑90 (MQ=255)
                                               ccATATTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAAATCAAAATAGTTGGGGGGGGTGTACTTGGCGGTCGGGGTTAAGAAACCag                                     >  2:177666/1‑90 (MQ=255)
                                                  tatTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGc                                  <  2:386898/90‑1 (MQ=255)
                                                  tatTTGGCGCGATGTTCCAGTAATATGGCAGGGAGAACTCAAAATAGTTGGGGGGGGGGTACTTTGCGTTCGGGATCAAGAAACCAGAGc                                  >  2:789283/1‑90 (MQ=255)
                                                        gcgcGATGTTCCAGTAATATGGCAGGGAGAACTCAAAAAAGTTGGTGGGGGGGTACTTGGCGGTCGGGGTCAAGAAACCAGAGCGACGtt                            >  2:87723/1‑90 (MQ=255)
                                                          gcgATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGGTCAAGAAACCAGAGCGACGTTTg                          >  1:910777/1‑90 (MQ=255)
                                                             aTGTTCCAGTAATATGGCAGGGAGAACTCAAAGGAGTTGGTGGGGGGGTACTTGGGGTTCGGGGTCAAGAAAACAGAGCGGCGGTTGTca                       >  2:27700/1‑90 (MQ=255)
                                                                        atatGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCaa            <  1:637166/90‑1 (MQ=255)
                                                                            ggCAGGTAGAACTCAAAGTAGTTGGTTGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGc        <  2:406063/90‑1 (MQ=255)
                                                                            ggCAGGGAGAACTCAAAATAGTTGGTGGGGGGGGACTTGGCGTTCGGGGTCAAGAAAACAAAGCGACGTTTGTCACCCACCGGCCACTGc        >  2:342676/1‑90 (MQ=255)
                                                                                  tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAg  <  1:719443/90‑1 (MQ=255)
                                                                                  tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAg  <  1:431803/90‑1 (MQ=255)
                                                                                  |                                                                                         
CATGATGTTGCCACGACGATGCATATAATGCGGCGTGATGGTGGCATCCATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAATAG  >  NZ_CP009273/56297‑56468

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