Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 49,765 C→T intergenic (+134/‑58) kefC → / → folA K(+) : H(+) antiporter KefC/dihydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091349,7650CT98.0% 179.2 / ‑5.3 50intergenic (+134/‑58)kefC/folAK(+) : H(+) antiporter KefC/dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (26/23);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGA  >  NC_000913/49703‑49826
                                                              |                                                             
ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             >  1:924134/1‑65 (MQ=255)
ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             <  1:2465180/65‑1 (MQ=255)
ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             >  1:2330130/1‑65 (MQ=255)
ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             >  1:3403084/1‑65 (MQ=255)
 gagagaGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgac                                                            >  1:1107266/1‑65 (MQ=255)
   gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgacga                                                          >  1:3445662/1‑65 (MQ=255)
    agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTgacgac                                                         >  1:2518516/1‑65 (MQ=255)
    agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTgacgac                                                         >  1:1826854/1‑65 (MQ=255)
      agCGTGGACTCGTCAGCAGAATATAAAATCTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:2476109/65‑1 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCGCCAGTTGACGACgg                                                       <  1:852665/65‑1 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       >  1:2133899/1‑65 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:3568109/65‑1 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:2062600/65‑1 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:876923/65‑1 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:1586325/65‑1 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       >  1:905725/1‑65 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAAGATTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:2453522/65‑1 (MQ=255)
       gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGANGACGGt                                                      <  1:2222536/65‑1 (MQ=255)
        cGTGGACTCGCCAGCAGACTATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGtt                                                     >  1:3486118/1‑65 (MQ=255)
           ggACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTAc                                                  >  1:924289/1‑65 (MQ=255)
            gACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACg                                                 <  1:2945424/65‑1 (MQ=255)
             aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                >  1:513921/1‑65 (MQ=255)
             aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                >  1:71573/1‑65 (MQ=255)
             aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                >  1:2262856/1‑65 (MQ=255)
             aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                >  1:1865004/1‑65 (MQ=255)
             aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                >  1:3203408/1‑65 (MQ=255)
             aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                >  1:1451633/1‑65 (MQ=255)
              cTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCt                                               <  1:2762054/65‑1 (MQ=255)
              cTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCt                                               <  1:136489/65‑1 (MQ=255)
              cTCGCCAGCAGAATATAAAATTGTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCt                                               >  1:2442179/1‑65 (MQ=255)
                 gCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTa                                            <  1:3715839/65‑1 (MQ=255)
                  ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGTTTTAc                                           <  1:2187515/65‑1 (MQ=255)
                  ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           <  1:2743213/65‑1 (MQ=255)
                  ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           <  1:573473/65‑1 (MQ=255)
                  ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           <  1:306120/65‑1 (MQ=255)
                  ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           >  1:3029424/1‑65 (MQ=255)
                      cagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtat                                       <  1:2886375/65‑1 (MQ=255)
                       agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:3137606/1‑65 (MQ=255)
                       agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:2318376/1‑65 (MQ=255)
                       agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:1800530/1‑65 (MQ=255)
                       agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:2173054/1‑65 (MQ=255)
                       agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:3647676/1‑65 (MQ=255)
                             taAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCg                                <  1:772811/65‑1 (MQ=255)
                               aaaTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGAc                              >  1:2941005/1‑65 (MQ=255)
                                 aTTTTCCTCAACATCATCCTCGCACCAGTTGGCGACGGTTTACGCTTTACGTATAGTGGCGACaa                            >  1:1963593/1‑65 (MQ=255)
                                 aTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACaa                            >  1:3470065/1‑65 (MQ=255)
                                 aTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGAGGCGaaaa                            >  1:211904/1‑65 (MQ=255)
                                 aTTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACaa                            >  1:2111297/1‑65 (MQ=255)
                                                  ccTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaaa           >  1:2128020/1‑65 (MQ=255)
                                                  ccTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaaa           >  1:2361068/1‑65 (MQ=255)
                                                   cTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAt          <  1:861640/65‑1 (MQ=255)
                                                    tCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtc         <  1:2713087/65‑1 (MQ=255)
                                                    tCGCACCAGTTGACGACGGGTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtc         <  1:1175250/65‑1 (MQ=255)
                                                      gCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtctc       <  1:513326/65‑1 (MQ=255)
                                                      gCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtctc       <  1:1548249/65‑1 (MQ=255)
                                                        aCCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTTaa     >  1:652439/1‑65 (MQ=255)
                                                        aCCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCaa     <  1:3475203/65‑1 (MQ=255)
                                                        aCCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCaa     >  1:2995550/1‑65 (MQ=255)
                                                           aGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGa  >  1:50520/1‑65 (MQ=255)
                                                              |                                                             
GGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGA  >  NC_000913/49703‑49826

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: