Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 475,640 T→G D375A (GAT→GCT)  fsr ← predicted fosmidomycin efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606475,6400TG91.7% 274.2 / 7.3 72D375A (GAT→GCT) fsrpredicted fosmidomycin efflux system
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base G (27/39);  total (29/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAA  >  REL606/475566‑475712
                                                                          |                                                                        
gcaggAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATc                                                                         <  2:1534961/76‑1 (MQ=255)
 caggAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCg                                                                        <  1:731502/76‑1 (MQ=255)
 caggAATATGGTCAACATCCCCAATAGTGGCACGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCg                                                                        >  1:235604/1‑76 (MQ=255)
  aggAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCgg                                                                       <  2:184206/76‑1 (MQ=255)
  aggAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCgg                                                                       <  2:1972699/76‑1 (MQ=255)
        atGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATaa                                                                 <  1:266815/76‑1 (MQ=255)
         tGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAg                                                                <  2:251859/76‑1 (MQ=255)
         tGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAg                                                                <  1:1706017/76‑1 (MQ=255)
         tGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAg                                                                <  2:366499/76‑1 (MQ=255)
          ggTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATaaac                                                               >  1:493675/1‑74 (MQ=255)
          ggTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGGGTGGGCGGGGATAAGc                                                               >  1:1457082/1‑76 (MQ=255)
             cAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGcccc                                                            <  2:542165/76‑1 (MQ=255)
              aaCATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCa                                                           <  2:411308/76‑1 (MQ=255)
              aaCATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCa                                                           <  2:933782/76‑1 (MQ=255)
               aCATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAg                                                          <  1:253482/76‑1 (MQ=255)
               aCATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAg                                                          <  1:1637120/76‑1 (MQ=255)
                  tCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAAc                                                       <  2:998962/76‑1 (MQ=255)
                  tCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAAc                                                       <  2:998952/76‑1 (MQ=255)
                    cccAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTg                                                     <  1:1569630/76‑1 (MQ=255)
                      cAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGcc                                                   <  2:79763/76‑1 (MQ=255)
                      cAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGGGTGATCGGCGATAAGCCCCAGAACTGcc                                                   >  1:92994/1‑76 (MQ=255)
                       aaTAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCg                                                  <  1:249984/76‑1 (MQ=255)
                       aaTAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCg                                                  <  2:1742829/76‑1 (MQ=255)
                            tGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTccc                                             <  2:346306/76‑1 (MQ=255)
                             ggCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCa                                            <  1:22952/76‑1 (MQ=255)
                               cAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGa                                          <  1:1267769/76‑1 (MQ=255)
                               cAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGa                                          >  2:916707/1‑76 (MQ=255)
                               cAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGa                                          >  1:70833/1‑76 (MQ=255)
                                aGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGAc                                         <  2:19838/76‑1 (MQ=255)
                                aGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGAc                                         >  2:762119/1‑76 (MQ=255)
                                aGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGAc                                         <  2:1202791/76‑1 (MQ=255)
                                aGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGAc                                         >  1:86335/1‑76 (MQ=255)
                                  gAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCt                                       <  2:104136/76‑1 (MQ=255)
                                    aaGCACAGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTcc                                     <  1:1546412/76‑1 (MQ=255)
                                       cacaGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCAt                                  <  2:1984064/76‑1 (MQ=255)
                                       cacaGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCAt                                  <  2:1958453/76‑1 (MQ=255)
                                       cacaGATTTTATAGACTAACTCGATGCTGGCGAGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCAt                                  >  2:1531285/1‑76 (MQ=255)
                                         caGATTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGc                                >  2:1136332/1‑76 (MQ=255)
                                            aTTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCaa                             >  1:275351/1‑76 (MQ=255)
                                            aTTTTATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCaa                             <  1:446534/76‑1 (MQ=255)
                                             ttttATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCaaa                            >  1:1085304/1‑76 (MQ=255)
                                             ttttATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCaaa                            >  1:1726238/1‑76 (MQ=255)
                                              tttATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCaaaa                           >  2:492256/1‑76 (MQ=255)
                                               ttATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAg                          >  2:1481221/1‑76 (MQ=255)
                                               ttATAGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAg                          >  2:1723710/1‑76 (MQ=255)
                                                tataGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGc                         <  1:336616/76‑1 (MQ=255)
                                                  taGACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCaa                       <  1:583016/76‑1 (MQ=255)
                                                   aGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCaaa                      <  1:263892/76‑1 (MQ=255)
                                                    gACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCaaaa                     <  2:1175219/76‑1 (MQ=255)
                                                    gACTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCaaaa                     <  1:1213692/76‑1 (MQ=255)
                                                      cTAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAAcc                   >  1:1746649/1‑76 (MQ=255)
                                                       tAACTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCg                  >  2:908458/1‑76 (MQ=255)
                                                         aCTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGaa                >  2:776703/1‑76 (MQ=255)
                                                         aCTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGaa                <  1:1044678/76‑1 (MQ=255)
                                                          cTCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAATCAAAACCGaaa               >  1:1137391/1‑76 (MQ=255)
                                                           tCGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGaaaa              >  1:356234/1‑76 (MQ=255)
                                                            cGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGaaaaa             <  1:825377/76‑1 (MQ=255)
                                                            cGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGaaaaa             >  1:1245658/1‑76 (MQ=255)
                                                            cGATGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCCCCCATGCCAAAGGCAAAACCGaaaaa             >  2:41129/1‑76 (MQ=255)
                                                              aTGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAaga           >  1:1527966/1‑76 (MQ=255)
                                                               tGCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAagag          >  1:1512448/1‑76 (MQ=255)
                                                                gCTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGt         >  2:1084847/1‑76 (MQ=255)
                                                                 cTGGTGTGAGCGGCGATACGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCCAAACCGACAAAGAGTc        >  2:749035/1‑76 (MQ=255)
                                                                 cTGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTc        >  2:1691617/1‑76 (MQ=255)
                                                                  tGGTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTcc       >  2:1008163/1‑76 (MQ=255)
                                                                   ggTGTGAGCGGCGATAAGCCCCGGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCa      <  2:1659607/76‑1 (MQ=255)
                                                                   ggTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCa      <  2:1812839/76‑1 (MQ=255)
                                                                   ggTGTGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCa      >  1:1220273/1‑76 (MQ=255)
                                                                    gtgtGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAg     <  1:1664854/76‑1 (MQ=255)
                                                                    gtgtGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAg     <  2:795813/76‑1 (MQ=255)
                                                                    gtgtGAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAg     <  2:1121809/76‑1 (MQ=255)
                                                                      gtgAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGaa   >  1:816348/1‑76 (MQ=255)
                                                                      gtgAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGaa   <  2:131597/76‑1 (MQ=255)
                                                                       tgAGCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGaaa  <  1:113303/76‑1 (MQ=255)
                                                                          |                                                                        
GCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAA  >  REL606/475566‑475712

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