Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 49,765 C→T intergenic (+134/‑58) kefC → / → folA glutathione‑regulated potassium‑efflux system protein KefC/type 3 dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,7650CT96.8% 105.5 / ‑4.3 31intergenic (+134/‑58)kefC/folAglutathione‑regulated potassium‑efflux system protein KefC/type 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (13/17);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

CATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGT  >  NZ_CP009273/49681‑49851
                                                                                    |                                                                                      
cATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgacga                                                                                   <  1:604674/90‑1 (MQ=255)
 aTTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTgacgac                                                                                  <  2:318897/90‑1 (MQ=255)
      aaTGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGttt                                                                             <  1:331591/90‑1 (MQ=255)
       aTGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTa                                                                            <  1:118444/90‑1 (MQ=255)
        tGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTAc                                                                           <  1:492583/90‑1 (MQ=255)
           ggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCt                                                                        <  1:351369/90‑1 (MQ=255)
            gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCtt                                                                       >  2:435388/1‑90 (MQ=255)
              gAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgacga                                                                                   >  2:595335/1‑76 (MQ=255)
              gAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgacga                                                                                   <  1:595335/76‑1 (MQ=255)
              gAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTa                                                                     <  1:292968/90‑1 (MQ=255)
              gAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTa                                                                     <  1:531962/90‑1 (MQ=255)
                      ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGt                                                             >  1:576302/1‑90 (MQ=255)
                         gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGc                                                          >  2:543054/1‑90 (MQ=255)
                          agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCg                                                         >  1:505599/1‑90 (MQ=255)
                           gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGa                                                        <  1:49649/90‑1 (MQ=255)
                             gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgacga                                                                                   >  1:442709/1‑61 (MQ=255)
                             gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgacga                                                                                   <  2:442709/61‑1 (MQ=255)
                             gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACa                                                      >  2:465164/1‑90 (MQ=255)
                              cGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACaa                                                     <  1:435388/90‑1 (MQ=255)
                                 ggACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAAttt                                                  >  2:165644/1‑90 (MQ=255)
                                     tCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAAttttttt                                              >  1:368510/1‑90 (MQ=255)
                                     tCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAAttttttt                                              >  1:136884/1‑90 (MQ=255)
                                      cGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAAtttttttt                                             >  1:413776/1‑90 (MQ=255)
                                             agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGa                                      <  1:543054/90‑1 (MQ=255)
                                              gAATATAAAATTTTCCTCAACATCATCCTCGCAACAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaa                                     <  1:345677/90‑1 (MQ=255)
                                                atatAAAATTTTCCTCAACATCTTCCTCGCAACTATTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAt                                   <  2:241148/90‑1 (MQ=255)
                                                atatAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAt                                   <  2:136884/90‑1 (MQ=255)
                                                   taAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtctc                                <  1:407166/90‑1 (MQ=255)
                                                     aaaTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCaa                              <  1:516827/90‑1 (MQ=255)
                                                        ttttCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGa                           >  1:40474/1‑90 (MQ=255)
                                                          ttCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc                         <  2:80741/90‑1 (MQ=255)
                                                                      atcCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcg             >  1:419261/1‑90 (MQ=255)
                                                                      atcCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGGGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcg             >  1:347072/1‑90 (MQ=255)
                                                                           cGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcgg            >  2:342939/1‑86 (MQ=255)
                                                                           cGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTgcgg            <  1:342939/86‑1 (MQ=255)
                                                                                 aGTTGACGACGGGTTACGCTTTACGTATAGTGGCGACAATTTGCTTGGTCGGGAAATGTCAATGAGCAGTCTGGTTGCGGCGTTAGCgga  >  2:508255/1‑89 (MQ=255)
                                                                                    |                                                                                      
CATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGT  >  NZ_CP009273/49681‑49851

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