Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,102,646 T→C G140G (GGA→GGG wzx ← O16 family O‑antigen flippase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,102,6460TC96.8% 95.4 / ‑3.4 31G140G (GGA→GGGwzxO16 family O‑antigen flippase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (14/16);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

TTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTA  >  NZ_CP009273/2102561‑2102727
                                                                                     |                                                                                 
ttGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCt                                                                               >  1:521684/1‑90 (MQ=255)
  gtgtTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGa                                                                             >  1:163874/1‑90 (MQ=255)
   tgtTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGaa                                                                            <  1:190059/90‑1 (MQ=255)
     ttCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAAc                                                                          <  1:712829/90‑1 (MQ=255)
     ttCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAAc                                                                          <  2:163874/90‑1 (MQ=255)
        acaAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGc                                                                       <  2:207663/90‑1 (MQ=255)
                 aatGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAggg                                                              >  1:333607/1‑90 (MQ=255)
                  atGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGt                                                             <  2:33782/90‑1 (MQ=255)
                         agGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTaa                                                      <  1:722686/90‑1 (MQ=255)
                            gATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAAtt                                                   >  1:173025/1‑90 (MQ=255)
                                    ccAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTa                                           >  2:124492/1‑90 (MQ=255)
                                     cAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTaa                                          >  1:446291/1‑90 (MQ=255)
                                     cAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTaa                                          >  1:449490/1‑90 (MQ=255)
                                     cAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTaa                                          >  2:292914/1‑90 (MQ=255)
                                     cAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTaa                                          >  2:549280/1‑90 (MQ=255)
                                      aaGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAAt                                         >  1:63812/1‑90 (MQ=255)
                                        gCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGc                                       <  1:331312/90‑1 (MQ=255)
                                                  aTACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAAc                             <  2:278702/90‑1 (MQ=255)
                                                   tactacTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACg                            <  1:1351/90‑1 (MQ=255)
                                                   tactacTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACg                            <  2:417156/90‑1 (MQ=255)
                                                   tactacTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACg                            <  1:124492/90‑1 (MQ=255)
                                                       gctTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCa                        >  1:172799/2‑90 (MQ=255)
                                                                 aGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACg              >  1:457287/1‑90 (MQ=255)
                                                                 aGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACg              >  2:365809/1‑90 (MQ=255)
                                                                  gCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACgg             >  1:540426/1‑90 (MQ=255)
                                                                  gCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACgg             <  2:449490/90‑1 (MQ=255)
                                                                  gCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACgg             <  2:691473/90‑1 (MQ=255)
                                                                       ttcatttTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGc        <  2:716271/90‑1 (MQ=255)
                                                                        tcatttTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCt       <  1:317182/90‑1 (MQ=255)
                                                                        tcatttTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCt       <  1:176428/90‑1 (MQ=255)
                                                                          atttTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTg     >  2:267976/1‑90 (MQ=255)
                                                                             ttttctttccCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTta  <  2:157064/90‑1 (MQ=255)
                                                                             ttttctttccCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTta  <  1:479469/90‑1 (MQ=255)
                                                                                     |                                                                                 
TTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTA  >  NZ_CP009273/2102561‑2102727

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