Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,110,952 T→C intergenic (‑6/+166) mdtG ← / ← lpxL multidrug efflux MFS transporter MdtG/kdo(2)‑lipid IV(A) lauroyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,110,9520TC92.9% 74.4 / ‑0.4 28intergenic (‑6/+166)mdtG/lpxLmultidrug efflux MFS transporter MdtG/kdo(2)‑lipid IV(A) lauroyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (16/10);  total (17/11)
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.38e-01

AGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATC  >  NZ_CP009273/1110880‑1111037
                                                                        |                                                                                     
aGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCa                                                                      >  1:98456/1‑90 (MQ=255)
        cAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGt                                                              <  1:363344/90‑1 (MQ=255)
          ggCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGa                                                            >  2:801466/1‑90 (MQ=255)
           gCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGaa                                                           >  1:256540/1‑90 (MQ=255)
           gCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGCGTCATTTTCCCAGGGAGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGaa                                                           >  1:1523410/1‑90 (MQ=255)
                       tttcgtttcCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGaaa                                               <  1:878864/90‑1 (MQ=255)
                       tttcgtttcCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGaaa                                               <  1:896046/90‑1 (MQ=255)
                            tttcCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATa                                          <  1:21412/90‑1 (MQ=255)
                               ccAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGc                                       >  2:265433/1‑90 (MQ=255)
                                 aGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGcac                                     <  2:728510/90‑1 (MQ=255)
                                 aGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGcac                                     <  2:1450781/90‑1 (MQ=255)
                                 aGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGcac                                     <  1:833795/90‑1 (MQ=255)
                                   tttATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCacac                                   <  1:1027260/90‑1 (MQ=255)
                                            gtgtCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTg                                                                           <  1:926476/41‑1 (MQ=255)
                                            gtgtCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTg                                                                           >  2:926476/1‑41 (MQ=255)
                                            gtgtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgc                          >  2:1392772/1‑90 (MQ=255)
                                              gtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgcgc                        >  1:587778/1‑90 (MQ=255)
                                              gtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgcgc                        >  1:1276019/1‑90 (MQ=255)
                                              gtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgcgc                        >  1:1220521/1‑90 (MQ=255)
                                                cATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGccc                      >  2:225447/1‑90 (MQ=255)
                                                      cacaGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTc                <  1:555917/90‑1 (MQ=255)
                                                        caGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCa              >  1:876872/1‑90 (MQ=255)
                                                          gggTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAgg            >  2:310536/1‑90 (MQ=255)
                                                                    gCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCg                     >  2:128929/1‑71 (MQ=255)
                                                                    gCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCg                     <  1:128929/71‑1 (MQ=255)
                                                                    gCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGGCTTCCAGGGGGCGCAATc  >  2:911885/1‑90 (MQ=255)
                                                                    gCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATc  >  2:419375/1‑90 (MQ=255)
                                                                    gCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATc  >  2:498256/1‑90 (MQ=255)
                                                                        |                                                                                     
AGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATC  >  NZ_CP009273/1110880‑1111037

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