Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 879,010 T→C intergenic (‑166/‑119) ybjG ← / → mdfA undecaprenyl‑diphosphate phosphatase/multidrug efflux MFS transporter MdfA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273879,0100TC84.4% 80.5 / 2.7 32intergenic (‑166/‑119)ybjG/mdfAundecaprenyl‑diphosphate phosphatase/multidrug efflux MFS transporter MdfA
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base C (14/13);  total (17/15)
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.00e+00

AAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTA  >  NZ_CP009273/878923‑879094
                                                                                       |                                                                                     
aaGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGCAg                                                                                     >  1:29465/1‑90 (MQ=255)
       aTCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGCAGCTATACt                                                                              <  1:579559/90‑1 (MQ=255)
          gTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGCAGCTATACTCGt                                                                           <  1:805160/90‑1 (MQ=255)
            cTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGCAGCta                                                                                  >  2:227979/1‑81 (MQ=255)
            cTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGCAGCta                                                                                  <  1:227979/81‑1 (MQ=255)
                  tCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGCAGCTATACTCGTAATAATGt                                                                  <  2:610549/90‑1 (MQ=255)
                   cccTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGCAGCTATACTCGTAATAATGTa                                                                 >  2:97671/1‑90 (MQ=255)
                    ccTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGCAGCTATACTCGtaat                                                                        <  2:189893/83‑1 (MQ=255)
                    ccTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGCAGCTATACTCGtaat                                                                        >  1:189893/1‑83 (MQ=255)
                                         tgggtgATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGttt                                           >  1:51012/1‑90 (MQ=255)
                                           ggtgATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCa                                         >  2:368783/1‑90 (MQ=255)
                                                  gTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACaa                                  >  1:404044/1‑90 (MQ=255)
                                                             aGGTTAAATAAAACTTAAAGAAAGCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTtc                       <  2:29465/90‑1 (MQ=255)
                                                              ggTTAAATAAAACTTAAAGAAAGCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTtct                      >  2:516973/1‑90 (MQ=255)
                                                                 taaataaaACTTAAAGAAAGCGTAGCTATACTCGTAATAATTGTAAGaa                                                             >  2:912425/1‑49 (MQ=39)
                                                                 taaataaaACTTAAAGAAAGCGTAGCTATACTCGTAATAATTGTAAGaa                                                             <  1:912425/49‑1 (MQ=39)
                                                                 taaataaaACTTAAAGAAAGCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGt                   <  2:250070/90‑1 (MQ=255)
                                                                  aaataaaACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGtt                  >  1:338731/1‑90 (MQ=255)
                                                                  aaataaaACTTAAAGAAAGCGCAGCTATACTCGTAATAATGTAAGaa                                                             <  2:210696/47‑1 (MQ=255)
                                                                  aaataaaACTTAAAGAAAGCGCAGCTATACTCGTAATAATGTAAGaa                                                             >  1:210696/1‑47 (MQ=255)
                                                                  aaataaaACTTAAAGAAAGCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGtt                  >  2:58499/1‑90 (MQ=255)
                                                                        aaattAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGTTACAAAATTCGCTTTCTCGTTAGCTgc            >  1:97381/4‑90 (MQ=255)
                                                                        aaCTTAAAGAAAGCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTgc            >  1:704474/1‑90 (MQ=255)
                                                                         aCTTAAAGAAAGCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTgcg           >  2:632359/1‑90 (MQ=255)
                                                                           ttAAAGAAAGCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCt         <  2:218794/90‑1 (MQ=255)
                                                                           ttAAAGAAAGCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCt         <  2:659020/90‑1 (MQ=255)
                                                                             aaagaaagCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCttt       <  2:266874/90‑1 (MQ=255)
                                                                             aaagaaagCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCttt       <  2:509264/90‑1 (MQ=255)
                                                                             aaagaaagCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCttt       <  1:632359/90‑1 (MQ=255)
                                                                             aaagaaagCCTAGCTATACTCGTATTAATGTAAGAATGTGCTTAAGCGTGGTTTCAGCTACAAAGTTCGCTTTCTAGTTAGCTGCGCttt       <  2:97381/90‑1 (MQ=255)
                                                                                  aagCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTttatta  >  1:812154/1‑90 (MQ=255)
                                                                                  aagCGCAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTttatta  >  1:462681/1‑90 (MQ=255)
                                                                                       |                                                                                     
AAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTA  >  NZ_CP009273/878923‑879094

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