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.4% 222.6 / ‑4.3 61intergenic (+134/‑58)kefC/folAK(+) : H(+) antiporter KefC/dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (24/36);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

GGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGA  >  NC_000913/49703‑49826
                                                              |                                                             
ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             >  1:628144/1‑65 (MQ=255)
ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             >  1:2352507/1‑65 (MQ=255)
ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             >  1:1484603/1‑65 (MQ=255)
ggAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtga                                                             >  1:287621/1‑65 (MQ=255)
  agagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgacg                                                           >  1:2938853/1‑65 (MQ=255)
  agagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgacg                                                           >  1:1233110/1‑65 (MQ=255)
    agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCACCATCATCCTCGCATCAGTTgacgac                                                         >  1:225383/1‑65 (MQ=255)
    agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTgacgac                                                         >  1:548620/1‑65 (MQ=255)
      ggCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:1504707/64‑1 (MQ=255)
      agCGTTGCCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACgg                                                       <  1:1388401/65‑1 (MQ=255)
      agCGTGGACTCGCCGGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       >  1:2296072/1‑65 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:1981195/65‑1 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:593699/65‑1 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       >  1:1045551/1‑65 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:1947758/65‑1 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       >  1:1900701/1‑65 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:2039355/65‑1 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:1837239/65‑1 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       >  1:1789909/1‑65 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       >  1:1722545/1‑65 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       >  1:2439580/1‑65 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:2789043/65‑1 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:680646/65‑1 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       <  1:1125585/65‑1 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       >  1:856019/1‑65 (MQ=255)
      agCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACgg                                                       >  1:932096/1‑65 (MQ=255)
       gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGt                                                      <  1:1652461/65‑1 (MQ=255)
       gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGt                                                      <  1:1395636/65‑1 (MQ=255)
       gCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGt                                                      <  1:311507/65‑1 (MQ=255)
            gACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACg                                                 >  1:1863989/1‑65 (MQ=255)
             aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGc                                                >  1:1250577/1‑65 (MQ=255)
             aCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCATCGCACCAGTTGACGACGGTTTACGc                                                >  1:2944825/1‑65 (MQ=255)
              cTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCt                                               >  1:612155/1‑65 (MQ=255)
              cTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCt                                               <  1:483118/65‑1 (MQ=255)
               tCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCtt                                              <  1:2934301/65‑1 (MQ=255)
               tCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCtt                                              <  1:1651066/65‑1 (MQ=255)
               tCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCtt                                              <  1:756486/65‑1 (MQ=255)
                  ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCGGTTGACGACGGTTTACGCTTTAc                                           <  1:184464/65‑1 (MQ=255)
                  ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           <  1:2982949/65‑1 (MQ=255)
                  ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           <  1:2945556/65‑1 (MQ=255)
                  ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           <  1:1493720/65‑1 (MQ=255)
                  ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           <  1:1352712/65‑1 (MQ=255)
                  ccAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTAc                                           >  1:2882522/1‑65 (MQ=255)
                    agcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGt                                         <  1:455998/65‑1 (MQ=255)
                    agcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGt                                         <  1:660498/65‑1 (MQ=255)
                    agcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGt                                         <  1:982570/65‑1 (MQ=255)
                    agcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGt                                         <  1:1741813/65‑1 (MQ=255)
                     gcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGta                                        >  1:2823971/1‑65 (MQ=255)
                       agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:1099133/1‑65 (MQ=255)
                       agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:663729/1‑65 (MQ=255)
                       agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:1965607/1‑65 (MQ=255)
                       agAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGtata                                      >  1:599716/1‑65 (MQ=255)
                          atatAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTg                                   <  1:2610730/65‑1 (MQ=255)
                               aaaTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGAc                              >  1:195201/1‑65 (MQ=255)
                                 aTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACaa                            >  1:2498722/1‑65 (MQ=255)
                                 aTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACaa                            >  1:1046929/1‑65 (MQ=255)
                                 aTTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACaa                            >  1:2057059/1‑65 (MQ=255)
                                    ttCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAAttt                         <  1:435581/65‑1 (MQ=255)
                                             atcatcCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCg                <  1:1682727/65‑1 (MQ=255)
                                                  ccTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACCATTTTTTTTATTGGGaaa           >  1:731451/1‑65 (MQ=255)
                                                  ccTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTTTCGGGaaa           >  1:1198135/1‑65 (MQ=255)
                                                  ccTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaaa           >  1:1432027/1‑65 (MQ=255)
                                                   cTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAt          <  1:2383029/65‑1 (MQ=255)
                                                    tCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtc         <  1:208965/65‑1 (MQ=255)
                                                    tCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtc         <  1:1222154/65‑1 (MQ=255)
                                                    tCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtc         <  1:764244/65‑1 (MQ=255)
                                                      gCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAAtctc       <  1:837310/65‑1 (MQ=255)
                                                         ccAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAAt    >  1:2339604/1‑65 (MQ=255)
                                                           aGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGa  <  1:2907834/65‑1 (MQ=255)
                                                           aGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGa  <  1:835776/65‑1 (MQ=255)
                                                           aGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGa  <  1:1021669/65‑1 (MQ=255)
                                                              |                                                             
GGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGA  >  NC_000913/49703‑49826

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