Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 49,736 A→T intergenic (+105/‑87) 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,7360AT80.5% 99.1 / 19.2 41intergenic (+105/‑87)kefC/folAglutathione‑regulated potassium‑efflux system protein KefC/type 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (2/6);  new base T (22/11);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 4.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAT  >  NZ_CP009273/49651‑49818
                                                                                     |                                                                                  
cacTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATtttt                                                                                >  2:135023/1‑90 (MQ=255)
    ttACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTc                                                                            >  2:195246/1‑90 (MQ=255)
    ttACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTc                                                                            >  2:167616/1‑90 (MQ=255)
        aGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAAca                                                                        <  1:227408/90‑1 (MQ=255)
          cTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAAcatc                                                                      <  2:50774/90‑1 (MQ=255)
          cTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAAcatc                                                                      <  2:137205/90‑1 (MQ=255)
          cTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAAcatc                                                                      <  1:397395/90‑1 (MQ=255)
            aaCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAAcatcat                                                                    >  1:774086/1‑90 (MQ=255)
            aaCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAAcatcat                                                                    >  1:827881/1‑90 (MQ=255)
                  ttgtttttgtttCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGc                                                              >  1:448350/1‑90 (MQ=255)
                     tttttgtttCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCAcc                                                           <  2:774086/90‑1 (MQ=255)
                              cATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTcgacga                                                  <  1:7917/90‑1 (MQ=255)
                               aTTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCgacgac                                                 >  1:674925/1‑90 (MQ=255)
                                       gcggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACg                                         <  1:282341/90‑1 (MQ=255)
                                       gcggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACg                                         <  1:237363/90‑1 (MQ=255)
                                         ggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCt                                       >  2:383050/1‑90 (MQ=255)
                                          gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCtt                                      <  2:76736/90‑1 (MQ=255)
                                          gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCtt                                      <  1:13116/90‑1 (MQ=255)
                                          gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCtt                                      <  2:81120/90‑1 (MQ=255)
                                           cgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCttt                                     <  1:270543/90‑1 (MQ=255)
                                            gAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTa                                    <  1:808711/90‑1 (MQ=255)
                                                        agagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTAc                                          <  1:573552/72‑1 (MQ=255)
                                                        agagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTAc                                          >  2:573552/1‑72 (MQ=255)
                                                        agagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCg                        >  1:637016/1‑90 (MQ=255)
                                                        agagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCg                        >  1:533104/1‑90 (MQ=255)
                                                        agagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCg                        >  1:306363/1‑90 (MQ=255)
                                                         gagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGa                       >  1:450019/1‑90 (MQ=255)
                                                         gagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGa                       >  1:400521/1‑90 (MQ=255)
                                                             gTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAAt                   >  1:761807/1‑90 (MQ=255)
                                                             gTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAAt                   >  2:63846/1‑90 (MQ=255)
                                                              tGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAAtt                  >  1:156189/1‑90 (MQ=255)
                                                               ggACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAAttt                 >  2:104234/1‑90 (MQ=255)
                                                                   tCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAAttttttt             >  2:355360/1‑90 (MQ=255)
                                                                   tCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAAttttttt             >  2:458436/1‑90 (MQ=255)
                                                                   tCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAAttttttt             >  1:529506/1‑90 (MQ=255)
                                                                   tCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAAttttttt             >  2:7496/1‑90 (MQ=255)
                                                                         gcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCgg       >  1:132352/1‑90 (MQ=255)
                                                                            gAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaa    <  2:544084/90‑1 (MQ=255)
                                                                            gAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaa    <  1:167616/90‑1 (MQ=255)
                                                                              atatAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAt  >  1:669630/1‑90 (MQ=255)
                                                                              atatAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAt  <  2:183830/90‑1 (MQ=255)
                                                                                     |                                                                                  
CACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAT  >  NZ_CP009273/49651‑49818

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: