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,7360AT95.7% 172.3 / ‑0.7 47intergenic (+105/‑87)kefC/folAglutathione‑regulated potassium‑efflux system protein KefC/type 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (26/19);  total (28/19)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

CACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGG  >  NZ_CP009273/49651‑49813
                                                                                     |                                                                             
cacTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATtttt                                                                           <  1:575060/90‑1 (MQ=255)
cacTTTACAGCTAACTGTTTGGTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATtttt                                                                           >  1:14712/1‑90 (MQ=255)
    ttACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTc                                                                       <  1:942748/90‑1 (MQ=255)
    ttACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTc                                                                       >  2:839259/1‑90 (MQ=255)
          cTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAAcatc                                                                 <  1:1422201/90‑1 (MQ=255)
            aaCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAAcatcat                                                               >  1:765957/1‑90 (MQ=255)
            aaCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAAcatcat                                                               >  1:1563392/1‑90 (MQ=255)
            aaCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAAcatcat                                                               >  1:452698/1‑90 (MQ=255)
            aaCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAAcatc                                                                 >  2:255967/1‑88 (MQ=255)
            aaCTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAAcatc                                                                 <  1:255967/88‑1 (MQ=255)
                   tgtttttgtttCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCa                                                        >  2:885057/1‑90 (MQ=255)
                     tttttgtttCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCAcc                                                      <  1:1007363/90‑1 (MQ=255)
                              cATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTcgacga                                             <  1:563965/90‑1 (MQ=255)
                              cATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTcgacga                                             <  1:110416/90‑1 (MQ=255)
                               aTTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCgacgac                                            <  2:452698/90‑1 (MQ=255)
                               aTTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCgacgac                                            <  2:1563392/90‑1 (MQ=255)
                                    aaTGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGttt                                       <  1:1093735/90‑1 (MQ=255)
                                     aTGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCCCCAGTCGACGACGGTTTa                                      <  1:1211735/90‑1 (MQ=255)
                                     aTGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTa                                      <  2:565762/90‑1 (MQ=255)
                                       gcggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACg                                    >  1:1033971/1‑90 (MQ=255)
                                        cggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGc                                   >  1:718311/1‑90 (MQ=255)
                                         ggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCt                                  <  2:731005/90‑1 (MQ=255)
                                            gAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTa                               <  2:149227/90‑1 (MQ=255)
                                            gAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTa                               >  2:60709/1‑90 (MQ=255)
                                               tCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGt                            <  2:140005/90‑1 (MQ=255)
                                                   gggAGAGAGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAg                        <  2:153241/90‑1 (MQ=255)
                                                     gagagaGCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTg                      <  1:255059/90‑1 (MQ=255)
                                                       gagagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTcgacga                                             >  1:1031321/1‑65 (MQ=255)
                                                       gagagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTcgacga                                             <  2:1031321/65‑1 (MQ=255)
                                                       gagagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGc                    >  1:1573698/1‑90 (MQ=255)
                                                        agagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACg                                    <  2:665495/73‑1 (MQ=255)
                                                        agagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACg                                    >  1:665495/1‑73 (MQ=255)
                                                        agagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCg                   >  2:1147036/1‑90 (MQ=255)
                                                        agagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCg                   >  2:81160/1‑90 (MQ=255)
                                                        agagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCg                   >  2:166818/1‑90 (MQ=255)
                                                        agagCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCg                   >  1:203595/1‑90 (MQ=255)
                                                        agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCg                   >  1:919665/1‑90 (MQ=255)
                                                        agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCg                   >  2:1183088/1‑90 (MQ=255)
                                                           gCGTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACa                >  2:397047/1‑90 (MQ=255)
                                                             gTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGGCGGTTTACGCTTTACGTATAGTGGCGACAAt              >  1:674502/1‑90 (MQ=255)
                                                             gTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAAt              >  1:796099/1‑90 (MQ=255)
                                                             gTGGACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAAt              >  2:89227/1‑90 (MQ=255)
                                                               ggACTCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACCACGGTTTACGCTTTACGTATAGTGGCGACAAttt            >  2:103144/1‑90 (MQ=255)
                                                                   tCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAAttttttt        >  2:306506/1‑90 (MQ=255)
                                                                   tCGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAAttttttt        >  2:185663/1‑90 (MQ=255)
                                                                    cGCCAGCAGAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAAtttttttt       >  1:1001311/1‑90 (MQ=255)
                                                                         gcagAATATAAATTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCgg  >  1:85366/1‑90 (MQ=255)
                                                                                     |                                                                             
CACTTTACAGCTAACTGTTTGTTTTTGTTTCATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGG  >  NZ_CP009273/49651‑49813

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: