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,9520TC93.8% 41.4 / ‑3.5 16intergenic (‑6/+166)mdtG/lpxLmultidrug efflux MFS transporter MdtG/kdo(2)‑lipid IV(A) lauroyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (6/9);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

GGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATC  >  NZ_CP009273/1110869‑1111037
                                                                                   |                                                                                     
ggTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTg                                                                                 <  2:423573/90‑1 (MQ=255)
                          cgaTCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATa                                                       <  1:11232/90‑1 (MQ=255)
                                        ttcCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATaa                                         <  1:149865/90‑1 (MQ=255)
                                          ccAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGc                                       <  1:351585/90‑1 (MQ=255)
                                          ccAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGc                                       <  1:59790/90‑1 (MQ=255)
                                            aGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGcac                                     <  2:134558/90‑1 (MQ=255)
                                              tttATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCacac                                   <  1:60441/90‑1 (MQ=255)
                                                       gtgtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATa                                                       <  1:105311/61‑1 (MQ=255)
                                                       gtgtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATa                                                       >  2:105311/1‑61 (MQ=255)
                                                       gtgtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgc                          >  2:8677/1‑90 (MQ=255)
                                                         gtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgcgc                        >  1:253073/1‑90 (MQ=255)
                                                                      ggTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAggg           <  2:253073/90‑1 (MQ=255)
                                                                               gCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATa                                          >  1:475790/1‑50 (MQ=255)
                                                                               gCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATa                                          <  2:475790/50‑1 (MQ=255)
                                                                               gCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATc  >  1:225096/1‑90 (MQ=255)
                                                                               gCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATc  >  2:147619/1‑90 (MQ=255)
                                                                                   |                                                                                     
GGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATC  >  NZ_CP009273/1110869‑1111037

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: