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,9520TC98.2% 202.8 / ‑3.7 57intergenic (‑6/+166)mdtG/lpxLmultidrug efflux MFS transporter MdtG/kdo(2)‑lipid IV(A) lauroyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (21/35);  total (22/35)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

CACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATC  >  NZ_CP009273/1110865‑1111037
                                                                                       |                                                                                     
cACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCGGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAAccc                                                                                     <  1:614690/90‑1 (MQ=255)
              cAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGcc                                                                       <  1:3000608/90‑1 (MQ=255)
               aGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTg                                                                           >  2:539467/1‑85 (MQ=255)
               aGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTg                                                                           <  1:539467/85‑1 (MQ=255)
               aGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCa                                                                      >  1:1488656/1‑90 (MQ=255)
                 ccTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAgg                                                                    <  2:2000655/90‑1 (MQ=255)
                    aGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGgcgc                                                                 <  2:2484785/90‑1 (MQ=255)
                        aGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTg                                                             <  1:1950440/90‑1 (MQ=255)
                         ggCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGa                                                            <  1:426214/90‑1 (MQ=255)
                          gCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGGGCGGTGaa                                                           >  1:1434873/1‑90 (MQ=255)
                          gCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGaa                                                           >  1:2597819/1‑90 (MQ=255)
                              cgaTCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATa                                                       >  2:1625897/1‑90 (MQ=255)
                                   aGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATaagaag                                                  >  2:2777185/1‑90 (MQ=255)
                                     gtttcgtttcCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAgaagaa                                                <  1:1266711/90‑1 (MQ=255)
                                        tcgtttcCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAg                                             <  2:558354/90‑1 (MQ=255)
                                          gtttcCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGAt                                           <  1:1349698/90‑1 (MQ=255)
                                           tttcCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATa                                          <  1:497375/90‑1 (MQ=255)
                                           tttcCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATa                                          <  1:311744/90‑1 (MQ=255)
                                            ttcCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATaa                                         <  2:753044/90‑1 (MQ=255)
                                             tcCAGTTTATAGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGc                                       <  2:2159820/90‑1 (MQ=255)
                                              ccAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGc                                       <  2:2948554/90‑1 (MQ=255)
                                              ccAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGc                                       <  1:1625897/90‑1 (MQ=255)
                                              ccAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGc                                       <  1:1583054/90‑1 (MQ=255)
                                                aGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGcac                                     <  1:236630/90‑1 (MQ=255)
                                                aGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGcac                                     <  1:2085311/90‑1 (MQ=255)
                                                aGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGcac                                     <  1:192395/90‑1 (MQ=255)
                                                  tttATAGGGGTGTCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCacac                                   <  1:1579780/90‑1 (MQ=255)
                                                           gtgtCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgc                          >  1:386127/1‑90 (MQ=255)
                                                           gtgtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTgcg                                                                         <  1:2097349/43‑1 (MQ=255)
                                                           gtgtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTgcg                                                                         >  2:2097349/1‑43 (MQ=255)
                                                           gtgtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGca                                      <  1:1943626/78‑1 (MQ=255)
                                                           gtgtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGca                                      >  2:1943626/1‑78 (MQ=255)
                                                           gtgtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgc                          >  1:1046319/1‑90 (MQ=255)
                                                           gtgtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgc                          >  1:1587121/1‑90 (MQ=255)
                                                            tgtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgcg                         >  2:655436/1‑90 (MQ=255)
                                                             gtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgcgc                        >  1:1627419/1‑90 (MQ=255)
                                                             gtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgcgc                        >  2:533072/1‑90 (MQ=255)
                                                             gtCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATgcgc                        >  2:2460409/1‑90 (MQ=255)
                                                              tCATTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGcc                       <  1:655436/90‑1 (MQ=255)
                                                                aTTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCg                     <  1:817673/90‑1 (MQ=255)
                                                                aTTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCg                     <  2:2447465/90‑1 (MQ=255)
                                                                aTTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCg                     <  2:2349693/90‑1 (MQ=255)
                                                                aTTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCg                     <  1:533072/90‑1 (MQ=255)
                                                                aTTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCg                     <  2:2046193/90‑1 (MQ=255)
                                                                aTTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCg                     <  1:2460409/90‑1 (MQ=255)
                                                                aTTTTCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCg                     <  2:192005/90‑1 (MQ=255)
                                                                   ttCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCg                     <  1:2798835/87‑1 (MQ=255)
                                                                   ttCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCg                     >  2:2798835/1‑87 (MQ=255)
                                                                   ttCACAGGGTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACt                  <  2:2484468/90‑1 (MQ=255)
                                                                         gggTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCg                     >  1:1626583/1‑81 (MQ=255)
                                                                         gggTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCg                     <  2:1626583/81‑1 (MQ=255)
                                                                         gggTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAgg            >  2:1289590/1‑90 (MQ=255)
                                                                         gggTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAgg            >  2:2695659/1‑90 (MQ=255)
                                                                         gggCGACATAGCAACCCGCTGTTGGGGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGCCTTCCAgg            >  1:1392335/1‑90 (MQ=255)
                                                                          ggTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAggg           <  1:1289590/90‑1 (MQ=255)
                                                                          ggTGACATAGCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAggg           <  2:1460966/90‑1 (MQ=255)
                                                                                   gCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATc  >  2:659196/1‑90 (MQ=255)
                                                                                   gCAACCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATc  >  1:1925861/1‑90 (MQ=255)
                                                                                       |                                                                                     
CACCGGTAAGAAAACAGCCTAGCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATC  >  NZ_CP009273/1110865‑1111037

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